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SYNOPSIS 


INVESTIGATION INTO THE DESIGN OF POWER SYSTEM ST/aBILIZERS 

With increasing size and complexity of modern power 
systems, stability investigations are assuming greater signi- 
ficance in system planning, design and operation. Transient 
stability of power systems is improved by incorporating fast 
acting and high gain excitation systems with synchronous 
generators. However, such excitation systems introduce nega- 
tive damping in the system impairing dynamic stability. 
Oscillations of small magnitude and low frequency, in the 
range of 0.2 to 2.0 Hz, can persist for long periods of time 
and in some cases present limitations to the power transfer 
capability. Power System Stabilizers (PSS) are developed to 
damp-out these oscillations. PSS are dynamic compensators 
which receive a feedback signal from rotor speed, angular 
position, frequency or electrical power and provide a corre- 
ctive input to the excitation system. 

The basis for the design of PSS is explained in the 
classic paper by deMello and Concordia [l]. The PSS is 
designed as a dynamic conpensator, using input derived from 



local control signal, with the objective of increasing the 
damping torque v^ile maintaining the synchronizing torque. 

The model given in [l] is a simplified one w^ere the power 
system is represented by a single generator feeding an 
infinite bus through a transmission line. The synchronous 
machine is represented by a third order model neglecting 
damper -windings. This simplified model of the power system 
permits application of classical control theory for the 
design of PSS. 

Since then, there has been much work on power system 
stabilization reported in the literature using classical 
and modern control theory. In the latter approach the PSS 
design is based on the objective of assignment of closed- 
loop poles or minimizing a performance index. Inspite of 
the extensive literature on PSS there are still some gaps 
which are described below, t 

(1) It is observed that a simple model of a power system 
(that of single machine connected to infinite bus) is not 
adequate in many situations and coordinated application of PSS 
using multimachine system model is essential [2], The design 
based on the concept of damping torque, utilizing classical 
control theory is not adequate here. However the direct appli- 
cation of modern control theory is also not feasible due to the 
practical constraint of having to use only local control 



signals, some theoretical work on decentralized control 

in large systems has been carried out, its application for 
design of PSS in a large power system is not reported. 

If PSS are to be equipped on the generators that are 
already installed in a system for the purpose of improving 
dynamic stability, there arises the question of selecting 
appropriate locations for PSS in order to provide an effective 
control. This question has not also been satisfactorily 
answered. 

(2) The objective of PSS design given in [l] is to increase 
the damping torque in the generators. It can be shown that 
this is equivalent to shifting the complex pair of eigen- 
values, corresponding to electromechanical oscillations, to 
the left in the complex plane [3]. This indicates that pole 
assignment is a suitable objective for PSS design. However, 
the design freedom (in the selection of parameters) is not 
fully utilized in the work reported in this area; for example, 
with a second order dynamic compensator, using a single input, 
it is possible to assign upto five closed-loop poles to their 
desired locations* With more complex compensators (multi- 
input) it may be possible to assign more closed-loop poles. 

Generally, for the purpose of stabilization, it is ade- 
quate to specify the region in the complex-plane where closed- 
loop poles should lie, rather than specifying exact locations. 



Ho'vever, unique compensators cannot be obtained unless exact 
locations are specified. This problem can be, ^voided if a 
suitable performance index is chosen as an objed'tive in con- 
junction with the application of optimal control 

(3) While there is some work reported on the choic#'^ 
control signals [4], a detailed study comparing their i^lilttive 
merits is not available. The comparison should be based on 
satisfactory performance of the controller under varying condi- 
tions of operating point and system strength. 

Interaction between two or more oarallel connected gene- 
rators in a power station complicates the choice of control 
signals for PSS. A detailed analysis is not available in the 
literature. 

(4) Third order synchronous machine model is 'widely used for 
PSS design. Adequacy of this model needs to be tested. 

The work reported in this thesis is directed at trying to 
fill the gaps mentioned above. The major contributions are as 
follows : 

(1) Study of decentralized control for the design of PSS in 
large power systems by considering application to a particular 
large power system, 

(2) Determination of effective PSS locations in a large multi- 
machine system through eigenvalue sensitivity analysis and 
validation of the results. 



(3) A detailed investigation on the effectiveness of various 
feedback control signals (viz. rotor speed, electrical power 
and bus frequency) , for PSS design, including consideration of 
interaction between adjacent machines in the choice of signal. 

(4) Application of modern control theory for the design of 
dynamic compensators (to be used as PSS) utilizing full design 
freedom. Comparison of the performance of PSS designed using 
pole assignment techniques, parameter optimization technique 
and classical control method. 

(5) Verification of the adequacy of using simplified model of 
synchronous machine in PSS design. 

An outline of the work reported in the thesis is given 
below : 

(1) Chapter 1 introduces various aspects of dynamic stability 
and design of PSS in large multimachine power systems and re- 
views the literature in this area. 

(2) Chapter 2 deals with design techniques for PSS. The tech- 
niques developed here are general enough to be applied to multi- 
machine systems. Pole-assignment in closed-loop system is a 
well accepted design objective for PSS and is used in the first 
two techniques presented. In the first technique, modal control 
theory is utilized for PSS design. Only accessible states are 
used as feedback signals to PSS. Pole assignment technique 
given by Munro and Hirbod [5], with dynamic output feedback, is 



utilized for PSS design in the second techniques. 

Electromechanical oscillations of machine rotors with 
respect to one another cause power surges between machines 
and limit power transfer capabilities* particularly in 
systems having long transmission lines. Damping of these 
oscillations can be achieved by minimization of power flu- 
ctuations in generator output. The third design technique 
is based on parameter optimization with the objective to 
minimize power fluctuations in generator outputs due to dis- 
turbances. The three design techniques presented in this 
chapter are compared with the classical design technique used 
in the power industry. 

(3) The main study in chapter 3 is concerned with the 
adequacy of simplified machine models, neglecting the effect 
of danger circuits, for the design of PSS. This is done by 
employing the PSS designed for lower order machine model to the 
the machine represented in detail, considering the effect of 
damper circuits. Conparison between the performances of the 
stabilizers designed using simplified and detailed models of 
the machine is presented based on time-domain analysis of the 
closed-loop systems. 

(4) Qiapter 4 presents a case study on the evaluation of the 
effectiveness of various feedback signals that can be used for 
PSS design. Three control signals, namely, rotor speed. 



electrical power output and terminal bus frequency are con- 
sidered. PSS for the three control signals are designed and 
the effectiveness of signals is judged by determining stable 
regions in the P-Q domain. The influence of ac system strength 
on the design of PSS utilizing the three control signals, is 
studied. Analysis of interaction between two generating units 
equipped with PSS, in a power plant is also carried out. The 
results indicate the importance of using average rotor speed 
instead of individual speeds as control signal. 

(5) Chapter 5 deals with the development of a state-space 
model for large power systems, suitable for dynamic stability 
analysis and the selection of proper locations of PSS. The 
system model of the multimachine system is obtained systema- 
tically by developing the models of various machines and their 
interconnections through the network. The network is represent 
ed by its Jacobian matrix. 

Tvi?o analytical techniques are presented for the choice of 
effective locations of PSS. A case study of a 13-machine sys- 
tem is presented for the selection of PSS locations. Loci of 
eigenvalues are obtained to validate the results obtained 
through the analytical techniques . 

(6j Chapter 6 presents the application of two design tech- 
niques for PSS in a multimachine system using decentralized 
feedback control. The first technique is based on the 



algorithm given in [6]. PSS for the selected machines are 
designed, neglecting the interconnection between machines, to 
assign closed-loop eigenvalues. Parameters of the stabilizers 
are there updated by introducing the interconnections in dis- 
crete increments, keeping the closed-loop eigenvalues unchanged. 
The algorithm is found to be not suitable vshen the constraint 
of providing PSS on limited number of machines is considered. 

In the second technique a perfomaance index, based on mini- 
mization of fluctuations in the power outputs of relevant gene- 
rators, is used to design PSS using parameter optimization 
technique. 

Two examples of power systems, one with 3 machines and the 
other with 13 machines, are considered to illustrate the tech- 
niques . 

(7) The concluding chapter summarizes the research results and 
indicates possibilities for future work. 
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CHAPTER 1 


INTRODUCTION 


1 . 1 GENEPu^L 

One of the important considerations in the planning, 
design and operation of a power system is the stability of the 
system. Stability of a power system is usually classified as 
’transient stability' and 'dynamic stability’. The forri^r 
refers to the stability of the system under large disturbance 
such as fault on a heavily loaded line which requires opening 
of the line or the tripping of a loaded generator or an abrupt 
drop in a large load. The latter refers to the stability of 
the system under small perturbations such as random change in 
loading conditions of the system. Transient stability of a 
power system is improved by the use of high speed protection 
gear and by incorporating fast acting and high gain excitation 
systems with the synchronous machines. Although the use of 
such an excitation system improves transient stability, it 
also introduces negative damping in the system, impairing 
dynamic stability. Another reason for an increasing tendency 
tovi/ards dynamic instability is decrease in the transmission 
strength of the system relative to the size of generating 
stations. Low frequency oscillations of a sustained or growing 
nature have been reported in the literature by powex utilities 
[l“5]. These oscillations typically occur in the frequency 
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range of approximately 0.2 to 2.0 Hz and insufficient damping 
of these oscillations may limit the power transfer capability 
in the system* These oscillations in power can be related to 
oscillations of rotors of synchronous machines in the system 
relative to one another. Power system stabilizers (PSS) 
have been developed to provide damping of these electro- 
mechanical oscillations. 

Power system stabilizers are auxiliary controllers 
which receive feedback signal from rotor speed, angular posi- 
tion, electrical power output or bus frequency and provide a 
supplementary stabilizing signal to the excitation system of 
synchronous machines. 

Design of PSS has received much attention in the power 
industry. Industrial PSS design is primarily based on the 
work of de Mello and Concordia [6] . The work reported in 
[6] provided an insight into the nature of damping and syn- 
chronizing torques which is utilized in PSS design. They 
used a simplified power system model consisting of a syn- 
chronous machine connected to an infinite bus through a 
transmission line. This model enables use of classical 
control theory for the design of PSS. It should, however, be 
noted that this representation of power system is not adequate 
in general and coordinated application of PSS using multimachine 
system model is necessary [l7]. Classical control theory is not 
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adequate in such cases. Modern multivariable system control 
theory has to be used for PSS design in multimachine systems. 
Two approaches would be useful for PSS design, in multimachine 
systems, namely those based on optimal control theory and pole 
assignment . 

The main theme of this thesis is to investigate into the 
design procedures of PSS for single and multimachine systems 
using modem control theory. Before discussing the work 
reported here, a brief literature survey on the design of PSS 
will be undertaken. 

1.2 DESIO^ OF POWER SYSTBi ST.^SILIZERS : STATE OF THE ART 
1.2.1 Design Objective 

The different methods available for PSS design can be 
broadly classified as : 

( 1) Improvement of damping torque. This objective of PSS 
design is used by most of the pow/er industries. A 
single machine infinite bus equivalent of the actual 
system is considered and usually Hef fron-Phillips 
model [?] is obtained to apply the frequency domain 
classical design procedure given by de Mello 
Concordia [6]. 
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(2) Minimization of a perfoimance index. Not much work has 
been done to develop design methodology based on this 
approach and its application is not adopted by the power 
industry yet. Lack of feel for performance index is 
another reason for its poor acceptance. 

(3) Assignment of closed-loop eigenvalues. This objective 
of PSS design is slowly being adopted by power engineers 
particularly where it is difficult to represent the 
system by single-machine infinite bus model and coordi- 
nated design of PSS is essential. 

1.2.2 Power System Model 

In most of the literature on PSS design, a single-machine 
infinite-bus equivalent of power system is considered [l-7] . 

In deriving the system model, the effect of turbine and gover- 
nor dynamics is usually ignored because of the large time 
constants associated with them. Damper circuits are aleo 
neglected because their effect is considered to be small. Also, 
since damper circuits provide some damping of rotor oscilla- 
tions, neglecting them, leads to a conservative design. Another 
assumption which is usually made is that the network is in 
steady-state. Dynamical equations of machine and excitation 
system are described as a set of first order differential 
equations. .Although, the system is nonlinear, for the purpose 
of PSS design, dynamical equations are linearized around an 
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operating point, to get the system model in the linearized 
state-space form as 

P X = [a] X + [b] u (1.1) 

where X and u are state and input variables. ’A* is known 
as the coefficient or system matrix. The single-machine 
infinite-bus equivalent model can be justified when design- 
ing PSS on a generator connected to a load centre by a long 
transmission line. Even in this case it may be inaccurate to 
neglect the interaction between adjacent generators in a 
power station. Baker et al . [Sj investigated the excitation 
system interaction in a two unit plant feeding an infinite 
bus. Their investigation reveals that dynamic behaviour of 
machines in isolation is much different from that \ftiien they 
are operating in parallel. 

When a single-machine infinite-bus equivalent model is 
not adequate, it is necessary to consider a detailed system 
model including the dynamics of ail the relevant machines. 

The linearized model of a multimachine system can still be 
obtained in the general state-space form given in egn. (1.1). 
To simplify the model it may be sufficient to represent the 
generators, on which PSS is not being designed, by classical 
models, neglecting flux decay. 



The above equation is general enou^ to consider any 
degree of detail in representing the machine, exciter and 
turbine-governor systems (if desired). The excitation system 
used on modern generators have been compiled into standard 
categories [9]. 

The formulation of system matrix is rather involved and 
recently many papers have been published on this [10-16] . 
Laughton [lO] obtained the system matrix from the algebraic 
and differential equations representing the entire power 
system by employing matrix reductions and inversions. 

Van Ness and Goddard [ll] and Undrill [12] have given methods 
of building up the system matrix from submatrices representing 
individual components. Methods of forming the system matrix 
which do not involve matrix inversion are suggested in Refs. 

[ 14-16] « In all the formulations given in Refs, [10-15], the 
power system network is reduced by treating loads as constant 
impedances. Padiyar et al, [16] have suggested a systematic 
method for the formulation of system model for large AC/DC povi«r 
systems. The system model is obtained by developing the models 
of various components and subsystems (such as machines, excita- 
tion and governor systems) and their interconnection through 
the network which is represented by its Jacobian, The effect of 
any noolinfear voltage dependent loads can be obtained by mcxiify- 
ing the corresponding diagonal blocks in the Jacobian. 
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1.2.3 Selection of PS5 Location 

In a large .aultimachine system there are many modes of 
rotor oscillations, some of '/^lich can pose problem in terms 
of the dynamic stability. This can be solved by providing PSS 
on selected generators. 

The choice of generators on v^hich PSS is to be equipped 
is complicated by the fact that not all machines are effective 
in damping a particular mode of oscillation. Furthermore , PSS 
on an arbitrarily chosen machine can introduce negative damp- 
ing on some modes while damping other modes. The selection of 
PSS location can be done by trial and error, but this is very 
time consuming. Some recent papers have addressed the ques- 
tion of fast analytical meth.ods to determine the best loca- 
tion of PSS to a particular mode [ 17-20] . de Mello et al« 

[17] used eigenvectors to develop the criteria for optimum 
location. But this is not the best approach. Recently a 
technique based on eigenvalue sensitivity has been given for 
the selection of PSS location [18,19J. This is a fast tech- 
nique as it avoids simulation, at the saif^ tine it can reveal 
the conflicts in the design of PSS (e.g., whether PSS can 
undamp some modes), fieference [20] also gives a similar idea 
but the selection of PSS is based on simulation. 

It is to be noted that case studies on large systems have 
not been reported in the application of tiie techniques given 
in [ 17 - 203 . 



1.2,4 Choice of Control Signal 


As the objective of PSS design is to produce a component 
of electrical torque in phase with rotor velocity [6], it is 
natural to use rotor speed as feedback signal to PSS. Accelerat 
ing power and frequency deviation are directly related to the 
rotor speed and, therefore, these three signals are widely 
used for the design of PSS [21-27j. Conceptually other signals 
can also be used as input to PSS, Signals such as armature 
current, field current and reactive power were investigated 
[2l] and it was found that the signals derived from rotor 
speed are more effective than these signals. It was observed 
[4,22] that speed input stabilizers can contribute to negative 
damping of shaft torsional oscillation mode. The use of acce- 
lerating power as control signal has received considerable 
attention [23,24] due to its inherent low level of torsional 
interaction. Rapidly fluctuating loads in the vicinity of 
generating plants restrict the use of frequency signal [25], 
Larsen and Swann [26] in their three part paper have presented 
an investigation into the effectiveness of three control sig- 
nals, namely, rotor speed, electrical power and frequency. 

They recommend that stabilizers using frequency signal should 
be tuned for weak systems and those using speed and power 
signals should be tuned for strong systems. Both speed and 
electrical power signals are used for PSS design given in 
Ref. [27], They claim that the resulting PSS removes the 
inherent limitations of speed stabilizers. 



9 


Schleif et al» [25] suggested that in case of a two- 
unit plant, average speed signal cancels any component of 
swings between units (local mode). They further suggest that 
such a signal is suitable v/hen rapidly fluctuating loads (like 
arc furnace) are close to generating plants. Concordia in 
a discussion of Ref. [26] suggested that the use of average 
power signal is advisable in comparison with individual power 
signal. 

1.2.5 Design Procedure 

1.2. 5.1 Design Based on Classical Control Theory 

The problem of undamped oscillations in po’/^er systems was 
analysed by de Mello and Concordia [6] through the concept of 
synchronizing and damping torques produced in the system at the 
oscillation frequencies. A design method for stabilizers using 
frequency response techniques was suggested with the objective 
of improving damping torque without reducing the synchronizing 
torque. This led to the choice of speed signal and provision 
of a phase lead compensator in the feedback path. Bollinger 
et al. [28] used the root-locuc method for design. Multi- 
variable frequency response plots were used by Hamdan and 
Hughes [29] for the analysis and design of PSS. A good analysis 
of damping and synchronizing torques has been presented by Alden 
and Shaltout [30]- Larsen and Swann [26] used the concept of 
’plant’ characteristic for the design of FSS. They suggest a 
tuning procedure for a practical PSS used in the industry. The 
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above methods are applicable mainly to single machine 
equivalent models of the system. 

Coordinated design of Po>S in multimachine systems have 
received rauch attention recently, de Kello et al. [17] were 
the first to present the need for coordinated design of PSS in 
multimachine systems. Crenshaw et al. [31j used the procedure 
given in [26] for the design of PSS in a two unit plant. 

Rudnick et al, [32] extended the concepts of synchronizing and 
damping torques for the analysis of dynamic stability in multi- 
machine systems. How'ever, this was not utilized for PSS design, 
Gooi et al. [33j presented an iterative technique for the 
’optimal’ setting of PSS parameters, by solving nonlinear alge- 
braic equations, in multimachine system. Their technique is 
based on the computation of damping and synchronizing torque 
coefficients v/ith the help of signal flow graph. The technique 
is restricted to smaller systems. 

Classical design of PSS based on damping torque analysis 
is devised mainly for single-machine infinite-bus equivalent 
model of the power system. Although, there have been some 
efforts to extend this technique to multimachine system, this 
is not easy and application for large systems is not reported 
in the literature. 
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1.2, 5. 2 Design Based on Minimization of Performance Index 

A distinct departure from the classical control techniques 
is found in the application of modern control theory to the PSS 
design. Design methods based on optimal regulator theory have 
been suggested by many authors[34-38] . By this, the stabilizer 
is designed as an optimal, state feedback controller which 
minimizes a quadratic performance index. Yu, Vongsuria and 
Wedraan [34] were the first to apply optimal regulator theory to 
pov;er system problems. Optimal controllers were designed for 
excitation and governor control. Anderson [35] presented design 
of an optimal controller for a wide range of operating condi- 
tions, Moussa and Yu [36,37] proposed methods where the 
weighting matrix in performance index is obtained in terms of 
shift of dominant eigenvalues to the left in the complex plane, 

A slight departure from the conventional optimal regulator 
theory is found in the method proposed by Kumar and Richards 
[38]. Optimal linear feedback law is obtained to achieve 
pre-assigned eigenvalues for a closed-loop system with much 
smaller computational effort. The resulting controllers 
obtained in Refs, [34-38] improve the closed-loop system^ 
response. Hov/ever, inaccessibility of all the states for 
feedback poses problems in practical situations. 
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Yu and Siggers [39] developed suboptimal controllers with 
the system represented by lower order models with measurable 
states as state variables* Davison and Rau [40] developed a 
design method for PSS using output feedback by minimizing an 
average performance index as proposed by Levine and Athans [53]. 

De Sarkar and Rao [4l] designed excitation control for stabili- 
zing a synchronous machine through output feedback by first 
considering the state feedback obtained through the application 
of optimal control theory and then getting a minimum norm sub- 
optimal control for output feedback. 

Recently use of Kalman filter to counteract the noise in 
the control signal output due to random fluctuation of load is 
given in [42], Doraiswami et al. [43] presented the design of 
decentralized controllers for multimachine system using feedback 
signals from terminal voltage, rotor speed and power. 

Industry acceptance of optimal regulator theory is poor 
because the choice of a suitable performance index is difficult to 
make. It is not easy to judge good system performance through 
a performance index. Another reason is that direct application 
of full state feedback optimal regulator theory is not possible 
for large multimachine systems because of the decentralized 
structure of PoS. Moreover, the design of optimal dynamic out- 
put feedback controllers is not reported. 
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1.2, 5. 3 Design Based on Pole Assignment 

Improvement of damping torque is equivalent to the shift- 
ing of eigenvalues corresponding to rotor oscillations [44], 
Thus, pole assignment techniques give better understanding of 
the system damping. State or output feedback can be used for 
pole assignment, however, output feedback is more practical. 
Constant gain feedback or dynamic feedback are two alternatives 
for pole assignment. From system response point of view dyna- 
mic feedback is attractive. Most of the literature available 
have considered PS3 design as that of dynamic output feedback 
compensator [45-50], The number of poles that can be assigned 
depends uoon the order of PSS. Usually a second order PSS is 
designed to assign dominant eigenvalues of the system, 

Padiyar et aL, [45] have presented three methods for the 
design of PSS by pole assignment with output feedback. The 
first method is based on modal control theory. The second and 
third methods involve static optimization for exact assignment 
of poles and assignment into a sector. Ajrcidiacono et al, 

[46] proposed a technique for PSS design which utilizes the 
residues of open-loop transfer function to get a specified set 
of closed-loop eigenvalues. Fleming et al. [47] proposed a 
sequential algorithm for tuning the parameters of PSS with 
fixed poles in multimachine power systems, A similar pro- 
cedure is given in [20] . Sequential addition of stabilizers 
disturbs previously assigned eigenvalues and hence is not 
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satisfactory. An algorithm for simultaneously tuning the PSS 
in a multimachine system is proposed by Lefebvre [48j. Abe and 
Doi [49,50] developed a design technique by combining the fre- 
quency response method and the pole- assignment method. 

1.2. 5, 4 Design Utilizing Adaptive Control 

In all the methods discussed so far, P3S are designed at a 
particular operating point to provide adequate damping to the 
system. However, the operating point is not fixed, but keeps 
changing and the performance of stabilizers, under changed 
conditions, may not be satisfactory. The need for adaptation 
has been felt since long. The availability of low cost micro- 
processors has made it possible to design adaptive PSS. The 
use of digital adaptive regulators, called self tuning regula- 
tors, for power system stabilization has been reported recently 
in the literature [51,52]. 

1.3 CONTROLLER DESIGN 

As the techniques for PSS design presented in this thesis 
are based on modern control theory, a brief survey of the 
relevant design techniques is presented in this section. 

Although the techniques based on optimal regulator theory 
are well knov/n, mention of a few points will help in the use 
of these techniques in subsequent chapters. The solution of an 
optimal output feedback regulator problem depends on the ini- 
tial values of state variables (x^) of the system. This poses 
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problems in controller realization. Levine and Athans [53] 
were first to eliminate this dependence by assuming as a 
random variable uniformly distributed on the unit sphere. 

They presented a design technique for optimal gain controller 
using output feedback in which solution of nonlinear simulta- 
neous matrix equations is required. Later, Hole [54] and 
Horisberger and Belanger [55] modified the technique to 
avoid solution of simultaneous matrix equations. Design of 
optimal dynamic output feedback controller is, however, not 
available in the literature. 

Controller design based on pole assignment techniques 
is of later origin. Much work has been reported in this area 
[56-64]. Pole assignment can be achieved either by state or 
output feedback, using gain controller or dynamic controller. 
ViJang and Davison [6C] have shoi,vn that min (n,m+r-l) closed- 
loop poles can be assigned with a gain output feedback, for a 
completely controllable and observable system, where n,m 
and r are respectively the number of state, input and output 
variables. Brasch and Pearson [61] and Ahmari and Vacroux 
[62] have shown that for a completely controllable and obser- 
vable system defined by plant, input and output matrices 
A,B and C respectively, a pth order controller can be designed 
to assign max (a+p, P+p) clused-loop poles to desired loca- 
tions where a and § are given by 



a = Rank[B AB A^B ... A^B] 

p = Ranlc[C^ aV (A^)V... (a'‘^)P C^J 

Techniques given in Refs. [57-64] suffer from the dravA)ack 
that the resulting feedback matrix is dyadic. Due to the 
rank deficiency of the feedback controller, the resulting 
closed-loop system has poor disturbance rejection properties 
[65]. Munro and Hirbod [66] presented a systematic design 
of full rank output feedback controllers which allow arbitrary 
pole assignment in the closed-loop system. Youla et al. [67] 
presented a necessary and sufficient condition for the stabi- 
lization of a dynamical plant described by the transfer fun- 
ction matrix P(s). The plant P(s) can be stabilized by an 
asymptotically stable controller C(s) iff real poles and zeros 
of P(s) in the right half of complex plane, have the inter- 
lacing property (i.e., poles and zeros of P(s) in the right- 
half of complex plane should lie alternatively). 

Very often, it is uneconomical and impractical to globally 
(centrally) control a large, geographically dispersed, inter- 
connected physical system. Control problems associated with 
large-scale systems led to the area of decentralized control 
[68-78]. The system is divided into various subsystems with 
local control stations. At each station, the controller ob- 
serves only local system outputs and controls only local in- 
puts, All the controllers are, however j involved in coh'trol.ling 
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the same large system. A brief review on decentralized control 
is presented here. 

V!?ang and Davison [68j introduced the fixed mode concept. 
Fixed modes cannot be assigned using ’local' dynamic controllers. 
The remaining modes can be so assigned. Corfmat and Morse [69} 
gave a more complete characterization of the conditions for 
stabilization of large systems. Sacks [70] used a decoupled 
model known as component connection model (CCM) to show that 
the fixed modes of the composite system [68] and the set of 
fixed modes of the individual components coincide in CCM. 
Basically this paper shows that, within the CCM framework, 
there exist local output dynamic controllers which will assign 
a prescribed eigenvalue spectrum provided each component or 
subsystem is controllable and observable. A good bibliography 
on decentralized control is given in a survey paper by Sandell 
et al. [71]. V'Jang and Davison [72] presented a design of local 
feedback controller when all the fixed modes are not stable. It 
is shown that by limited exchange of information between local 
stations, it is possible to stabilize the entire system. Sezer 
and Huseyin [73] have given sufficient conditions based on graph 
theoretic concept for an interconnected system to be stabiliza- 
ble using only local state feedback. Xingolas et al. [74] pre- 
sented '^a two— level hierarchical computational scheme for the 
determination of decentralized gains. At lower level, two 
Liapunov equations are solved and at higher level a predictive 
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corrective routine is used. Ramakrishna and Viswanadham [75] 
have developed sufficient conditions for the decentralized 
stabilization of a class of large systems viith dynamic inter- 
connections between subsystems. These conditions are given in 
terms of the subsystem coefficient matrices. It is shown that 
if interactive subsystems are stable with mild restriction on 
the subsystems, then the large system can be stabilized with 
decentralized feedback. Viswanadham and Ramakrishna [76J have 
presented the synthesis of decentralized observer based on 
local measurements to generate maximum possible linear combina- 
tion of local states for the stabilization of the interconne- 
cted system. 

Basically the above papers provide existence theorems for 
a system to be stabilizable with decentralized control. No 
clear cut methodology is available for the design of such a 
controller. Lefebvre et al. [77] have given an algorithm to 
assign closed-loop poles by dynamic decentralized output feed- 
back. Two-level design is presented. In the first level 
closed-loop poles are assigned by local outputs neglecting the 
interconnection between subsystems. Interconnection is then 
introduced in discrete steps to update the parameters of local 
controllers v/ith the objective to retain the closed-loop eigen- 
values at their specified location. 
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1.4 CSJECTIVE AMD SCOPE OF THE THESIS 

The design of PSS can be based on pole assignment or 
minimization of a performance index provided the index is 
properly chosen. Although much work has been reported on the 
use of these techniques for PSS design, the following points 
have not been properly considered : 

(1) In pole assignment techniques full design freedom has not 
been utilized in the selection of P^^S parameters. 

(2) Design of dynamic optimal output feedback compensator has 
not been reported. Moreover, use of parameter optimiza- 
tion technique to large systems taking into account the 
decentralized nature of controller is not reported, 

(3) Although methods are available for the selection of PSS 
location in multimachine systems, their application to large 
systems is not available in the literature. 

(4) Though some work has been reported on the choice of feed- 
back control signals, a detailed study is not available, 
particularly for plants having two or more units operating 
in parallel. 

(5) A third order synchronous machine model is widely used in thj 
design of PSS. Adequacy of the model has not been investi- 
gated. 

The work reported in this thesis is directed at removing 
some of the gaps in the asoects of PSS design mentioned above. 
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The major emphasis is on evolving design procedures in large 
multimachine pov.er systems subject to the practical constraints 
of using decentralized dynamic compensators. The design obje- 
ctives are chosen as either pole assignment or minimization of 
a perfoimance index in the framework of optimal control theory. 

As adaptive control is yet to be introduced in power 
system, scope of the thesis is limited to the application of 
deterministic control. 

1.5 SUtjJARY 0? THE WORK 

A chapterwise summary of the work done in the thesis is 
given below. 

Chapter 2 discusses the problem of dynamic stability of 
present day power systems and the need for its improvement 
through the use of power system stabilizers. Different tech- 
niques based on modern control theory are developed for PSS 
design. The techniques developed here have better methodology 
and are general enough to be applied to multimachine systems. 
Two design objectives are used. Pole-assignment is used in 
the first two techniques presented. Modal control theory is 
applied in the first design technique. Only accessible states 
are used as feedback signals to PSS. The PSS dyi^araics is 
chosen in such a way that zeros of PSS transfer function are 
specified, Pole -assignment technique developed by Tiunro and 
Mirbod [66] with dynamic output feedback is utilized for PSS 



21 


design in the second technique. This technique puts no 
restriction on the parameters of PSS and their values are 
directly obtained in terms of the desired locations of closed- 
loop eigenvalues. For complete pole assignment a two-input 
dynamic compensator is designed. In the techniques developed, 
PSo parameters are directly computed, unlike trial and error 
method used in classical design. 

l.linimization of a performance index is another objective 
which is used in PSS design. Electromechanical oscillations 
of machine rotors with respect to one another cause fluctua- 
tions in power output of machines and limit power transfer 
capabilities of the systems. This effect is of great concern 
in systems having long transmission lines. Damping of these 
oscillations can bo achieved by minimizing the fluctuations in 
power output of generators. The third technique is based on 
parameter optimization using output feedback with the objective 
to minimize power fluctuations in generator outputs following 
a disturbance. These techniques are used to design PSS for a 
single-machine infinite-bus system. Comparison is made of the 
performance of these PSS with the PSS designed using the classi- 
cal technique used in the power industry. 

Chapter 3 presents a detailed investigation into the eff- 
ect of synchronous machine modelling on the design of PSS. 

Two different mcdels are considered for the investigation. 

The simplified or the lower order model considers only the 
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effect of field flux dec 3 y in addition to the equation of 
mechanical motion. Effect of amortisseur windings are 
neglected in this representation. The detailed or higher 
order model considers the damper circuits in addition to the 
lower order representation [78]- PSS are designed, for both 
representations of synchronous machine, at different machine 
loading, inertia constant and voltage regulator gain, to 
determine the effect of machine representation on PSS para- 
meters. Adequacy of the simplified model for PSS design is 
analysed by incorporating PSS, designed for lovijer order model, 
to the machine represented in detail. Performances of the 
stabilizers, designed for both representations of machine, are 
compared by eigenvalue and time-domain analyses. 

Chapter 4 presents a detailed analysis on the effective- 
ness of different feedback control signals that can be used 
for PSS design. Rotor speed, electrical power output and 
terminal bus frequency signals are considered. A generating 
plant, having a number of generating units, connected to a 
large power system is considered for the analysis. First, 
the effectiveness of control signals is studied for a plant 
having one machine. Then, the interaction between machines 
present in the plant on the effectiveness of control signals 
is investigated by considering another machine operating in 
parallel with the existing machine. Effectiveness of a PSS 
utilizing a particular feedback signal is judged by determin- 
ing the stable P-Q region for the machine equipped v/ith that 
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PoS. The influence of ac system strength on the design of P3S 
utilizing the three control signals is studied. The suitabi- 
lity of average speed and average power signals is also presen- 
ted in the case of plant v/ith more than one machine. 

Chapter 5 presents the development of a state-space model 
for large power systems, suitable for dynamic stability analy- 
sis, and the investigation into the selection of effective 
locations of P33 in such a system. The development of the 
system model proceeds systematically by developing the models 
of various machines and their interconnection through the net- 
work model. The network is represented by its Jacobian matrix. 
In the development of system model it is not necessary to eli- 
minate the non-generator buses and the identity of each bus of 
the system is retained. Changes in the system matrix caused by 
the changes in the system configuration can be easily accommo- 
dated. 

A proper choice of Poo location is essential for the 
coordinated operation of multimachine systems. Two analytical 
techniques are presented for the choice of effective locations 
of Poo. The first technique is independent of controller stru- 
cture and the second is based on eigenvalue sensitivity, A 
case study of a 13-machine system is presented for the selection 
of PSS locations. The results obtained through analytical 
techniques are validated from the loci of eigenvalues. 



24 


Chapter 6 presents the application of decentralized 
control with output feedback for the design of P3S in multi- 
machine pov^er systems. Techniques based on eigenvalue assign- 
ment and minimization of a performance index are presented for 
the design of PSS, 

The technique based on eigenvalue assignment utilizes the 
algorithm given in Ref. [77j. PSS for the selected machines 
are designed, neglecting the interconnections between machines, 
to assign closed-loop eigenvalues. Parameters of the stabili- 
zers are then updated by introducing the interconnections, in 
steps, keeping the closed-loop eigenvalues unchanged. The 
advantage of this itertive method is that stabilizers are 
designed for smaller subsystems, easy to handle in comparison 
with the design of stabilizers for the entire system considered 
together. In the second technique a performance index, based 
on minimization of power fluctuations in the output of rele- 
vant generators, is used to design PSS. Steeoest descent 
method is used to update the parameters of stabilizers in the 
optimization routine, A three machine and a thirteen machine - 
examples are considered to illustrate the proposed techniques 
and to compare the results obtained. 

A brief review of the major contributions of this thesis 
and possibilities for further work are given in the concluding 
chapter. 
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DYNA’ilC COI.iPEKSATOR DESIGN FCR POWER SYSTEK! STABILIZATIOI'i 

2.1 INTRODUCTION 

Since early 1950, most of the nev/ generating units added 
to electric utility systems have been equipped with continuously 
acting voltage regulators. As these units increased in number, 
it became apparent that the voltage regulator action had a 
detrimenttil. effect on the dynamic stability of power systems. 
Oscillations of small magnitude and lov/ frequency often persis- 
ted for long periods of time and in some cases presented limita- 
tions on power transfer capability [l ,2 ]• Damping of these 
oscillations can be achieved by the use of power system stabi- 
lizers (PSS). PSo are auxiliary feedback controllers which 
receive signal from rotor speed, electrical power output or 
terminal bus frequency and superimpose a corrective signal on 
the voltage error signal of the normal voltage regulators. 

Design of PSS has been conventionally carried out using 
classical control theory [ 6,26j. This approach has its limi- 
tation in that it is restricted to systems which can be repre- 
sented as a single machine connected to an infinite bus. This 
model of power system is not adequate in many situations and 
coordinated application of PSS using multimachine system model 
is essential [17]. Thus, it is possible that PSS designed on 
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the basis of single machine infinite bus system may introduce 
negative damping for some modes of oscillations while damping 
some others [l7]. 

ilodern control theory provides a better tool for the design 
of PSS in large multimachine power systems [34-48] than classi- 
cal control theory. Modern approach to PS3 design is based 
either on pole assignment techniques [44-48] or on optimal 
regulator theory [34-41] . Pole assignment can be achieved via 
state feedback, but inaccessibility of state variables necessi- 
tates the use of output feedback. Industry application of 
optimal regulator theory is somewhat poor because of lack of 
feel for performance index. Pole assignment techniques, on the 
other hand, give better understanding of system damping. It 
is to be noted, however, that in the works reported, design 
freedom has not been fully utilized. 

The basic aim of this chapter is to provide techniques, 
with better methodology, for the design of PSo as a dynamic 
compensator, fully utilizing the freedom available in the 
choice of its parameters. Three design techniques are presen- 
ted based on pole assignment and minimization of a performance 
index. The first technique utilizes partial state feedback and 
is based on modal control [56] to assign dominant eigenvalues 
corresponding to rotor oscillations. The second technique 
utilizes output feedback and follows the algorithm of Munro and 
Hirbod [66] with the objective of assigning closed-loop poles. 
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Design objective in the third technique is minimization of 
fluctuations in pov/er output of generators. The technique 
is based on the works reported in Refs. [53»54j* 

Although the PJS design techniques presented can be 
applied to multimachine power systems, their application is 
restricted to single machine infinite bus system in this 
chapter. Chapter 6 deals with their application to multi- 
machine power systems. 

2.2 PROBLBl OF Pd^ER SYSTBl CTABILIZATIOI^J 
2.2.1 General 

The problem of power system stabilization can be explained 
with the help of a single-machine infinite-bus system, shown 
in Fig, 2.1. This representation of power systems has been 
considered by many authors [ 1-7 j for the design of Poo, As 
discussed in Chapter 1, for the purpose of dynamic stability 
analysis, linearized state space model of the system is used. 

In this section a state space model of the system given 
in Fig, 2.1 is developed and with the help of eigenvalue 
analysis, dynamic behaviour of the system is analysed. The 
need of a PS3 to stabilize the system and the objective of its 
design is discussed in teims of the concepts used by power 
engineers. 




FIG. 2.1 SINGLE MACHINE INFINITE BUS SYSTEM. 



FIG. 2.2 SINGLE-PHASE EQUIVALENT CIRCUIT OF THE 
MODEL . 



FIG. 2.3 BLOCK DIAGRAM OF EXCITATION SYSTEM. 
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2.2.2 Development of System Model 

The following assumptions are made for developing thc. 
state-space model of the system given in Fig, 2,1 : 

1, The synchronous machine has no amortisseur windings. 

2, Governor and turbine dynamics can be ignored. 

3, The machine stator and the external network are in 
quasi*-steady state. 

2. 2, 2.1 Machine ?>iodel 


A three winding model of the machine is considered. 
Stator of the machine is represented by a dependent current 
source in conjunction with a constant impedance [78], DerivS' 
tion of the machine model is given in Appendix D. The system 
dynamical equations are as follov/s ; 


pw 


A fT ~T ) 


^0 

M” 




( 2 . 1 ) 


P5 = (w - (2,2) 

and 

pid = vpYo 

where 

• I = 2! 


= "I « 







q 


(2.4) 


§ 4 (L^ - L^)/L^ 


(2.5) 
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The non-state variables i ,, i and I in eqns. (2.1) to 

(2.4) are eliminated by representing them in terms of the state 

variables and 6 as given below. Figure 2,2 represents the 

single phase equivalent circuit of the power system shown in 

A 

Fig, 2.1. From Fig. 2.2, the armature current phasor i^^ can be 
written as t 


A 

1 




® 2eq-"j='d 




V, 




Substituting the values 

5 ^ / o *4- *1 o \ T *1 T ^ 


^ A 

of I, Z and , we get 
eq eq 




eq 


( 2 . 6 ) 


where 


a,+jao ^ 




2 ~ R- 


and 


R+ j (X+x^T 


Substituting the value of from eqn. (2.5) into eqn, 
(2.6) and equating the real and imaginary parts, we get 


®2 ’^ec 

iq = ~ T-ly ^^1 ^ ^2 ^ 


(2.7) 


and 


Bn § 


[aoln - cos6 + hn sine )3 


Vd ■ l^ITl '-=‘2^d T "eq'^l “2 


+ (bo cos6 - b Sins ) 
eq ^ 2 i 


( 2 . 8 ) 
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2. 2,2. 2 Excitation oystem Model 

A single time-constant excitation system, IEEE Type IS 
[9 ], is considered. Figure 2.3 gives the block diagram 
representation of the excitation system. The state equation 
of the excitation system as obtained from Fig, 2,3 is 

K- 


p E 


fd 


1 

T^ ^fd ^ T^ ^\ef - ’-'t ^ 


(2.9) 


where = (v^ + 


, 2 
'd 


t ‘ 'd ■ 'q' 

The Park's components of terminal voltage, v, and v are 
obtained from the equivalent circuit of Fig. 2,2 as ; 


v^ = VV„ - JV,) e-^^ 


w. 




Equating the real and imaginary parts, v»e get 


and 


>^d<l -5)iq 


''q = ’'d^-d - ^d) 


2. 2, 2, 3 Linearized Model 


The linearized state space model is derived by lineari- 
zing eqns. (2.1) to (2,3) and (2.9). The non-state variables 
Ai, and M are eliminated using the foilov«;ing expressions 

tpA 'Ll 

obtained by linearizing eqns. (2.7) and ( 2 . 8 ), 

Aiq = Kj_A + K2A6 

~ ^3^ ^d ^4^^ 


(2.10) 
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A 

wli0ir0 1C -I ~ 

1 l-aj^T- 


and 


V. 


^2 = 1^1 ^o " 4 


K ^ a 

^3 - - I-IYI 


32 S 


K4 - '^ eqC -(‘=2 + *>1 TTa '- s ) *0 ~ 
^‘’1 ■ ‘’2 r?55' > ®oi 


The resulting fourth order model of the system is des- 
cribed by the state equation 


P - [a] X + _b u (2.11) 

where 


X 

A 

[ AW 

A6 AE^^] 


u 

4 

Au' 






^11 

^12 ®13 

0 

[A] 



1 

0 0 

0 




0 

^32 ®33 

®34 




0 

®42 ®43 

^44 


b = [O O O b^j^ 
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‘11 


Kq/a 


X 


12 “ m' ^ ^4 " ^^do“ ^ 


^d 


*13 


p [( Sko-I) iqQ + ( ^ iclo'”^do^^l-‘ 


*32 


(>^ci ~ ""d^ 


X’, T', 

Q GO 


Xj + ( X j ~ Xj) 


d' ‘^3 


‘33 


■^do 


‘34 




>42=?^'^ tv K4 - Vj^d-nKj] 

% '^to 


‘43 


h.lL 

Tr ^to 




^44 “ "■■^/^p 


^4 =%/T^ 


2. 2. 2 .4 Numerical Example 

A single-machine infinite-bus system as shO¥m in Fig. 2.1 
is considered. The following system datas given in Ref. [l3], 




are taken for the study. The machine and transmission line 
parameters are given in per unit on machine base. 

Generator parameters • ~ “ 1*63 

= 6.3 sec; H =4,0 sec. 

£. 0 ^ = 3i4rad/sec. 

Transmission Line parameters : Rj. = 0.024 ; = 0.6 

G 3 = 0 ; Bg = 0.066 

Voltage Regulator parameters : = 50. ; Tr, = 0.02 sec. 

i'k It 

Operating data : Generated power = 1.0 p.u. at 0.9 power 

factor lagging 

Infinite bus voltage = 1.0 / 0 ° 

The matrices A and b of the system at the operating point 
are as follows : 


~ 0 
1 

[A]= 

0 

- 0 

and 

b = [0 


-35.95933 -15.34320 O 

0 0 0 

-0.38083 -0.34625 0.35273 

240.33139 -544.65848 -50 

0 0 2500 ]^ 


Open-loop eigenvalues of the system matrix A are given in 
Table 2.1. The complex pair of eigenvalues, shown in Table 2.1, 
corresponds to the rotor oscillation mode. As the real part of 
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this eigenvalue is positive, the system mil be dynamically 
unstable and there is a need to stabilize the system by employ- 
ing a PS3. 

Table 2.1 Open-Loop eigenvalues 


0.20788 + j 6. 11383 
-5 .00647 
-45.75554 ’ 


2.2.3 Power System Stabilizer 

The concepts used by power engineers for the design of PSS 
can be explained with the help of the pioneering work of 
de iVlello and Concordia [6 ]. Stability analysis given in [6 ] 
is based on the concept of synchronizing and damping torques. 
Synchronizing torque is defined as the component of electrical 
torque which is in phase with the rotor angle phasor. Damping 
torque is representeo. by the out of phase component, namely, the 
component in phase with rotor velocity. Synchronizing torques 
are basically provided by the electrical ties between machines 
and are functions of machine operating conditions and network 
characteristics. Damping torques are provided from different 
sources : machine danper windings, field windings, excitation 
system (AVR controller) and speed governors [32j. Damping 
torques produced by first t \¥0 sources are positive vAile 



36 


contributions of .aVFI controllers and speed governors may be 
negative. According to stability criteria given in [6 ], both 
the damping and synchronizing torques should be positive for 
the system to be stable. It is showri in [6 ] that excitation 
system has conflicting effects on the synchronizing and damping 
torques. Under heavy loading conditions the effect of excita- 
tion system is to increase the synchronizing torque in the 
system and reduce the damping torque. The reduction in dancing 
torques, particularly in systems with long transmission lines, 
results in extremely oscillatory behaviour of the system. The 
oroblem is solved by providing a torque component, through 
external means, in phase with rotor velocity such that the 
reduction in damping torque is compensated. This forms the 
basis for the design of PS3 and is widely used in po\^r industry. 

It can be shov^a [44] that the damping and synchronizing 
torques are approximately related tcc real and imaginary parts 
of the eigenvalues corresponding to rotor oscillations. The 
design objective of PSS given in [ 6 ] can, therefore, be trans- 
lated into that of shifting the real parts of the dominant 
eigenvalues corresponding to rotor oscillations. Thus , the 
design of PSS can be treated as a problem of pole-assignment. 

Before discussing the design techniques for PSS, the fun- 
ctions of various components of a practical PSS using speed 
signal is explained with the help of the block diagram given 
In Pig. 2.4* 
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Vash o ut Sta ge 

The washout stage is provided to prevent offsets in ter- 
minal voltage for the steady-state change in system frequency. 
The washout time constant T^ should be large enough to pass 
signals at the required stability frequencies relatively 
unchanged . 


G ain 


The amount of damping associated with rotor oscillations 
depends on stabilizer gain The damping increases mth 

increase in K upto a certain value, beyond which, further 

o 

increase in gain results in a decrease in damping. The sta- 
bilizer gain should be set at the value corresponding to the 
maximum damping of the rotor oscillations. 


Lead/La q S tages 


To provide damping torque for oscillation of the machine 
with respect to the power system, it is necessary to develop an 
electrical torque which is in phase with speed changes. Lilith 
speed as input signal, it is necessary to provide a phase lead 
to compensate for lags between regulator input and the result- 
ing electrical torque. The function of this stage is to 
provide the required phase lead over the desired range of 
frequencies. 
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Xo^L?ijon_a _ Fi jU e r 

Torsional filter is used to attenuate the stabilizer gain 
at high frequencies to limit the impact of noise and to mini- 
mize torsional interaction between generator and turbine. 

Stabi li zer Limit 

To allow maximum contribution of the stabilizer, a large 
positive limit of about 0.2 p.u. is used. On the negative side, 
a limit of 0.05 to 0.1 p.u. is used to allow sufficient control 
range and to reduce the probability of a unit trip for failure 
of a stabilizer component. 

As the influence of washout, filter and limiter stages is 
very small in the range of rotor oscillation frequencies, the 
design of a practical P33 can be considered as the determina- 
tion of gain and time constants of lead/lag stages. In general, 
the design of PSS can be viewed as the design of a dynamic 
compensator for stabilization of the system. 

2.3 DHSIGiT OF DYAAviIC C0,;i^EM3AT0R 

As discussed in the previous section, pov/er system stabi- 
lizers are basically dynamic compensators of a fixed order. The 
control signals to the PSS are appropriate measurable outputs 
of the system such as rotor speed, electrical power output or 
bus frequency. Thus, the design of PSS can be considered as 
the design of output feedback dynamic compensator. The aim 
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of this section is to develop various techniques, based on 
modern control theory, for the design of dynamic compensators. 
These techniques v/ill be subsequently used for the P5S design. 

The first technique utilizes modal control [56 j with 
partial state feedback to assign closed-loop eigenvalues. 
Compensator dynamics and control lav/ are assumed in such a way 
that the zeros of compensator transfer function are fixed. 

Thus, gain and poles of the compensator are to be determined. 

The second technique utilizes the algorithm suegested by /nunro 
and Kirbod [66]. The design objective is to assign closed-loop 
eigenvalues with output feedback. Compensator structure is 
fixed and the parameters of the compensator, namely, gain, 
poles and zeros are obtained to meet the design objective. The 
third technique is based on parameter optimization to minimize 
a performance index using output feedback [53,54]. 

2.3.1 Pole Assignment with Partial State Feedback Compensator 

2. 3. 1.1 Formulation of Control Problem 

Consider the single input linear time invariant system 

p x = [A] X + b u (2.12) 

where x e and u eR^. Let _x' represent the vector of acce- 
ssible state variables, then 


X* = [Cj X 


(2.13) 
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where x‘ £ R^, r < n, and C is an r x n matrix. Without loss 
of generality C can be considered to be of the form 

[C]= [G I I^] (2.14) 

where is an r x r identity matrix and C represents an 
r x(n-'-r)null matrix. 

Consider a kth order dynamic controller and the control 
law of the form 


p z = [Dj z + e 


u 


and 


u = kJ X' + k| 2 


(2.15) 

(2.16) 


where [D'j = diag [d^^ d 2 


and 


= [11 ... 1 ]' 


The structure of D and e fixes zeros of the controller at 
di,d 2 , ... and dj^ (Appendix A). and K 2 are respectively 

r X 1 and k x 1 matrices, representing feedback gains. Com-- 
bining eqns. (2.12) and (2.15) and substituting for x’ from 
eqn. (2.13) into eqn. (2.16), v;e get 


A 
p X 


r V ^ V 

= [Aj X + b u 


and 


u = 


45 


(2.17) 

(2.18) 


where 




0 


0 D 


A 

b 


A 

X 


X 
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k! = [o 


^ : 4 ] 


(2.19) 


Eqn. (2.17) represersts the dynamical equations of the 
composite system of order n^^, where n^^ = (n+k). in eqn. 
(2.18) is an n^xl matrix of state feedback gain which from 
modal control theory [56] can be obtained as 


K„ = [W] £ 


( 2 . 20 ) 


where ¥? is n^^ x n^^ matrix forraed by the reciprocal eigen- 

A 

vectors of matrix A corresponding to the open-loop eigen- 
values expression for the jth element 

of G is given as [56 J 


(X,- p.) 


( 2 . 21 ) 


where 




and P, , ...» P are the desired locations of the closed- 

_L ^ ^ X 

loop eigenvalues. From eqns. (2,19) and (2,20), it is clear 
that entries of a should be such that first (n-r) elements of 
are equal to zero. Now, as the elements of a are functions 
of open-loop and closed-loop eigsenvalues and as they have to 
satisfy the above constraint, all the clos-ed-loop eigenvalues 
cannot be assigned arbitrarily. In other words, if some of 
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the closed-loop eigenvalues are to be assigned to specified 
locations, there is no control on the location of remaining 
closed-loop eigenvalues. Let n 2 eigenvalues be assigned and 
t = - 02 ) eigenvalues are allowed to assume arbitrary 

locations. Expanding eqn. (2,21) 


~ ^2^ ••• ^ ^t+1^ ^t+2^^ 

1 2 ... t 12 


... (Xj - A )] 

^2 

Pi, P 2 * are unknown while Pt-f2' **** ^n 

specified, a. can be written as 


G. = S'. + q, 

J 1 J J 




( 2 . 22 ) 


where ~ P^^2) “* - P,. ) 


n. 


and qj^, q 2 , ..., are coefficients of polynomial in and 
functions of Pj^, P 2 , ..., P^. In vector matrix notation eqn. 

(2.22) can be written as 



(2.23) 
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q = [1 q 2 ... 

and 

~ - 1 2 

oubstituting the value of a from eqn. (2.23) into eqn. (2.2C) 


q (2.24) 

where 

[3j = [V/j [P] 

Partitioning of eqn. (2.24) results in 

q (2.25) 

considering the first set of (n-r) equations in eqn. (2.25), 
the order of the submatrix. 3^ will be (n-r) x (t+1). These 
equations can be rearranged as 


^0 1 


s 


0 


= 

^1 

_K2 


_ ^ 2 _ 


S’ 



where 


""-ol 




-ot+l^ 


q’ 




(2.26) 


and s^ 2 » t+1 

(2.26) is a set of (n-r) 


are columns of the submatrix 3^. Eqn. 
simultaneous linear equations in t 
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unknowns ¥/hich can be solved for qg . . . q^ and in turn for 
Pf, P 21 P^, the unassigned closed-loop eigenvalues. The 

necessary and sufficient condition for the existence of the 
solution of eqn. (2.26) is [66 ]• 


Rank [o' j = Rank [o' s'] (2.27) 

*-o 0 0 

oubstitution of the solution vector q in eqn. (2.24) gives the 
feedback gain vector The parameters of the compensator can 

be obtained in terms of and matrix D (see Appendix A). 

2. 3. 1.2 Algorithm 

Parameters of the compensator of a fixed order can be 
obtained as follows : 

Step 1 i Select zeros of the compensator by defining the 
matrix D in eqn. (2.15) and augment the system 
equations as given in eqn. (2.17). 

Step 2 i Construct matrix i of eqn. (2.24) for different 
values of t and from the condition given in eqn. 
(2.27) find out minimum value of t, the number 
of unassigned closed-loop poles. 

Step 3 t Solve eqn. (2.26) for q (i.e., 0^2 ... q^) and 

find out the locations of unassigned closed-loop 
poles. 

Step 4 ; From eqn. (2,24) evaluate the feedback gain vector 

Step 5 i Calculate the parameters of dynamic compensator: gain 
and poles in terms of and D as given in eqn, (A. 9) 
oi" Appendix A, 
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has a desired set of poles. The kth order compensator has 
the form 

F(s) = [M^(s) ... M^(s)] (2.33) 

where iVij^^Cs), i = 1, r are polynomials of degree k and 

A (s) is the characteristic polynomial of F{s). 

If is the closed-loop characteristic polynomial 

with the desired set of poles, then the basic equation which 
relates the parameters of compensator (i.e., coefficients of 
A_(s) and M. (s), i = 1, ..., r) and the desired locations of 

O JL 

closed-loop poles is [66 j 

A^(s) = ^o(s) \(s) + 2^ U^is) :i^(s) (2.34) 

If q is the number of specified closed-loop poles ^*/ith 
0 < q <_ n-fk and t =(n+k-q)is the number of unspecified poles, 
then ^c[(s) can be expressed as 

= (s'^ + d^ s'-*"^ + ... + d^,)(s^ + ej_s^"^+...+ e^) 

(2.35) 

where d- , i = 1, ..., q are known and e- , i = 1, ... t are 

i i- 

unknown. Equating the coefficients of like power of s in 
eqns. (2,34) and (2.35), v/e get 

[Xj^-J (2.36) 

and 6^ are defined in eqn. (B.23) of Appendix 3. 


where X^, 
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Hqn. (2.36) is a set of (n+k) equations in [r(k+l) + k-rtj 
unknovms. Vector contains the parameters of F(s) and the 
coefficients of the polynomial of unassigned closed-loop 
poles. The solution of eqn. (2.36) gives the desired compen- 
sator and the location of unassigned closed-loop poles. Two 
questions can be asked at this stage. 

(1) For the system defined by eqns. (2.28) and (2.29) j what 
should be the minimum order of the compensator (i.e., the 
minimum value of k) , v;hich assigns all closed-loop poles 
(i.e., q= n+k)? 

(2) For a given value of k, how many poles can be assigned? 

The answer to these questions can be obtained from the 
necessary and sufficient condition given in eqn. (B.15), 
Appendix B, for existence of solution to the pole assignment 
problem, whether it is partial or full assignment : 

Rank[X^] = Rank[X^, 6-^] (2.37) 

2. 3. 2. 2 Algorithm 

The following steps are to be follov/ed for the design of 
output feedback compensator to assign closed-loop poles. 

Step 1 ; Construct the (r x 1) open-loop transfer function 
matrix G(s) 

Step 2 t If k is specified, determine the maximum va'iue of q 
for which eqn. (2.37) is satisfied by defining dj^, 
d2, (2.39). 
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Step 3 5 If q is specified, determine minimum value of k 
for which eqn. (2.37) is satisfied. 

Step 4 s Solve eqn. (2.36) for p^ to get the parameters of 

F(s) and locations of unassigned closed-loop poles. 

2.3,3 Design of Optimal Controller for Cutout Feedback Systems 

2. 3. 3.1 rormulation of Control Problem 

Consider the linear, time-invariant, controllable and 
observable systems 

p _x = [a] X +[b]u (2,38) 

and = [*-J (2.39) 

where x e R^, u and y e are respectively vectors of state 
input and output variables. The problem is to design a con- 
troller using output feedback v^ith control law of the form 

u = -[K] y (2.40) 

such that a quadratic cost function 

^ , 

J = / Q X + H ^ jy) (2.41) 

o 

is minimized, where Q is a constant n x n positive semi- 

definite matrix and R is a constant m x m positive definite 

matrix. From eqns. (2,38) to (2.40), the closed-loop system 

and the cost function can be expressed as 

P “ 


(2.42) 
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and 

“ t + -t- 

J = / (Q + R K C) X dt (2.43) 

o 

where 

= A - 3KC 

As the solution of eqn. (2.42) depends on the initial 
state x(o), the cost function J given by eqn. (2.43) is a 
function of x(o). The dependence of cost on x(o) can be 
eliminated by assuming the initial state to be a random 
variable uniformly distributed on the surface of n~ 
dimensional unit sphere [ss]. With this assumption, the 
problem is reduced to the determination of the optimal 
feedback matrix K which minimizes the cost function 

J = trace [Oj (2.44) 

where 6 is the cost matrix and is the unique positive defi- 
nite solution of Liapunov equation 

aJ 3 + 3A^ + (Q + cVrKC) = O (2.45) 

Derivation of eqns. (2.44) and (2.45) is given in Appendix C. 
Defining the Harailtonian corresponding to eqns. (2.44) 

and (2.45) as 

H = trace [Sj + trace[p'^(A^o t SA^ + Q + cV RKC)] (2.46) 

the following necessary conditions are obtained for minimum 
value of j[54 ] : 

(1) Zqn. (2.45) 

(2) P.lJ + A^P +1=0 


(2 .47) 



(3) 


(2.48) 


II -= 2(RKCPC^ ~ B^SPC^) = 0 

Aopenciix C gives an iterative procedure, utilising the 
above necessary conditions for optimality, to determine the 
optimal feedback matrix K. 


2,3.4 Comparison of Methods 

The design objective in the first two techniques 
presented here is the sane t to assign poles of the closed- 
loop system. In both techniques, relationship between 
compensator oararaeters and closed-loop poles is used to 
obtain compensator parameters. No trial and error method is 
to be applied as in the case of classical control theory. 

In the first technique zeros of the compensators are 
specified. For a second order single input corapensator , (i.e. 
with one accessible state), three parameters, namely, gain 
and two poles of the compensator are to be obtained. Since 
zeros of corapensator are specified apriori, the design does 
not utilize full freedom and is therefore restricted. 

In the second technique no restriction is placed on the 
parameters of compensator. For a second order single input 
compensator all the five parameters - gain, poles and zeros 
are to be determined. This results in a less restrictive 
design and more closed— loop poles can be assigned in com- 
parison with the first technique. 
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The third technique is based on parameter optimization 
vdnere parameters of compensator are determined to lainimize a 
quadratic performance index- The selection of proper loca- 
tions of closed-loop poles, vjhich is a difficult task in a 
practical system, is not required here. Closed-loop poles 
are assigned to stable locations in the process of optimiza- 
tion. If it is desired to define a sector for closed-loop 
poles, then the technique can be modified to solve the 
constrained optimization problem. 

2.4 nmEBlCAL EXAMPLE 

The example of a single-machine connected to an infinite 
bus, given in Sec. 2. 2. 2. 4, is considered. PS3 are designed 
for this system with the three techniques suggested in 
Sec. 2.3. Neglecting the washout and filter stage in Fig, 
2.4, the transfer function of a practical P3S used in power 
industry can be given as : 

(1+sT, )(l+sT-,) 

F(s) = K (2.49) 

^ (1+sT2)(1+sT4} 

which is equivalent to 

0_^s^ + 9, s + Q, 

F(s) =-■% 

s + + Y 2 


(2.50) 
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where 


T T 




T2 T4 


l2l.^4 

12 “ t; 


^2 T- t; 


The corresponding output feedback matrix for the augmented 
system (machine and PSS) can be given as (see Appendix A). 


CK]= 


0 0 


i -Yc 


where 


9i = - S^Yi 


BX = 0^-8 
"^2 2 ( 


2.4.1 Design of P3S with Partial State Feedback 

Change in rotor speed Aw, which is one of the state 
variables of the system is used as feedback si^al* Zeros of 
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the P3S transfer function are specified as -1.5 and -2;5. 
A maximum of three closed"loop eigenvalues can be assigned 
by this technique to satisfy the condition given in eqn . 
(2,27). The desired locations of these three poles are 
selected as -1 + j7.0 and “2,0. The feedback matrices 
and K 2 and the corresponding PSS transfer function are 
obtained as 

[kJ j k|] = [0.09146 I -2.37908 -2,10896] 

F(s) = Q- Q9146 - r 0.36584 s + 0 .34298 
s^ + 8.4S804 s + 12.86115 

The remaining three closed-loop eigenvalues alongwith the 
assigned three are given in the first column of Table 2,2, 


Table 2,2 Closed-Loop Eigenvalues of the System 



'Tith Partial 

’(Jith Output 

Feedback 


State Feedback 

q = 5 

q = 3 

Desired 
Location of 
Eigenvalues 

-1 + j7.0 

-2.0 

-1 + 37.00 

-2.0 

-1.0 

-6.0 

-1.0 + 37.0 

-2.0 

Location of 

Remaining 

Eigenvalues 

-4.18011 + 
j2. 49997 

-46.47408 

-46.26379 

—4 . 18012 + 
32.49997 “■ 

-46 ,47409 
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2 ,4.2 Design of PSS with Output Feedback 

2. 4. 2.1 Single Input Compensator 

Change in rotor speed is used as outout feedback signal. 
It can be seen from the condition of eqn. (2.37) that, with 
a second order PSS, at the most five closed-loop eigenvalues 
can be assigned. If zeros of PSS are specified as in 
Sec. 2.4.1, then a maximum of three closed-loop eigenvalues 
can be assigned. The following tvro cases are considered here : 

(a) Five closed-loop eigenvalues are assigned at -1 j; j7.0, 
-1.0, -2.0 and -5.0 and no restriction is placed on the 
zeros of PSS . 

(b) Zeros of PSS are specified at -1.5 and -2.5 as in 
Sec. 2.4.1 and three closed-loop eigenvalues are 
assigned at -1.0 + j7.0 and -2.0. 

The transfer function of PSS and the output feedback 
matrix for the above two cases are obtained as follows ; 


Case (a ) 


s ) 


0.06501 
' ^2 


s^. + 0 .171 755 - 1.2 2838 

+ 6.91754 s + 3.23809 


[K] = 


0.06501 


0 


>1 .43940 


0 


- 0.27798 


1 


-3.23809 


-6.91754 
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p(s) = 0.0914 6 ^ 0.36584 s + G. 34298 

+ 8.48806 s + 12.86118 

0,09146 -0.83332 

[k]= 0 0 

_1 -12.86118 

The remaining closed— loop eigenvalues alongv^ith assigned 
closed-loop eigenvalues for the above two cases are respecti- 
vely given in the second and third columns of Table 2.2. 

2.4.2 .2 Two Input Compensator 

In the previous section it was observed that all the six 
closed-loop eigenvalues cannot be assigned with a single input 
second order PSS. At the most, five eigenvalues can be assign- 
ed and there is no control on the location of the remaining 
eigenvalue. All the eigenvalues of the system can be assigned 
either by using a compensator of higher order or by using a two 
input compensator. In order to avoid the complexities of 
increasing the size of closed— loop system, two input compensa- 
tor is designed. Rotor speed was used as the control signal 
for PSS design in the previous sections . The second signal for 
•j-j-jg two input compensator is selected as the terminal voltage 
of the generator. The block diagram of the system with two 
input compensator is shovffi in ^ig * 2.5, 


■O.4i048~l 


■8 .48306 
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Synchronous machine is represented by the model discu- 
ssed in Sec. 2.2.2. The excitation system used here is shown 
in Fig. 2.5 with the value of as 0.005 sec. Machine and 
line parameters and system operating conditions are the sar* 
as in the previous section. The matrices A, 3 and C obtained 
are as follows : 


[A]= 


0 

1 

0 

0 


-35.36912 

0 

”0.13571 

0 


-26 .43598 


-0.35072 


0 


O 


0.15873 


200.0 


b, = 

[0 

0 

0 200.0]"*^ 

[C]= 

~0 

0 .20863 

-0.54971 0 - 


- 1 

0 

0 0 _ 


The open-loop eigenvalues of the system are given in the 
first column of Table 2.3. The desired locations of closed- 
loop eigenvalues are given in the second column of Table 2.3. 


Table 2.3 Eigenvalues of the System with Two Input 
Compensator 


Open-loop Eigenvalues 

-0.05063 + 35.94486 
-0.249466 
- 200 .00 


Closed-loop eigenvalues 


-1.05 + j5.945 
-1 .25 
- 200.00 
- 5.00 
-10.CX3 
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The two input compensator is designed using the 
technique given in Sec. 2.3.2 with r = 2. The parameters 
of the two input compensator to be determined are : 

®T[^Q» ^20’ ^21' ^22’ ^2* this case eqn. 

(2.36 ) is a set of six equations in eight unknown parameters 

of the compensator. To satisfy the condition given in eqn. 
(2.41), two parameters are to be specified. Let both the 
poles of the compensator be placed at -10. The transfer fun- 
ctions of the two input compensator are obtained as 

p (s) = 0-95 406 s^ 4- 197. 9891 s + 39 4.4889. 

^ s^ + 20.0 s + 100 

p /g) = Q-*0 5621 s^ + 20.28633 s + 277.2835 
^ s^ + 20.0 s + 1CX).0 

The results are discussed in Sec. 2.4.5. 

2.4.3 Design of PSS with Pararaeter Optimization 

Low frequency electromechanical oscillations have detri- 
mental effect on the dynamic stability of power systems, 
particularly in systems having long transmission lines. This 
is due to the continuous variation of rotor angles which in 
turn results in fluctuations in power output of the generators 
leading to fluctuation of power in transmission lines. Thus, 
daraping of rotor oscillations which is the basic function of a 
PSS can be viewed as the minimization of fluctuation in power 
output of machine. 
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PS3 of the form given in eqn . (2.50) is designed using 
the parameter optimization technique given in Sec. 2.3.3. 

2. 4. 3.1 Representation of Fluctuation in Povjer Output of 
Generator 

The electrical power output of machine in per unit can be 
obtained from eqn. (2.4) as 


Linearizing eqn. (2.51) and substituting the values of 
and Ai^ from eqn. (2.10), we get 


(2.51) 
A i 

q 


Ap = K. AI, + ICA6 (2.52) 

e 5 d D 

where 

% ^d ^^qo %^^do ■* ^ ^qo^ 

and 

^6 ~ ^d ^do^ “ K4 § I^q] 


Eqn. (2.52) can also be written as 
APg = [TJ X 

where [t] = [O O] 


(2.53) 


For the example considered matrix 1 is given as 
T = [O 0.91570 0.40344 O] 

The cost function, defined by eqn. (2,41) will take the 


form as 
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J = j dt 

o ® 


(2.54) 


or 


where 


J 


4 - 

J _x g 2f, 
o 


Q = [T]"^ [t] 


To start the optimization routine, PSS obtained in Sec. 

2 .4. 2.1 case (a) is selected. The value of performance index 
obtained with this PSS is 0.88535. The optimum value of perfor- 
mance index is obtained as 0.26852 and the corresponding PSS 
transfer function and feedback gain matrix are obtained as : 


F(s) = 0»256 48 s^ 0.6 7 334 s 4- 0.48338 


s^ -i- 6.81412 s + 3.23453 


[K]: 


0.25643 

C 

1 


-0.34571 


0 


-1 .07434 


-3.23453 -6.81412 


The eigenvalues of the closed-loop system with the above PSS 
are given in the first column of Table 2.4. 
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Table 2.4 Closed-Loop ’Eigenvalues of the 
System 


.iith Parameter Optimization '.Jith Classical Control 

Technique Theory 


-0.48070 

-2.02667 + j 9. 90360 
-2.43789 + j 2. 37399 
-47.65066 


-3 .20609 

-3.57405 ± j 12. 64363 
-1.42775 + j 3. 05105 
-49.13654 


2.4.4 Design of PS3 with Classical Control Theory 

The performance of the stabilizers designed in Secs. 

2.4.1 to 2.4.3 and the stabilizer used, in power industry is 
compared by designing PSS using classical control approach of 
Larsen and Swann[26J • As per reconsaendations given in [26j» 
PSS utilizing speed signal should be designed at full load 
conditions and for strong systems (i.e., with low value of 
transmission line reactance). Thus, Xj = 0.2 p.u. is con- 
sidered. 

The design of PSS involves determination of gain and 
time constants of tvTO lead/lag stages in Fig. 2.4. Usually 
two identical lead/lag stages are used with T^ = and 

= T^. The ratio n.^. = T^/T 2 » known as band spread, varies 
from 2 to 10 depending upon the requirement of phase compen- 
sation. Thus, PSS design reduces to the date ii.dnat ion. of gain 
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Kg, band spread and the time constant Tj^. The values of 
these parameters are (see ;!^pendix E) • 

nt = 3 

= 0.5 sec. 

K = 0.05 

s 

and the corresponding PSS transfer function is 

p(s) = 0*45 s^ + 1. 8 s +1.8 
s^ + 12 .0 s + 36 .0 

Eigenvalues of the closed'-loop system with this PSS are given 
in the second column of Table 2.4. 

To compare the performance of stabilizers, time res- 
ponse of A6, and Aod are obtained for unit step change 

in the reference voltage when machine is equipped with the 
follovjing stabilizers ; 


(1) 

P3S-1 , 

obtained 

in 

Sec. 2.4.1 

(2) 

PSS-2 , 

obtained 

in 

Sec. 2.4. 2.1 case (a) 

(3) 

P3S— 3 , 

obtained 

in 

Sec . 2 .4 .3 , and 

(4) 

pDo— 4 , 

obtained 

in 

See. 2.4.4. 


The responses are obtained for full load and half full 
load conditions at 0.9 power factor lagging and are given in 
Figs. 2.6 to 2.8. 




FIG. 2.6 [RESPONSE OF 

(q) at full LO 






in rad /see 
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2.4.5 Discussion of Results 

The first design technique, based on modal control 
theory with partial state feedback is restricted, since PSS 
zeros have to be specified apriori and consequently full 
design freedom is not utilized. However, the techniqi;© is 
able to assign the dominant modes of rotor oscillation with- 
out significant shift of non-dominant eigenvalues. 

The second technique with output feedback is less 
restrictive as the number of closed-loop poles that can be 
assigned is more in comparison with the first design tech- 
nique. If zeros of PSS are specified, the sarae results are 
obtained (see first and third columns of Table 2.2 and PSS 
transfer functions given in Secs. 2.4,1 and 2, 4. 2.1 
case (b)). '.Vith two input compensator it is possible to 
assign all closed-loop eigenvalues. 

The third technique, based on parameter optimization, 
minimizes the fluctuations in power output of the generator. 
Responses of AP^, A6 and Aw of the power system, with the 
various PSS designed, are shown in Figs. 2,6 to 2.8. It 
should be noted that there is no justification for comparing 
the responses to judge the superiority of any one design over 
the others, since the design objective is different in 
different cases. However, we wish to draw attention to the 
aoprop ri atene ss of PSS— 3, if the main c^^^oern xs xeduction of 
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power fluctuations following system disturbance, as is clear 
from the responses of Ap^ shown in the figures. 

Responses of Ap^,A6 and Ato with 'FSB— 4 is better than 
those obtained with PSS—l and PSS— 2. But this does not reflect 
the superiority of the PSS designed using classical control 
theory over PSS designed using pole assignment techniques, 
because one can always assign poles to the same locations by 
first tm techniques as obtained mth PSS-4. The result 
shows the importance of selecting proper locations of closed” 
loop poles. 

2.5 CONCLUSIO?x!S 

In this chapter a basic system model, in the state-space 
form, of a single machine infinite bus power system, is deve- 
loped for the design of PSS. Three techniques for dynamic 
compensator design are presented. These techniques are 
successfully applied to a single machine infinite bus system 
example and the relevant numerical results are presented. 

Pole assignment is directly related to the design objective, 
of providing damping torque, used in power industries [6]. 
However, it may sometimes be difficult to choose proper pole 
locations which result in satisfactory system perforniance. 

The technique based on parameter optimization to minimize a 
given performance index has the advantage that it permits 
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minimization of power fluctuations at the output of the 
generator, following system disturbance. The method, thus, 
provides direct control on an important aspect of system 
performance. Indeed, the problem of power system dynamic 
stability is associated primarily with swinging of machine 
rotors and the consequent power fluctuations at the output of 
generators. 


CHAPTER 3 


INFLUEMCE CF GEMERATC.R MODELLING ON PSo DEoIGK 

3.1 INTRODUCTION 

In the previous chapter, a siniplified model of the 
synchronous machine, neglecting the effect of aniortisseur 
circuits, was considered for the design of Pod. Most of the 
literature available on the design of PSS utilizes the siinpli— 
fied model of synchronous machine neglecting the effect of 
damper windings, A detailed study on the adequacy of this 
simplified representation of synchronous machine, however, 
is not reported. The main objective of this chapter is to 
investigate the adequacy of the simplified machine model. 

A single machine infinite bus system is considered to 
analyse the effect of the machine modelling on the design of 
power system stabilizers. A sixth order model of synchronous 
machine [7S],with three amortisseur .circuits, is considered 
for the detailed representation of machine. A comparison of 
simplified and detailed representation of synchronous machine 
is m.ade by evaluating the effect of modelling details on 
(a) open-loop system behaviour and (b) characteristics of PSS 
designed for both representations of the machine. Aaequacy of 
the simplified representation of machine is analysed by employ- 
ing PSS, designed for lower order model, to the machine 
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represented in detail. Eigenvalue and time-domain analysis 
of the closed-loop system are carried out for this purpose. 

3.2 FOmiULATIOK OF SYSiei MODEL 


The single-machine infinite-bus system of Fig. 2,1 is 
considered. For detailed representation of synchronous 
machine, three amortisseur windings, as shown in Fin. D.2, 
are considered. The current source model of Ref. [ 73 ] with 
some modification, discussed in Appendix D, is utilized for 
the representation of synchronous machine. A single time 
constant excitation system of Fig. 2.3 is considered. Tur- 
bine and governor dynamics are ignored and it is assumed 
that the machine stator and the external network are in 
quasi-steady state. Dynamical equations of machine and 
excitation system (whose derivation is given in Appendix D) 
are as follows ; 


P“ = K ( - r 


P6 = 

P^d ^ ^l^d‘^^3^ f'^Sl^d'^S^fd 
Piq = ^k-^Cl2iq 

P^f ^5^d'^^7^f'’'^13^d’^^10^fd 






(3.1) 

(3.2) 


(3.3) 


(3.4) 


(3.5) 

(3.6) 


(3.7) 
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where C 2 , are functions of machine inductances 

and resistances; and are defined in Appendix D. Electrical 
torque, T^, is given as 

(3.8) 

where 

~ ^ ^q (3.9) 

si 1 - 

Eqn. (3.9) represents subtransient saliency. The equivalent 
circuit of the system is shovm in Fig. 3.1. The non-state 
variables and in eqns. (3.1) to (3.8) are eliminated 

by representing them in terms of state variables. 

A 

From Fig. 3.1, the armature current phasor i can be expressed 


A 


m 


3 


A 

I 


Zeq-^j^'d 


A 

V 


eq 


(3.10) 


A /V 

Substituting the values of I, v„„ and from Fig. 3,1 in 

’ eq eq ^ 

eqn, (3.10) and separating real and imaginary parts, we get 




Bry t a^ I Br 


= (^r r--Tiq;)^d - rriij i=i-j v., "-Vi"® 


+ Vgq(b 2 cos 6 -b^ sin 6 ) (3,12) 
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where 


aj_+Ja2 


k . Jli.. ™ 

R+ j (X+x^) 


. cosY+jsinv 
b^+jb^ = ~~ 

- R+J(X+x^) 

1 he terminal voltage can be expressed as : 
Vt = 


(3,13) 


Park’s components of terminal voltage can be obtained from the 
equivalent circuit of Fig. 3.1 as ; 


~ ^d^^d”^d^ (3.14) 


The state-space model is derived by linearizing eqns.(3.1) 
to (3.8). The non-state variables &i^, A i^ and A are 
eliminated using the following expressions which are obtained 
by linearizing eqns. (3.11), (3.12) and (3.13) 

= ^1 ^-^d + ^ 

Aid = Aid + 1% Alq + K^A6 
AV^ = AI^ -f Kg Alq T K^AS 


whore 


K 4 2 

*^1 - 


v i 1 
- i~?jr 


•^3 - e^-bgcos 6„) 

•i 


f 
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V ^ f ^ ^ 2 \ ^ ^ 2 

1% - - x;3]r 


% - ''eq[*( V^l T?53 l5^ri>=°® *0^ 

1 1 


,1 

A ^d 


'^- vf: t''dc5(i-5)+V'i-K4)3 


to 


v’* 

A ^d 


'% ^to ^ “’^qo 


and . x;| 

K^ = VT” ^^do 
to 


qo 6- 


The resulting seventh order model of the system is 
described by the state equations 


p X = [a] _x + b u 
whoro 


and 


x=[Aa) A6 AI^ ^%d^^ 


u = Au' 


(3 


The matrix A and the vector b are as follows 


[A] = 


b = 


^11 

®12 

®13 

®14 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

c 

®32 

^33 

®34 

®35 

0 

837 

0 

®42 

®43 

®44 

0 

®46 

0 

0 

®52 

®53 

®54 

^55 

0 

®57 

0 

®62 

®63 

®64 

0 

%6 

0 

_ 0 

®72 

®73 

®74 

0 

0 

877 

[o c 

G 0 

0 0 

by]* 





.16) 
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where 

^11 

®12 


-Kq/M 


^13 M ' ^~^qo‘*’Vqo"^l^do'*’^l ^ ^do^ 


*14" 'T^^do'^%^qo“^^do‘^^ ^ ^do^ 


^32 

= ^6^11 ^ 

®3o 

- C3+K4C 

®35 

= C3 

®37 

= s 

®42 

^3^12 ’ 

®43 

“ ^1^12 

^52 

— K 

” *^6^13 ' 

^53 

= C^+K^C 

^55 

= C^ ; 

®57 

= ^10 

®52 

= 

®63 

= K3C34 


*34 


%^11 


*72 


T 


Kq » 


R 


- ^R V 
^73 ■ " Tj^ ^ 


^74 


Tr. ^ 


•R 


®77 " “ T 


R 


Kr 


^7 “ 


R 
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3.3 EIGEImVALUE ANALYSIS 

In order to investigate the effect of modelling details on 
dynamic stability, eigenvalues of system matrix are calculated 
when the machine is represented by (i) a siimplified model and 
(ii) a detailed model. The simplified model is obtained by 
neglecting amortisseur windings. This results in a third order 
machine raodel as developed in Chapter 2. The detailed model 
considers three amortisseur windings and results in a sixth 
order model as discussed in previous section. For both models 
a single time constant excitation system, as shown in Fig. 2.3, 
is considered. 

3.3. 1 Numerical Example 

The single machine infinite bus system example of Chapter 2 
is considered for analysis, machine parameters, with usual 
notations in per unit on machine base, for the detailed repre- 
sentation are as follows : 


^d 

= 1.72 

f 


= 0.45 ; 

X" = 0.33 

x 

= 1.68 

m 

x’ 

= 0.59 t 

x« = 0.33 

q 


7 

q 


q 

^c 

= 0.33 





T ' 
^do 

= 6.3 : 

sec. 1 

T 1 

^do 

= 0.033 sec. 


T' 

= 0.43 

sec . ; 

T W 

= 0.033 sec. 


qo 


f 

C|0 




H = 4.0 sec. 
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is the axtnature leakage reactance. 

Excitation System Parameters ; = 50 j “ 0.02 sec. 

System Parameters and Operating Data (in p.u.) : 

Transmission line (Fig. 2.1) : = 0.024 j = 0.6 

(on machine base) ^ ^ 

Ge = 0. ; = 0.065 

Infinite bus voltage =1.0 /0° 

Generation = 1.0 p.u. at 0.9 lagging power factor (i.e. 

3^ = 0.9 + jO. 43589 ) • 

Data of Equivalent Circuit : R = 0.026C2 ; X = 0.62470 

V__ = 1.04123 /--0.0945° 
eq ^ — 

Eigenvalues of the system matrix, for both lower and 
higher order models, are computed for variation in the follow- 
ing three parameters of the system i 

(1) Machine Loading - Machine loading is varied from 0.25 p.u. 

to 1.0 p.u. at constant power factor of 
0.9 lagging. 

(2) Inertia constant - Following three values of H are taken ; 

(a) 2.0 sec., (b) 4.G sec, (nominal 
value) and (c) 5.0 sec. 

(3) Voltage regulator — Following thr®© valijes -of are taken; 
gain 

(a) 25, (b) 50 (nominal value) and 

(c) 100. 
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Loci of the complex eigenvalue corresponding to the rotor 
oscillation mode are given in Fig. 3.2. Remaining eigenvalues 
of the system are given in Tables 3.1 to 3.3 for variations in 
the machine loading, inertia constant and the voltage regulator 
gain. 

Table 3.1 Open-Loop Real Eigenvalues for Variation 
in Machine Loading 


















Table 3.3 Open-loop Real Eigenvalues for Variation 
in Voltage Regulator Gain 















Simplified model 
Detailed model 




Variation in voltag® regulator gain. 
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3.3.2 Discussion of Results 

Inspection of Fig. 3.2 shows that the nature of eigen- 
value loci is same in both representations of the raachine for 
the variations in the parameters considered. More damping of 
rotor oscillation mode is observed in the detailed represen- 
tation. This is to be expected because of the inclusion of 
amortisseur windings in the detailed model. 

The effect of increasing the machine loading is to reduce 
the damping of the rotor oscillation mode. Reducing the iner- 
tia of the machine or increasing the voltage regulator gain 
have similar effects. 

The frequency of the rotor oscillation mode is related to 
the synchronizing torque coefficient and the inertia, the 
frequency increases with reduction of inertia. Increasing the 
voltage regulator gain has the effect of increasing the syn- 
chronizing torque coefficient thereby increasing the frequency 
particularly at higher values of gain [6]. Increasing the 
loading on the machine has not much effect on the synchronizing 
torque coefficient because the increase in the internal vol- 
tage of the generator counteracts the effect of rotor angle. 

It is observed that the synchronizing torque coefficient is 
less with the simplified model. This is because of the fact 
that the synchronizing torque contributed by the amortisseur 
windings is neglected. This contribution is seemed to be 
significant particularly at light loads- 
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Non-oscillatory modes are not of much concern as long as 
they are sufficiently away in the left half of the complex 
plane and their variations are, therefore, not important. How- 
ever, an inspection of Tables 3.1 to 3.3 shows that they are 
not affected much, in both representations, by load variation 
and variation in inertia constant. They vary with the varia- 
tion in the voltage regulator gain as shown in Table 3.3. 

High voltage regulator gain in the detailed representation of 
machine shows an additional oscillatory mode. As this mode 
is highly damped, it will not effect much the stability of the 
system. 

3.4 PSS DESIGN 

To investigate the effect of modelling detail on the PSS 
characteristics, PSS are designed for both the representations 
of machine. Pole assignment technique, given in Sec. 2. 4. 2.1, 
is used for the design of PSS utilizing rotor speed as feedback 
signal. Variations in the following throe parameters of the 
system are considered for PSS design ; 

(1) Machine Loading - 1.0, C.5 and 0.25P.U. loading at 

0.9 lagging power factor 

(2) Inertia constant - 2.0, 4.0 and 6.0 sec. 

(3) Voltage regulator- 25 , 50 and 100. 
gain 
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A second order PSS with transfer function of the form 

e„(s+z^)(s+Z2) 

rs+FpTsTSJi 

is designed for both representations of machine. The objective 
is same for each P3S design, namely, to assign five closed-loop 
poles to the specified locations. These are -l.C+j7.0,-l.C, 
-2.0 and -6.0. 

PSS parameters obtained for different machine loadings and 
for different values of inertia constant and voltage regulator 
gain are given respectively in Tables 3.4, 3.5 and 3.6. Fre- 
quency plots of the PSS transfer function for each case are 
given in Figs. 3.3 to 3.5. 


Table 3.4 PSS Parameters for Variation in Machine Loading 


Machine 

Loading 

Model 

Order 

®o 

^1 

^2 

P- . — — 

Pi 

L ... 

P2 

• 

o 

Lower 

Higher 

0.06501 

0.11148 

-5.86475 

-7.67585 

3.22300 

3.17181 


-0.50496 

-0.65124 

, 0.50 

Lower 

0.11286 

-6.06547 

0,47719 

-5.80205 

-0.69254 


Higher 

0.05374 

-6.21340 

1.45426 

-7.91252 

-0.78866 

0.25 

Lower 

Higher 

0.22691 

0.13097 

-6.10750 

-10.71719 

-0.61723 

-1.57561 ' 

-5.46609 

-20.02707 

* 

-0.87463 

-1.06966 





















85 


Table 3,5 PSS Parameters for Variation in Inertia 
Constant 


Inertia 

constant 

iModel 

Order 

®o 

^1 

^2 

Pi 

^2 

2.0 

Lower 

Higher 

0.03197 

-0.05000 

-11.27945 

- 7.99723 

-6.18399 

-3.49388 

-5.77025 

-4.36127 

-0.27530 

-0.40832 


Lower 

0.06501 

-5.86475 

3.223CC 

-6.41258 


4.0 








Higher 

C. 11148 

-7.67585 

3.17181 

-21.02139 

Hum 


Lower ! 

0.15799 

-5.88597 

-0.15379 

-6.41940 

-0.76093 

6.0 


I 

1 






Higher ; 

0.17414 

-6.38861 

-C. 09334 

-11.02523 

-0.81156 


Table 3,6 PSS Parameters for Variation in Voltage 
Regulator Gain 


V.R. 

Gain 

Model 

Order 

o 

CD 

"l 

^2 

Pi 

P2 

25 

Lower 

Higher 

0.18051 

0.16373 

-3.25016 

-2.34043 

C. 65649 

1.31371 

-8.50293 

-10.53550 

-0.77156 

-0.82940 

• 

bo 

Lower 

Higher 

0,06501 

0.11148 


3.22300 

3.17181 

■ 

-0.50496 

-0.65124 

100 

Lower j 

Higher 

0.04863 

1 

0.01952 

-C. 21972 

-7.97605 

i jl. 97548 j 

9.07333 

-1.03975 

-6.82988 

+30.62952 

-C. 5270 6 
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3.4.1 Discussion of Results 

An inspection of the gain and phase characteristics of 
stabilizers, given in Figs. 3.2 to 3.5, shows that in most 
of the cases they are similar for the simplified and the 
detailed representations of machine except for the case when 
the voltage regulator gain is equal to 100. The typical gain 
and phase characteristics of PSS transfer function for 
~ 100, with the simplified machine model, shown in 
Fig. 3.5, is due to the presence of a complex pair of zeros 
with the oscillation frequency of about 2 rad. /sec. 

It is interesting to observe that the PSS characteristic 
vary widely for the different cases considered. This is due 
to the fact that no limitations regarding the locations of 
zeros and poles of the PSS are imposed in the design, in 
contrast with the classical design of PSS where only lead 
compensators are considered. This does not provide full 
freedom that can be utilized to assign as many closed-loop 
poles as possible. There is no need to restrict the locations 
of zeros and poles, except for the requirement that the poles 
be in the left half plane. In all the cases of PSS design 
given here, poles of the PSS are in the left half plane. 

Bode plots of PSS transfer function, given in Figs, 3.3 
to 3.5, show that in most of the cases, the stabilizer has non 
minimum phase characteristics. This is a desirable feature 
which reduces the effect of torsional interaction [26 J . 
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3.5 ADEQUACY OF LCV/ER ORDER REPRESENTATION OF aACHINE 

Adequacy of lower order representation of machine for PoS 
design is analysed in this section. This is carried out by 
designing PSS for both the representations of machine, with 
the same objective of placing the five closed-loop eigenvalues 
at -1.0+j7.0, -1.0, -2.0 and -6.0. PSS^and PSSj^ are res- 
pectively the stabilizers obtained for the lower and the 
higher order representation of machine. Closed-loop eigen- 
values of the system, with detailed representation of the 
machine, is then obtained with PSS^ and PSSj. The analysis is 
carried out for three different loading conditions, namely, 

1,0, 0.5 and 0.25 oer unit at 0.9 lagging power factor. 
Eigenvalues of the closed-loop system, for each case, are 
given in Table 3.7. 

The performances of the two stabilizers and PSE^ , 

designed for full load conditions, are compared by observing 
the responses of S and m for unit step change in the 
reference voltage. The responses at full load and at half 
full load are shov/n in Figs. 3.6 and 3.7. 

3.5.1 Discussion of Results 

An inspection of Table 3.7 reveals that the damping of 
rotor oscillation mode, obtained v/ith stabilizer designed for 
simplified machine model {PSS^^) is comparable with that 
obtained with PSS designed for higher order machine model (PSS^^) . 
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Table 3.7 Closed-loop Eigenvalues with Detailed 
Representation of Machine 


^?i(Iachine 

Loading 

1.00 ( p . u . ) 

0.50 ( p . u . ) 

0.25 ( p . u . ) 






— 1*0 ^ j7*0 

-1.0 + j7.0 

-l.C + j7.0 


-1.0 

-1.0 

-l.C 


-2.0 

-2.0 

-2.0 


-6.0 

-6.0 

-6.0 


-14.023724 

-6.806135 

-16.827955 


-31.380179 

-28.797788 

-26.791239 


-47,080218 

-41.306949 

-46.949659 


-54.760952 

-56.005314 

-56.100309 

1 

-0.687814 + 
j 7. 628617 

-0.718301 + 
j 8. 130632 

-0.716801 -r 
j8. 395452 ” 


-1.140974 + 
j 0.314649 

-1.081705 + 
jO. 304101 

-1.100936 + 
jC. 181090 "■ 


-5.C503C9 

-4.547683 

-4.273976 


-5.552861 

-6.111549 

-6.1S3122 


-26.822258 

j 

-26.035408 

-29.911529 


-47.055388 i 

-46.980672 i 

-46.959130 


-55.351604 

-55.791721 

-55,949948 




rad /see. , ir 


3 ? 



FIG. 3.6 RESPONSE OF A 8 AND Aoo AT FULL LOAD 
FOR AVr^i rl.O p.u. 



degrees 


13 


0 With PS5 designed Tcdp! 
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ihis is true for all the three loading conditions 
considered. It is observed that with the application of PSS^^ 
another oscillatory mode of low frequency and high damping 
ratio is introduced v/hich will not affect much the stability 
of the system. 

Responses given in Figs. 3.6 and 3.7 show that, for the 
disturbance considered, performances , of the closed-loop system, 
when machine is equipped with P53^ and PSS^ are comparable. 

The initial peak and settling time obtained with PSS^ a^re less 
than that obtained with PSS^^, however, the difference is not 
significant. 

Thus from eigenvalue analysis of the closed-loop system 
and time responses given in Figs. 3.6 and 3.7, it can be 
concluded that simplified representation of machine for P33 
design is adequate. 

3.6 CONCLUSIONS 

A detailed investigation into the effect of generator 
modelling on PS3 design has been presented in this chapter. 

A sixth order machine model including representation of amor- 
tisseur circuits is compared to the third order machine model, 
described in Chapter 2, in terms of the system behaviour and 
the characteristics of P3S. 



93 


The results from a numerical example indicate that the 
open-loop system behaviour is predicted from the simplified 
and detailed machine models are similar. The characteristics 
of the PoS designed with the same objective are also similar 
in most of the cases. 

Adequacy of the simplified machine model is analysed by 
incorporating stabilizer, designed for simplified model to the 
machine represented in detail. Eigenvalue and time-domain 
analyses of the closed-loop system are used to demonstrate the 
validity of using simplified machine model in PSS design. 



CHAPTER 4 


INVESTIGATION OF CONTROL SIGNALS 

4 . 1 INTRODUCTION 

For improving dynamic stability through PSS , it is 
important to select a proper feedback signal. Since the 
purpose of the stabilizing signal is to damp— out electro- 
mechanical oscillation? rotor speed, electrical po’wer and 
frequency deviation have been considered as input signals to PSS 
[6, 21-27]. Other control signals derived from terminal 
voltage, armature current, active and reactive powers and 
field current were investigated [21] and it was found that 
none of these signals is as effective as signal derived from 
generator speed. Larsen and Swann [26], in their three part 
paper, analysed the effectiveness and tuning of stabilizers 
using speed, pov/er and frequency signals. The analysis 
presented in Ref. [26] is based on gain and phase characteri- 
stics of the transfer function of the plant consisting of 
the generator, the excitation system and the power system. 

The plant transfer function is strongly influenced by vol- 
tage regulator gain, generator power level and ac system 
strength. They recommended that speed and power stabilizers 
should be tuned for strong ac systems and stabilizers using 
frequency signal should be tuned for weak ac systems. These 
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recoiTOnendations are based on classical design of PSS. The 
design approach in this thesis is based on modem control 
theory as given in Chapter 2, where complete design freedcro 
is utilized. The validity of the recomiTiendations given in [26] 
for PSS design in general needs to be examined. 

Interaction between two or more parallel connected gene- 
rators in a generating plant [25, 3l] makes the choice of con- 
trol signal difficult. A detailed study in this regard is not 
available in the literature. 

The principal aim of this chapter is to present a case 
study on the effectiveness of various feedback control signals 
that can be used to design PSS for improvement of dynamic sta- 
bility. Three signals, namely, rotor speed, electrical power 
output and terminal bus frequency are considered. A generat- 
ing plant, having a number of generating units, connected to a 
large power system is considered for analysis. First, the 
effectiveness of control signals is obtained for a plant having 
one machine, that is, for the single machine infinite bus . 
system described in Chapter 2. Influence of interactions between 

machines present in the plant on the effectiveness of control 
signals is then investigated by considering anotter machine operat- 
ing in parallel with the existing machine. A general multi- 
machine system model is developed for this purpose in which 
dynamics of the generators are described, by linearized ^tate 
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equations and the network is represented by the reduced 
admittance matrix which retains only internal buses of 
the generators* The non— state variables in the system model 
are eliminated by the use of power flow equations. The 
important feature of this formulation is that the elements of 
system matrices are directly evaluated from the knowledge 
of generator loading and reduced bus admittance matrix. 
Although the system model developed here is applied to a 
simple system of a generating plant, having a number of 
machines, connected to a large power system , it can be used 
for any multimachine power system v/here the network is 
represented by the reduced admittance matrix retaining 
generator internal buses. 

PSS are provided to achieve increased damping of electro- 
mechanical oscillations. This enables enhancement of the 
dynamic stability limit. It is necessary to evolve a cri- 
terion to judge the effectiveness of PSS, designed for 
different feedback signals. A convenient criterion is to 
base the judgement on the stable region of operation in the 
P-Q plane for the machine equipped with the PSS. The chapter 
also includes study of the influence of system strength on 
the design and perforipance of PSS. The suitability of average 
speed and average power signals- has also been s-tudied in the 
case of plant with more than one machine. 
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4.2 SINGLE MACHINE INFINITE BUS SYSTB4 

The single machine infinite-bus system of Chapter 2 is 
considered. For simplicity, transient saliency of generators, 
resistance and line charging of transmission line are neglected, 
that is, it is assumed that = x^* , = 0 and Yg = 0. Other 

parameters of machine and excitation system are the same as 
given in Chapter 2. 

4.2,1 Design of PSS 

Effectiveness of the three control signals, namely, rotor 
speed, electrical power output and bus frequency are obtained 
with the designed PSS. Stabilizers are designed for full load 
conditions, that is, machine delivering rated power at 0,9 
power factor lagging (P+jQ = 0.90CX3+j0.4359 p.u.). The follow- 
ing three cases of system strength are considered in the de- 
sign of PSS for each control signal : 

(1) X^=0.4p,u. (Strong system), 

(2) Xj = 0.6 p.u. (System of moderate strength), and 

(3) X^ = 0,8 p.u. (Weak system). 

System matrices obtained for X^ = 0.6 p.u. are given 

0 

0 I 

0.35273 
-50,0000 — ^ 


below ; 


[A] = 


0 
1 
0 
»- o 


-35.36908 

0 

- 0.30157 
521.57320 


-11.89618 

O 

- 0.35072 
-618.42114 



ICO 


b = [0 0 0 2500]^ 

£ = 0 0 O] for speed signal 

= [O 0,90067 0,30293 O] for power signal 

= [0,66252 -0,01909 -0,02220 0,02233] for frequency 

signal 

Open-loop eigenvalues for the above three cases are given in 
Table 4,1. 


Table 4,1 Open-loop Eigenvalues 



= 0,6 p.u. 

Xj = 0.8 p.u. 

0. 26105+ j 7. 21497 

0. 35281+ j6. 19653 

0. 46974+ j 5. 06393 

-4.86833 

-5.95422 

-6.79143 

-46.04967 

-45.10213 

-44.46803 


Eigenvalue assignment with output feedback as given in 
Section 2. 4, 2.1 is used to design stabilizers with each of 
the three control signals and for each case of system strength. 
As discussed in Chapter 2, with a second order PSSj at the 
most five closed-loop eigenvalues can be assigned. The loca- 
tion of the remaining closed-loop eigenvalue is obtained 
alongwith the paremieters of PSS. The desired locations of 
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closed-loop eigenvalues and location of the resriaining eigen- 
value obtained for each case are given in Table 4.2. The 
parameters of PSS vvith transfer function of the form 

_ 8o(srZj^)(s+Z2) 

F(s) 

are obtained using each of the three control signals and are 
given in Table 4.3, where and Z2 are the zeros and and 
P2 are the poles of PSS. 

Frequency response of PSS transfer function, designed at 
= 0.6 p.u. is given in Fig. 4.1 

Table 4.2 Closed-Loop Eigenvalues 



Signal 

• Ui • 

Xj=0.6 p.u. 

Xj == 0*8 ] 



- 2.14359 

- 1.87364 

- 1.60770 



+j8.0CXX)0 

+j7. 00000 

+36.00000 

Desired 


- 4.86833 

- 5.95422 

- 6.79143 

Location 


- 1.00000 

- 1.00000 

- l.OOOCXD 



- 2.(XXXXD 

- 2.COOOO 

- 2.00000 


Speed 

-46.90553 . 

-45.89207 

-45.21988 

Remaining 

P ower 

-42 .02226 

-41.65311 

-41.86691 

Eigen- 





values 

Freq. 

-35.87726 

-35.99440 

-35.44419 



102 


Table 4,3 Parameters of P33 


Signal 


©O 



Pi 

P2 


0.4 

0.12954 

7.83588 

-4.36832 

-3.34385 

*""0 #32131 

Speed 

0.6 

0.12084 

5.25912 

-5.95422 

-7.89546 

-0.35140 


0.8 

0.11513 

2.78286 

-6.79144 

-7.50220 

-0.40442 


0.4 

-0.69783 

-5.91149 

-4.86835 

-2 .97439 

-0.80750 

P owe r 

0.6 

-0.59832 

-6.43534 

-5.95450 

-3 .22773 

-0.78017 


0.8 

-0.45801 

-8.18041 

-6.79142 

-3.81797 

-0.73576 


0.4 

0.16589 

9 ,08061 

-4.86833 

-4.95754 

-0.41363 

Frequency 

0.6 

0.13263 

6.00710 

-5.95423 

-5.34637 

-0.40702 


o 

• 

OD 

0.11228 

3.22053 

-6.79141 

-5.51648 

-0.43110 


4.2.2 Determination of Stable P-Q Region 

With each PSS of Table 4.3, stable P* Q region is obtained 
by varying the active power output in steps and determining the 
limiting values of reactive power output for closed-loop system 
stability. To study the effect of ac system strength on the 
performance of stabilizers, stable P-Q regions are obtained 
with each PSS of Table 4.3, using appropriate signal, at 
different values of transmission line reactances, that is, at 
equal to 0.2, 0.4, 0.6 and 0.8 p.u. These regions are 
shown in Figs. 4.2 to 4,5. The stable F-S region in these 
figures is the region above the appropriate boundaary lines. 
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Figure 4.6 compares the stable P-Q regions? for speed 
signal feedback, obtained with PSS designed using classical 
control theory (sec 2.4.4) and that obtained ¥/ith PSS of 
Table 4,3, at different ac system strengths (i.o. at Xj equal 
to 0.2, 0.4 and 0.5 p.u.). In case of PSS of Table 4.3, the 
PSS which gives maximum P-Q region at a particular ac system 
strength is chosen for comparison, 

4.2.3 Discussion of Results 

From Fig. 4.1 it is clear that the gain required with 
power signal is higher than those for speed and frequency 
signals. Speed signal requires minimum gain. In the fre- 
quency response plot it is seen that gain requiren^jnt reduces 
with increase in frequency for all the three signals. Phase 
characteristics of speed and frequency signals are the same. 
Phase requirement reduces to zero for higher frequencies in 
case of speed and frequency signals. For power signal the 
minimum phase is required at 0.5 Hz. 

Sys tem Stre ngth for PSS De sign 

Speed Signal : 

From Figs. 4.2 to 4.5, it is evident that PSS utilizing 
speed signal should be designed for weak system as the stable 
P"Q region obtained is maximuBa with such a PSS designed at 
= 0.8 p.u. Difference in stable r^egioas attained is 
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significant for strong system, (see curves i(a), 1(b) and 
l(c) of t-ig . 4.2). The difference reduces with the reduction 
in system strength, (see curves 1(a), 1(b) and 1(c) of 
Figs . 4.3 to 4.5). 

power Signal : 

Curves 2(a), 2(b) and 2(c) of Figs. 4.2 to 4.5 show that 
the stable P-Q regions obtained utilizing povier signal, with 
PSS designed at X-j. equal to 0.4, 0.6 and 0.8 p.u. coincide, 
suggesting that PSS can be designed for either v/eak or strong 
system. 

Frequency Signal : 

In case of frequency signal the following typical 
behaviour is observed in curves 3(a), 3(b) and 3(c) of 
Fig. 4.2 : for strong system, i.e., at Xj = 0.2 p.u. stable 
P'-Q region obtained with PSS designed at Xj = 0.6 p.u. is 
more than that obtained with PSS designed at Aj = 0.4 p.u* 
but again it reduces with PSS designed at = 0.3 p.u. For 
system of moderate strength and for weak system, stable P-g 
region obtained with PSS designed at Xj = 0.3 p.u, is maximum, 
curves 3(a), 3(b) and 3(c) of Figs. 4.3 to 4.5. Thus, for 
system of moderate strength and for weak system, it appears 
that PSS utilizing frequency signal should be designed for 
weak system. 
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Compa rison of Sign als 

From i-igs. 4.2 to 4,5 it is clear that stabilizers 
utilizing power signal give maximum stable P—Q region in 
comparison with PSS utilizing speed or frequency signal for 
all cases, irrespective of the system strength at v/hich it is 
designed and irrespective of the system strength at which it is 
operating. In fact, with pov/er signal, network limit is reached 
first for all values of Xj except for = 0.2 p.u. where 
instability is observed before network, limit. 3vsn for this 
case, = 0.2 p.u., the stable P~Q region is very close to 
network limit as shown in Fig. 4.2. For stronger system, 

Xj = 0.2 and 0.4 p.u., the difference in stable P-Q region 
obtained with power signal and that obtained with speed or 
frequency signal is significant. The difference reduces with 
the reduction in system strength. Figure 4.2 shows that 
speed signal is most inefficient for strong systems for all 
loading conditions. Figures 4.4 and 4.5 show that for v#eak 
systems and at higher loading condition speed signal gives 
more stable P-Q region than frequency signal but the difference 
is not significant. 

From Fig, 4,6 it is observed that PSS designed by classi- 
cal control theory utilizing speed signal gives larger stable 
P—Q region than that obtained by PSS of Table 4.3 for strong 
system, curves 1(a) and l(a')* For w®ak system it is vice-versa 
( see curves 1(b), iCc) , l(b* ) and l(c’)/* 
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4.3 GHNHRATIMa PLANT INFINIT3 BUS SYSTai 

In order to investigate the interaction of other ’■."achines 
present in the plant, on the suitability of control sirnals, 
another machine is considered operating in parallel ■■vith the 
existing machine of section 4,2. general system model is 
developed as follows : 


4.3.1 Development of System liodel 

Consider a generating plant having n number of generating 

units supplying power to an infinite bus through a transmission 

line of impedance as shown in Fig. 4,7. represents line 

charging and local loads of the generating plant, if any. Fig. 

4.8 is the equivalent circuit of the systera shown in Fig. 4.7, 

where Z and V /y are the Thevenin* s equivalent of the 
eq eq ^ 

transmission system and can be expressed as 


A 

e q (T+Yg Zj') 


(4.1) 


’^eq L-i = rr+^T 


v/o 


(4,2) 


H. /6- for i = 1,2, ...» n are the voltages of internal buses 
of alternators. 

4. 3. 1.1 Development of Dynamical Equation 

The following assuaptions are made in the development of 


system model. 




FIG. 4.8 EQUIVALENT CIRCUIT OF SYSTEM 
GIVEN IN FIG. 4.7. 
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(1) Synchronous machines have no amortisseur winding 

(2) Transient saliency is ignored for all machines 

(3) Governor and Prime Mover dynamics are ignored for all 
machines . 


The system model is developed by combining machine and 
network models. The network is represented by linear algebraic 
equations utilizing the reduced bus admittance matrix which 
retains only the internal buses of machines. 2ach machine of 
the plant is represented by the circuit model given in Fig. 4.9 
(Appendix D) . Parks components of variables shown in Fig. 4.9 
can be given from Eqn. (D.16) of /^pendix D as 


A 

I 


(I* 

q 




A 

i 


m 


= ( 




q 


jl(j) eJ® 

3Hd) 


A 

V 


m 


= (v. 


" jVn) 


where 6 represents the angular displacement of machine q-axis 
with respect to reference frame (infinite bus). From Figs. 
4.9(a) and 4.9(b) v/e have 

A ^ 

E = jx^ I 

Separating real and imaginary parts, ^ get 

E‘ = Xj I - 
q d a 


^d 


(4.3) 
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Since transient saliency is neglected, I' v;ill be zero as 

M 

Iq = §iq (Appendix D) when § = 1-x^/x^. Thus, there ’.vill be 
no voltage source resulting from I’ and the voltage source ^ 

M 

will be given as 



(4.4) 


From Eqn. (4,4) and Fig. 4,8, the voltage of internal bus of 
ith machine can be expressed as 




(4.5) 


The linearized dynamical ecjuations of the ith machine and 
the excitation system can be given as 


“i = IT 5^” ^^ei " ^ ^“i^ 

pA6j^= (4.7) 

pAE^ = ” [ AE^^^ - fiE^ - ^^di^ 

GDI 

Efdl = k- [%( ] (4.9) 

Rr 3.. 

The non-state variables ^Pg-i » ^^di ^^Ti ecjns. 

(4.6) to (4.9) can be eliminated by expressing them in teras 
of state variables as follows. 

Representation of A P^ 

The active and reactive power output of the machines at 
their internal buses. Fig. 4.8, can ba expressed as 
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n+1 

P. = E. 2 E, C. . 
^ ^ j=l J 13 


( 4 . 10 ) 


n+1 

= "i 3=1 °ij 

for i = 1,2, n * 


( 4 . 11 ) 


where 


cos( 6^- 6j) + Bj^j sin( 6 6j) 


°ij = °ij «i- '5j) - Bij «1- 6j) 


Vij S G,j ^ 


and 


^n+1 -^*^0+1 ^eq^^ 


Yij is the ijth element of reduced bus admittance matrix and 
can be expressed as follows. The star point at the terminal 
bus. Fig. 4.8, is eliminated by converting the star into an 
equivalent delta from the relation 

n+1 , 


= Z.Z. 2 j— 
^ J K=1 


(4.12) 


for i and j = 1,2, n+1. in Eqn. (4.12) is the impe- 

dance of ith machine, Z^^^^ is equal to Z^^ and 2 1“ is the sum 
of reciprocals of all impedances connected to the terminal bus 
including Z^^. The element Yy of the reduced bus admittance 
matrix is then obtained as 
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Y. 


IJ 


JL 

z. 


IJ 


for i j 


and 


Y.. = 2 


11 


^IK 


where subscript K represents the bus numbers connected to ith 
bus. Linearizing eqns. (4.10) and (4,11), we get 


-2 X A ^ _ .,2 V 1 


APi = .Ej^) ^6 i+{Pi„+G.,E^ ) 


10 11 lo'^ E« ^ 1 


10 


n 


^1° j=i 


(4.13) 


and 


AQ. 

1 


(P. -G-,Ef^) A6. + (Q. -B..E? ) i— AE 
^ 10 11 lo'^ 1 '^^lo 11 10^ E. 


10 


n 




(4.14) 


Representation of A i^ 

Referring to Fig. 4, 9(b), power output at the internal 
bus of ith machine can be given as 


A Ai^ 

= Pi+jQi = E.. i 


m 


(4.15) 


Separating real and imaginary parts from eqn, (4.13) and re- 
arranging, we get 




^di ~ E 


1 

P. 


1 . = 
qi E 


and 


(4,16) 
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Linosrizing egn* (4»16) and substituting th© ¥alije of 

x 

and AQ^ from eqns. (4.13) and (4.14), we get 

= (Pio-'^hEL* - Bu iEj 


10 


n 


and 


*V°ij 


*iql = l-ae.+G.. AE. 


(4.17) 


lO 


n 


Representation of AV^ 


* jfl ‘Ej) 


(4.18) 


The magnitude o'f terminal voltage can be expressed as 
^Ti “ ^"^di (4.19) 


Linearizing eqn. (4.19), we get 


Av_. = TT — (vj. AVj.+v . Av .) 
Ti ^Tio 


(4.20) 


Now, from Fig. 4.9(b), the Park's component of terminal voltage 
can be expressed as 


Vj- 

di 


and 


V 


qr 


x’ i . 
^di ^qi 


^di 


(4.21) 


(4.22) 
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Linearizing eqns. (4.21) and (4.22) and substituting Ai^. 
and Ai^. from eqns. (4.17) and (4.18), we get 


A6.^=C. , G,, AE. 


Y * 

^di 


n 


’'ii j‘i ''Ej) 




and 


XL. 


Av- = - ^ eL) 


qi 


(4.23) 


n 


’‘dl 

Vi 


(4.24) 


From eqns. (4.20) to (4.24), we get 


= AE.+ 2 (pjj A6.+p|^ AEJ (4.25) 


where 


r 

ti 

= A 6 

ii 

X* 

^di 

1 

Vt-. 

Tio 

ii 

1 

2 

Vt- 

Tio 


= - 


P 




CVi(hoBii+«io°ii) + (Pio-°ii4n 


Ct- '"di ‘'“lo°ii-%o®ii)-<5io»Eio(l+Xi.B^.)] 


EE ■ ^ ^ 

[ di (P.^D. .-Q^^C. .)+x^^C. -1 


^Tio E 


lO 


Po-^ 


XV? 


Tio E. 


lO 


and 
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and 


where 


Si 


^2 = 

= 

®13 

a^^ — 
®32 


[0 0 0 b^]* 


-Kdi/Mj 




-<ho+°li ^io’ 

-(x^i-x^l) (Pio-3iiE2^)/{Ti^. Ej^) 


^33 = B„-1]/T 


di '^di' “ii 


doi 


^34 


= 1/T 


cioi 


®42 “ ^Tio^ 


...2 


^43 


= -^Rit^di ^"di^^ioGir^ioBii)“'^ioJ/Eio^HioU+^kSii)3/ 


^Tio> 


®44 “ 

for j 7^ i 

aij - 
^12 ~ 

aij _ 
®13 ~ 

= 

®32 


-1/T 


Ri 


-E. E. D. ./M. 

10 JO ij^ 1 


■EipCij/«i 


(x ,--xi.)E. C. ./T’ . 
'“^di di'' JO ij^ doi 
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)D,j/T 


f 

dio 


= -KRiEioEjo 


ij"' io' 




and 

4, 3. 1.2 Development of Output Equation 

The output equation, in general, can be written as 


X = [C] X 


where the elements of matrix C depend upon the output signal. 
Here, three signals are considered, namely, speed, power and 
frequency and matrix C is defined for each case. In case of 
speed and power signals, individual value for each machine and 
average value of all machines are considered. 

Speed Signal 

Change in rotor speed of machines is utilized as control 
signal and can be expressed as 

(a) Individual Speed 




P Ol 

C2 


*" ^1 

Xg 

• 

• 


0 ■ c 


• 

Am 

— , . n,-J 




x„ 

u a —i 


1 


(4.27) 
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where C^ = [l 0 0 O] 

and X. = [Aw Aa. AE. AE-..]^ 

1 111 fdi^ 

(b) Average Speed i 



where ~ n ^ ^ 

Power Sign al 

As the armature resistance of machines are neglected, 
active power available at the terminal bus of a generator 
will be same as the active power developed by the generator 
at its internal bus and is given by eqn. (4.13) as shown 
below. 



and 
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(b) Average Power 


Feedback signal for each machine will be 


n 


7=^2 AP, 


i=l 


Freouency Signal 

Frequency signal can be derived from the voltage 
the terminal bus as 


Af. = pA8i 

jl. a. 

The voltage angle 0^ can be expressed as (Appendix D) 

0. = 6.-6. 

1 i *^1 • 

where Pi tan“^ (r™) 

qi 

Linearizing eqn. (4,30) and substituting the values of v^^, 
V -5 Av,. and Av . respectively from eqn. (4.21), (4.22), 

C|JL ClX 

(4.23) and (4.24), we get . 


^®i " ■*^1^ A6 ^+tiJ^ AE^+ A5^+t)|^ AE.) (4.31) 




where 


'll" = :V t'^lo-^dl *E— (Pfo-^«L>+='di'Plo'^ii-''ioBil) 

. V .x*" « 1 o 

Tio 


„ii - 
^2 


. E . 
Tio lo 
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ij 

1 


and 





E. E. 

9 

Tio 





E. 

lO 

V^. 

1 lO 


x‘2 

■ di 


iO 


(Q. C. .-P. D. .). 
ij ^lo 


■xL-C. .1 
di ijJ 


Substituting the value of from oqn. (4.31) in oqn. (4.29). 
we get 


n 

Af. =2 C. X. (4.32) 

1 ~i --x 


v/here is a row— vector and can be expressed as 




n 

i 7 ] a^o 

j=l 2 


n 

2 

j=l 






3 


4, 3. 1.3 Calculation of Operating Point 

The development of system model presented in the previous 
section is general and can be applied to any multimachine 
system where the network is represented by reduced admittance 
matrix retaining generator internal buses. Elements of 
system matrix can be directly evaluated, as shown, from the 
knowledge of active and reactive powers and voltages at the 
internal buses of generators. For any multimachine system 
this information can be obtained from load flow analysis. For 
the system considered. Fig. 4.7, bus quantities, i.e., P,Q,E 
and 6 at the internal buses of generators can be obtained 


as follows. 
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For a given plant output and fixed infinite bus voltage, 
first the magnitude and angle of terminal bus voltage are 
obtained and then from the knowledge of load distribution 
factor, internal bus quantities are obtained. From Fig, 4.8, 
plant output can be expressed as 


where 


Pj+jQj = 





(4.33) 


(4.34) 


Substituting the values of i, and /y respectively 

from eqns. (4.34), (4.19) and (4.2) in eqn. (4,33) and separa- 
ting the real and imaginary parts, we get 


where 

and 


Vj V cos® = a-u 

V sine = b-v 

If- 

U + jv = (Pj + 

a + jb = (1 -4- Y^) 


(4.35) 

(4.36) 


Solving eqns. (4,35) and (4.36) for Vj and 6 , as 

Vt = [ fe - '-‘>1 

^ V&- V 

9 = tan ) 

YjB - u 


and 


(4.38) 



128 


2 2 

v/here aj^ = a + b 

= -(V^ + 2ua + 2vb) 

2 2 

= u + V 


Consider the load distribution among various units of the plant 
as 

Pxi = a? Px 
Qxi = aj Qx 


af and a? are distribution factors such that £a? = 1 

X X 1 

Za? = 1. Current supplied by the ith unit will be 

i = 

i ' ’ 

and the internal bus voltage, from Fig. 4.8, will be 


and 


(4.39) 


EiZ^i = v^/e + jx'. i. 


(4.40) 


Power developed at the internal bus of the generator can be 
expressed as 


P. + iQ. = E./6. . i. 

X X — 1 1 


(4.41) 


4.3.2 Numerical Example 

The system described in Section 4.2 is considered for the 
analysis. In order to investigate the effect of interaction 
between machines in the plant, an identical machine is 
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considered to operate in parallel with the existing machine 
of the plant. PSS of Table 4,3, designed at = C.6 p.u,, 
is employed on both the machines with appropriate signals. 
Stable P— Q regions are then obtained for different system 
transmission strengths, i.e., for values of of C.2, C.4, 

C.6 and 0.8 p.u. In obtaining the stable P~Q regions it is 
assumed that both the machines are equally loaded. Figures 
4.10 to 4,14 represent the stable P-Q regions obtained in the 
two machine case. P^ and Qj in these figures represent the 
per unit plant output. 

Sp eed Signa l 

Figures 4.10 and 4.11 compare the stable P-Q regions 
obtained with speed signal in case of single machine and two 
machine systems. In Fig. 4,10 stable P~Q regions for tv/o 
machine system are obtained using individual speed signals 
while in Fig. 4.11 the regions are obtained using the average 
speed signal. Figure 4.12 compares the effectiveness of 
individual and average speed signals for the two machine case. 

Power Signal 

Picture 4,13 represents the stable P— Q regions obtained 

with power signal in case of single and two machine systems. 

In the. two machine case, individual and. average power signals 

giv© same stablo P-Q regions. Curves with bold, linos roprosont 
tho notv;ork limits. ... 





P-j- in p.u. 


FIG. 4.1 2 STABLE P-Q REGION Will 
AVERAGE SPEED SIGNAL 
SYSTEM. 
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Frequency Signal 

Figure 4.14 represents the stable P-Q regions obtained 
with frequency signal in case of single and two machine 
systems, 

4.3.3 Discussion of Results 
Speed S ignal 

Figure 4.10 shows that stable P-Q region obtained in the 
case of the single machine system is larger than that obtained 
in the two machine system, for all values of , when the 
signal fed to PSS is the individual machine speed. The 
difference in the stable P-Q regions is significant for 
stronger system and reduces with the reduction in system 
strength. This shows that the presence of other machines in 
the generating plant reduces the effectiveness of PS3 utili- 
zing individual speed signals. If, instead of individual 
speed signal, average speed signal is used in the two machine 
system, the stable P-Q region obtained is larger than that 
obtained in single machine system as shown in Fig. 4.11. Again, 
the difference in stable P-Q region is significant for stron- 
ger systems. Thus, in case of plants having more than one 
machine, average speed signal increases the effectiveness of 


PSS. 
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? ov/er Signal 

As in the case of single machine system with power 
signal, the network limit is reached before the stability 
limit for individual and average power signals, except for 
X|. = C.2 p.u. Curves with full lines in Fig. 4.13 represent 
network limit. For = 0.2 p.u. the stable P-Q region in 
the two machine system is the same for individual and average 
pov;er signals. The stable P-Q region for Xj = 0.2 p.u. in 
the two machine system is less than that obtained in the 
single machine system as shown in Fig. 4.13. This shows 
that effectiveness of PJ3 utilizing power signal reduces 
with the presence of oxher machines in the plant. 

Freq uen c y S ignal 

Figure 4.14 shows that stable F~Q regions obtained with 
frequency signal for the single machine and two machine 
systems coincide for all values of Xj except for 'Jij = 0.2 p.u. 
For Xy = 0.2 p.u. the stable ?-Q region obtained in two 
machine system is more than that obtained in single machine 
system. Thus, presence of other machines in the generating 
plant improves the performance of PSS utilizing frequency 
signal. 

An inspection of closed-loop eigenvalues (which are not 
reproduced here) at the stability boundary reveals that 
instability in every case is due to the shifting of exciter 
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node tc the right half of the coniplex plane. Eigenvalues 
corresponding to the rotor oscillation node remain in the 
left half of the complex plane. 

In the two machine system, with speed and power signals, 
where individual or average signal can be used, average speed 
signal is nore effective than individual speed signal as shown 
in Fig. 4.12, But individual pov/er and average power signals 
give the same P-Q region and are thus equally effective. 

In two machine case with plant signals, that is, average 
speed, average power and bus frequency signal, it is found 
that (with identical machines and with equal loading) six 
of the closed-loop eigenvalues are unaffected even though PSS 

m 

is provided on both the machines. This suggests that the 
system is not fully controllable. This situation is analysed 
as given in Appendix F and proof is obtained in the support 
of these results. 

4.4 COiiCLUSIONS 

A case study on the effectiveness of the various feedback 
control signals that can be used in the design of PSS for the 
iuprovement of dynamic stability has been presented. Rotor 
speed, electrical power output and bus frequency are the three 
signals considered for the study. The effectiveness of signals 
is judged from the stable regions obtained using the above 
signals as input to PSS. A comparison of the- effectiveness of 
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Pia desigr.Gd using classical control theory and pole assign- 
nent technique has been made. The effect of the prosarxe of 
other TiachinGs in a generating plant on the syster. perfor- 
nanco with different P3o is studied. 

From the results it can be concluded that ?3o utilizing 
speed and frequency signals should be designed for weak sys- 
tems. In case of power signal the results show that it is 
immaterial what value of ac system strength the Poo is 
designed for, in terms of stable P-Q domains. 

Povf'er signal is most effective for all operating condi- 
tions in comparison Vsfith speed and frequency signals. Effe- 
ctiveness is more significant for strong systems. For strong 
systems speed signal is most inefficient. For weak systems 
and higher loading conditions speed signal is bettor than 
frequency signal, however the difference in stable P-Q 
regions is small. 

In the example considered it is observed that with 
speed and frequency signals, effectiveness of Fio increases 
with the presence of other machines in the generating plant. 

But with power signal effectiveness decreases with the 
presence of other machines. For the two machine case 
average speed signal is better than individual speed signals 
but individual and average power signals arc equally efficient. 



CH.APT5R 5 


/,iODHLLI.:3 OF . :ULTL«ACKIN5 SYST3A AMD OiOICE OF P3S 

LOCATXO!’ 

5 . 1 INTRODUCTION 

In the previous chapter, development of system nodal is 
presented for a generating plant having n number of machines, 
connected to a large power system. The formulation is general 
enough to be applied to any multimachine system vjith the condi- 
tion that (a) the network is represented by reduced admittance 
matrix which retains only the internal buses of the machines 
and (b) the transient saliency of the machines is neglected. 

The formulation assumes that the loads are represented by con- 
stant impedances. It is not always possible to reduce the 
network to the generator internal buses, oarticularly when the 

identity of the- loads is to be retained, i.’.oreover, v/hen the 

# 

saliency of machines is considered, it is not oossible to 
represent the deviations in active and reactive powers as 
functions of state-variables only, as given in eqns. (4.13) 
and (4.1'^). Each machine, alongwith its excitation system in 
the formulation, is represented by a fourth order model. No 
freedom is available to increase or decrease the order of 
machine model. /»part from the raethod proposed in Chapter 4, a 
number of methods are available for the formulation of system 
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matrix in a multimachine power system [10-15]. In all the 
methods, the power system network is reduced to eliminate the 
non-generator buses by representing loads as constant impe- 
dances. 

In this chapter a general formulation, similar to that of 
Ref. [l6], is presented for the development of system model in 
a multimachine system. In the formulation presented here, the 
identity of the entire power system network is retained by 
utilizing the network Jacobian which is calculated during the 
load-flow analysis by Newton’s method. This approach helps in 
the representation of nonlinear voltage dependent loads vAich 
otherwise are represented as constant impedances. Freedom is 
available for the representation of each machine either by a 
detailed model or by simplified models as per the requirement. 
The system model of the multimachine system is obtained sys- 
tematically by developing the models of various machines and 
their interconnection through the netv/ork Jacobian. 

The second aspect which is discussed in this chapter is 
the choice for proper locations of PSS in a large multimachine 
system. A proper choice of PSS location is essential to get 
the maximum benefit from the point of view of improving dynamic 
stability in an interconnected power system. In a large scale 
power system the selection of effective locations of PSS is a 
complex problem and is normally decided by the economical and 
technical considerations. Incorporation of PSS on all the 
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'•.achines of a .Tiultimachine power system is not essential as it 
v/ojld bs uneconomical and also it may not be desirable from the 
point of view of coordination between various machines in a 
large system [17], For satisfactory system damping it is 
usually sufficient to equip relatively few machines with ?33. 
The optimal location of PSS has been investigated in various 
references [17”2C3 . de hiello et al. [17] suggested a method, 
based on eigenvalue-eigenvector analyses, for the selection 
of the machine which is the most effective candidate for 
stabilization. No conclusive results were given about the 
effectiveness of the method. Reference [20 ] oresents a sequen- 
tial procedure for the selection of machines to be equipped 
with PSS. The selection is based on the sensitivity of cri- 
tical eigenvalues to the increase in the coefficient of a 
danping term v/hich is inserted in each equation of motion, 
in succession. 

Two analytical techniques are presented hers for the 
selection of PSS locations in multimachine power systems. 

The first technique is independent of the feedback controller 
structure and the second one is based on eigenvalue sensiti- 
vity analysis [18]. A case study of a 13-aachine system is 
presented fox the selection of PSS. Root-loci of the closed- 
loop eigenvalues are obtained to validate the results obtained 
through the proposed techniques. 
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5.2 3Z73L0P'/,£:JT OF SYSiai MCDHL 


Tha fornulation of system matrices A,B and C in a nulti- 
■T.achinc- aowor systera is based on identifying the interconne- 
ct icnc bef-veen various machines through the pomi system net- 
work. each generating unit comprises of synchronous machine, 
excitation system, turbine and governor systems. For dynamic 
staoility analysis, dynamics of turbine and governor systems 
are neglected because of the large time constants associated 
ydth them. Network is represented by its Jacobian matrix, J 
as 


where 


and 


hS = 



AS = 

[ 

AS| ... 


^4 = 

[ 'SZJ 

A^ ... 


AS, = 
1 

[ 

AQil^ 


AZ- = 
1 

[ A 8. 

X 




( 5 . 1 ) 


The matrix [j] is of (2N, 2N) dimension, vdiera N is the number 
of buses in the system. The effect of any nonlinear voltage 
dependent load can be considered by modifying the coxresponoing 
diagonal block of J. 

In order to utilize the network equations given in the 
form of eqn. (5.1), the machine model has to be compatible. 

The linearized machine model is derived in a form such that 
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A? anc Aq are the inputs from the network, end the machine 

y 

terminal quantities A 8 and are expressed in tc-rciS of thc- 

machine state variables and Ap , Aq . a general synchrcncus 

9 9 

machine model [ 73 ] is utilized here with some modification 
where the Park's components of stator currents Ai. and Ai arc 
input variables. Thus, in order to make the machine model 
compatible vdth the network modal, the stator currents Ai^ and 
A i^ are expressed in terms of Ap^ and ^Qg* This representation 
eliminates the use of transformation of the variables from 
machine reference frame to the coaraon reference frame because 
APg and A remain invariant with change of reference frame. 
The only assumption made here is that the network is in quasi- 
steady state, which is well accepted in stability studies. 

TI:le development of the system model is based on the formu- 
lation of individual machine models and combining them through 
the network equations given in (5.1). 

5.2.1 Generating Units 

Synchr on ous Machine 

The general synchronous machine model [78 ] with some 
modification (see Appendix D) is used to represent the syn- 
chronous machine. The linearized dynamical equations of the 
machine, (eqn. (D.25)), are given by 



p X, = [a ] ^ + [3,„3 + b_„ ^ (5.2) 

^ — m ^ ti?-* —Cl '- ap** —g ~-ine id —30 m v * 

A A 

where x„ is the vector of state variables. Crder of x_ Depends 
upon the details used in the representation of synchronous 
machine deoending upon the availability of data anc the need 
for suitable accuracy, usually used in dynauic stability 
analysis are 

Detailed representation : 

Sinplified representation: 

Classical representation • 3^ 

The generator terminal bus voltage phase angle and magnitude 
are expressed, (eqn. (0.24)), in the form 

AZ„ = [D^j + Eb^3 ^3^ (5.3) 

-~g s ”*01 p —g 

w^ere 

= [AO AV.]'^ 

— g ^ 

The derivation of eqns. (5.2) and (5.3) is given in Appendix D. 
Excitation System 

The representation of various excitation systems used, 
is compiled in IEEE Committee Report [9 ]. The state-space 
model of the excitation system is represented, (eqn. (0,32))* 
in the form 

px = fAlx + b AVr 4-b AV x'"*'b .u 

^ ‘"6 ^ e"^ “-ev t — eu ref —eu 


(5.4) 



1^3 


and 




?/r,ar2 x^, u and are xes’Dectiveiy the state, input ano cut“ 
put variables of the excitation system* can be expressed 

in terius of state variables of machine and A3 as 

y 




ab-* -in 


(5.6) 


v^here 71^^ and are submatrices obtained respectively fro® 

A , 

D and D of eon. (5.3). 

p s 

Co mbined Model o f Generat ing Uni t 

Combining eqns. (5.2), (5.4), (5.5) and (5.6), the 
state~space model of the generating unit is obtained as 


p 4! = -‘2 * '’'ref + ^ (5-7) 

where 



‘■Vfl 

me e 



” B 1 

ap 

^ ~ 

b . '1 , 

L_ ev ao 

A 

6 -J 

» 


fev*''^aj 


b = Cb"^ 0]^ i b = [0 b^,j^ 

and the terminal condition of the generating unit is given as 

= [D^] + [D_] AS^ 

— g s — m p “"'g 

whare 


( 5 . 8 ) 





The values of the various variables, apoearir.T in tne 
eleaents of and b at the operating point are 

calculated from load-flow results as discussed in npoaneix li. 

5.2.2 Network Representation 

The network is represented by its Jacobian Jiatrix in the 
polar form as given in eqn . (5.1) . For an M-bvs power systei* 
network, the Jacobian matrix is of (2N x 2M) dimension. The 
procedure for evaluating the elements of Jacobian matrix is 
given in Apoenclix G, '/ithout loss of generality, it is 
assumed that the elements of the Jacobian are ordered and 
arranged in such a form that the first 'n* nodes are genera- 
tor buses and the last ’m’ nodes are load buses. The network 
eqn. (5.1) can be rearranged as 




iJa Cj-L 



ASl 


'^LG ^Ll 




where Jq-. Jql* Jlq and Jj^j^ are matrices of order (2n x 2n), 
(2n x 2in) , (2m x 2n) and i2m x 2m) respectively. 

5.2.3 Load Representation 

The usual constant power, constant currant or constant 
impedance type loads or any other voltage dependent nonlinear 
load can bo represented in the general fortii 
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= Kp 7jP (5.10) 

Si = Vj'! (5.11) 

v/here the coefficients I<p and and the exponents and n„ 
denend upon the type of load under consideration. Linearizing 
sqns. (5.10) and (5.11), we get 


*f.£ 


(5.12) 


where 


*2. = [*®. 


and 


[A 1 = 

£ 


n 


v"p"^ 




0 n K 

q q *■ 


Hqn. (5.12) is used to eliminate ASji corresponding to the load 
buses in the network eqn. (5.9). ASf in eqn. (5.12) represents 
the change in load power at a load bus and as per convention in 
load-flow analysis, load power is considered as negative inje- 
ction. Thus, I eqn. (5.12) is one of the entries of ASj^ 
in eqn. (5.9) with a negative sign. 


5.2.4 State-Space Model of the Overall System 

For the analysis of large scale power systems consisting 
of large number of generating units and loads, models of the 
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inoividual components are assembled together to get the 
state-space model -of the entire system, liith 'n* number of 
generator buses and 'm* number of load buses, eqns. (5.7), 
(5.3) and (5.12) can be written respectively as 


P X.; 

* |s*i 

II 

[Bp] A_S^ + [; 





(5.13) 

AZ-, 

II 

m 

[Dp] ^3q 

(5.14) 

and 




^§.L 

= AZj^ 


(5.15) 


where 


X-j = [x^, 


t it 




ml '■m2 •** 




... 




C'^Vrefl "• 


u 


[uj. 


n** 


AT,, = C At , At o *.• at_ 1 

^ ml m2 mn-^ 


AZG = [AZgl^2 ••• 
^-S-L = [ 

~ f *^11 *-t2 ■■• 


l.atrices Aj^» Bp, B, and Bp ail have block-diagonal 

structure with *n* diagonal blocks as i 
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= 31oc4 diag. 

[3pj = ilock diag. Bp 2 * •••» 

etc., while is also a block-diagonal inatrix v/ith *. 7 .* diagonal 
blocks. !'iov/, oliiTiinating ^2^ from eqns. (5.9) and (5.15), we 
get 

(5.16) 

O kH' mI' \3 


v/here 


'GG- 




From eqns. (5.14) and (5.16), ^Zq and AS^ can be expressed in 
terms of systera state variables as 


A 7^ 

- O 


= c[ij - [Dp] [j.TjJ) 4i 


(5.17) 


and 




(5.18) 


v^here I is an identity matrix of (2n x 2n) dimension. Substi- 
tuting the value of As^ from eqn. (5.18) into eqn. (5.13), the 
state-space model of the overall system can be expressed as 


P 4i = [Ai 4, + [3j y + [B] + [Bj,;] At. (5.19) 

W h 0 iE*'© 

Ca] = + [Bp] ([l] ■” “^3*^ 




5.5 ZrLlZl OF ?53 LCCATIOi 

It is desired to provide adequate damping to the electro- 
raechanical oscillations of the rotors of synchronous machines 
'ey proper choice of ?33 location. For an n machine system, 
the^e arc in general (n-i) codes of rotor oscillation. For 
planning studies it is sufficient to use a simolified systew 
model containing only relevant features. Classical model 
of machines, that is a constant voltage source behind its 
transient reactance, retains the information regarding electro~ 
mechanical oscillations and is thus considered to he adequate 
for the development of systera model. '.»ith this assumption, 
the systera model given in eqn. (5.19) will take the form, 


v/here 


and 


p 'ft = [a] X + [b] y 



(5.20) 


^'ormally, rotor velocity of the machines is used as feacback 
control signal for P3S. Thus, the output equation of the entire 
system is 


Y = [Cj X 


(5.21) 



I5i 

■ f.'jie tl'.e output natrix C has a block diagonal forr.i .‘dth the 
"'iartrial block 

q = [1 0] 

iitY the systeni defined by eqns. (5.2C) anci (5.21), the 
orocle.j for the choice of P33 location can bo stator’ as : 
select the smallest set of locations for PSS, froit the n 
oos sibiii ties that exist, such that adequate damping is oro- 
vidsd, v/ith economical control effort, to the critical nodes 
of the system 

n 

p X = [a] X 1 b, U. (5.22) 

i=l ^ ^ 

by the incorporation of PSS at these locations. 

The two methods presented here are based on two different 
assuraptions given below ; 

1. The controller structure is not decided apriori, 

2. The controller structure is decided in advance* 

In the first method, to get the maximam damping to a 
critical mode X., the scalar products of V/. and b. are 

J J 3. 

obtained for all values of i , v^ere is the reciprocal 

eigenvector of . the system matrix A corresponcing to j th 

eigenvalue. The value of i for '.#iich magnitude of the scalar 

product, l<b’.,b. >j, is maximum decides the PSS location. 

J ^ 

Justification of this statement is given in Apoendix H. 



1:2 


In the second method, the controller structure is 
cecidad as : 


“i = h fi 


(5 .23) 


ax. 


and ths eigenvalue sensitivity, , ic obtained for all 

^ i 

values of i. The value of i for '.vhich the real part of 
fiX /dp has the maximum magnitude decides the PSS location. 

J A 

Justification of this statement is given in .=\po3ndix H. 
eigenvalue sensitivity is given as 


S . . 


dX./dF. 
3 1 


<b. 

—1 


V . , ii . > 

"1 .rol rx.. 




(5.24) 


where C- is the ith row of output matrix C, V. and /. are 

^ J J 

respectively the eigenvectors of matrices 4. and A corres- 
ponding to the jth eigenvalue. 

5.4 rlLIdllrllCAL uXAliPLE 


.1 case study of a large-scale power system consisting of 
13 machines, 71 buses and 94 lines [18] is presented here. 

It consists of 5 ther'sl units and 8 hydro units. A single- 
line diagram of the system alongwith the systera data and 
operating conditions are given in .Appendix I. The follcvAng 
assumptions are made in deriving the system model. 

1. The loads are represented by constant ir.ipecancas • 

2. For the sake of sisrpiicity, camper circuits and transient 
saliency of all the aiachines are ignored. The damping 

^ coefficient is considered to bo zero for all the 


machines 
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Ths cperating point is obtained fror^ load flov st^oies 
and f'-e Jacobian of eqn. (5.9) is constructed as given in 
Ape anti n 3. 

5.^.1 ii gen value .Analysis 

The systera matrix ’a’ is developed, using tne techniguo 
presented in Section 5.2, for the foilowinn three cases and 
eigenvalue analysis is carried out in each case- 

1. All the machines are represented by second order classi- 
cal model. 

2. .All the machines are represented by the third order 
laodel by including field flux decay. 

3. Machines 1,7 and 11 are reoresented by a third order 
modal along with a single time constant excitation 
system (I3E3 type 13 » --.ppendix D) . The gain and time 
constant of voltage regulator for each machine are 
respectively selected as 50 and 0.02 sec. Best of the 
machines are represented by the classical model. 

eigenvalues are computed for the three cases and are 
given in Table 5.1. 

The first twelve pairs .of comolex eigenvalues in Table 
5.1, for each case, represent the modes corresponding to 
electromechanical oscillations. In the first case, each 
machine is represented by the classical model and as damping 
has been neglected, pure oscillatory modes are obtained. 



Table 5.1 


Open-Loop Eigenvalues for the i2-.*lachir.2 
Oystem 


Oa s e X 


Case 2 


O.C 

0.0 

0.0 


p Pf 


O.C 

0.0 

0-0 

0.0 

0.0 

0.0 

0.0 

0.0 


T ji *. 

± J12 

± 3^^ 
-r jlO 
+ jlO 
T jlO 

± 

± jio 

.t 3 9 
±38 

t 3 7 
±36 


.05998 

.82911 

.83337 

-80104 

.57432 

.46123 

.34287 

.17931 

.40440 

.83647 

.37834 

.11920 


-0 . 0:^30 + 
-O.C3661 + 
-0.01742 
-€■ .00729 i 
-0.00962 j; 
-O' .00645 ;> 
“0.00512 + 
-0.01293 + 
-O .00231 .-r 
-U .01212 j, 
-0.01261 + 
—L' .0309 1 p 
-0.23862 
-C. 29501 
-0.30725 
-0.33450 
-C.378S1 
-0.40155 
-0.41425 
-0.42302 
-0.47232 
-0.53329 
-0.65403 
-0.78346 
—0 .9*+4i5 
O.OCXXX) 


jl^v .05321 
j 12. 32702 
jll. 33252 
j 10. 30095 
j 10. 57362 
j 10 .46102 
j 10. 34252 
jlO. 17328 
j 9.404-21 
j 8.33592 
j 7.37752 


j 6.11762 


0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0. 
0 i 
0 . 
0 
0 , 
-5 
—6 
-6 
-4-3 

"*'*“f*T 

0 


Case 3 

02892 + j 14 .05537 
C4537 + jl2.S3124 
C0S52 + jl 1.33472 
00988 .t j 10. 3027 3 
CC045 4 - j 10. 57437 
0067^ ^ jlO,-^‘6237 
01693 i j 10. 34604 
01-563 + j 10. 18523 
0OC40 + j 9-40696 
01546 + j 8.34131 
03X-2 j 7.39500 
05525 + j 6.15929 
78423 
0933C 
33359 
87653 
.58668 
.99036 
.00000 



3igsr.\’ai'je3 sccained in case 2 indicate that tha system is 
cync.'.'iic sliy stable and tliat the field-windings contribJte a 
s:;.all ar.cjnt of da.nping. The addition of fast acting excita- 
tion ays tens on the selected machines introduces nogcxive 
'■J'li'.oiny to all the nodes of rotor oscillcticn as seen fro.n the 
eigar. values obtained in case 3. 

5.4.2 Choice of 233 Location 

The system model described in eqn. (5.20) is formulated 
according to the procedure given in Section 5.2, representing 
each machine by the classical model, dotor soeed is consi- 
dered as the feedback signal for the calculation of eigenvalue 
sensitivity in eqn. (5.24). All the machines, taken one at a 
time, are considered for the location of P53 vath the objective 
to damp the electromechanical oscillations. In the first 
method, magnitude of the scalar product of the reciprocal* 
eigenvector Uj and the vector is calculated for each pair 
of the complex eigenvalues, and for different locations of PSS, 
that is, for all values of i. The scalar products for different 
locations of PSS are given in Table 5.2. In the second method, 
eigenvalue sensitivity is calculated from eqn. (3.24) for each 
pair of the complex eigenvalues and for different locations of 
PSS. The eigenvalue sensitivities are given in Table 5.3. 

Inspection of Tables 5.2 and 5.3 shows that the tvjo 
;nethods, suggested for the determination of effective PSS 


T3..1c 5.2 ..isgnituda of Scalar Product <ijj b. > for 
jiffsrent Locations of PSS 


?33 
L oca 


tion 

(i) 

w=14 .05998 

12.82911 

1 

0.4C331 

5.82663 

2 

1 .81758 

5.58953 

3 

0.12794 

0.21663 

4 

0.31270 

0.43132 

5 

0.55664 

0.69737 

6 

31.80959 

32.54492 

7 

26.94005 

27.17332 

8 

0.06791 

0.59849 

9 

0.11791 

1,52749 

10 

1.58906 

0.69276 

11 

0.33487 

2.05883 

12 

0 • 1007 1 

0.61253 

13 

G.C7338 

0.38369 



for j=l 

,6 


11.83387 

10.8Ca04 

10.57432 

10. 

8.71262 

12.94110 

5.91329 

10 

20.30348 

0.53353 

1.62319 

2 

1.S3405 

1.13433 

2.14303 

15 

2.37522 

0.63677 


6 

2.31093 

27.71291 

1.13545 

1 

0.50501 

1 .25031 

1.24167 

0 

1.46061 

0.33290 

2.31675 

0 

0 .02936 

0 .05923 

0.21050 

0 

0.08793 

0.12517 

0.54344 

0 

0.19437 

1.01525 

11.79566 

6 

0.15234 

0.13377 

5.36513 

3 

0.04435 

0.03587 

0.62303 

0 

0.05373 

0.27381 

20 .28680 

3 


46123 

7709S 

78C22 

23010 

81476 

.15367 

76357 

29446 

.12716 

.376C© 

.64598 

.020C-3 

.11939 

.23166 
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Table 3.2 ^ continued ...) 


?3S i<lV. , b. >, for j = 7,12 

Loca- -j "* 



'-)=iC. 3^237 10.17931 9.40440 3.33647 7.37334 6.11920 


1 

1 ,CCC39 

10.08257 

14.54410 

12 .40288 

C.C92C5 

13.17896 

2 

1 .C9C20 

6.68621 

6.54314 

1.99559 

0 .41993 

14.25138 

3 

1.17253 

3.38395 

0.30858 

1 .74902 

C -19675 

3.CC733 

4 

0.39402 

4.52835 

8.57733 

16.31328 

0.98267 

10.64746 

5 

0.34826 

2.79837 

3.57630 

3.74053 

0 « x. X X so 

4.41169 

6 

0.65268 

3.23808 

2.10063 

5 .14494 

0.71179 

1.81535 

7 

0.91773 

5.03313 

0.58178 

5.50522 

0.53124 

1.55011 

8 

0.06055 

0.69631 

11.08748 

C. 85931 

0.53172 

0.97262 

9 

0.10742 

1.1767S 

6.11219 

0.93190 

12.90535 

8.83553 

10 

4 .08823 

25.29693 

3.72950 

8.52706 

3.27959 

7.38256 

11 

13.53257 

5.09114 

4.76912 

3.66365 

12.60648 

13.21443 

12 

12.99066 

2.99232 

1.55000 

1.21674 

3.93132 

4.12291 

13 

0 .43437 

2.85502 

0.30512 

0.28569 

2 .27224 

2.99928 
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Table 


Eigenvalue Sensitivity for Different 
Locations of PSS 







Table 5.3 (continued) 


PSS 

Loca-_ 


% 

for j = 

tion 

(1) “ 

=10.34287 

10.17931 

9.40440 

1 

-0.00429 

-1 .07824 

-0.42762 

2 

-0.00797 

-0.54539 

-0.09461 

3 

-0 ,07449 

— 1 mQ6&)7 

-0.00131 

4 

-0 .01442 

-0.64578 

-0.36242 

5 

-0.00346 

-0.35382 

-0.08417 

6 

-0 .00785 

-0.36217 

-0.01051 

7 

-0.01520 

-0.83706 

-0 .00377 

8 

-0.00180 

-0.30414 

-10.96793 

9 

-0.00071 

-0.10177 

-0.37832 

10 

-0.23157 

-15.36443 

-0.11936 

11 

—4 • 10^32 

-0.77198 

-0.10811 

12 

-12.23194 

-1 .01880 

-0 .04056 
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7 to 12 


8.83647 

7.87834 

6.11920 

-3.41157 

*C3CXX)8 

-1 .70626 

-0.16075 

-0 .00332 

-2 .32014 

-0.34524 

-0 .CX)481 

-0.71375 

-10.30540 

-0.03856 

-2.98398 

-0.78071 

-O.CXX)58 

-0.59447 

-1 .02919 

-0.01523 

-0.03630 

-1 .02836 

-0.00876 

-0.00031 

-0.37378 

-0.10379 

-0.2^13 

-0.13488 

-9.82789 

-3.98238 

-1.73276 

-0.17001 

-1.12578 

-0.70465 

-4.85514 

-5.50703 

-0.25145 

-1.54218 

-1.74274 

-0.02246 

-0.57739 

-1 .07597 



1 a zl 3 


Ranking of 


5.4 


Mode of 
Oscillation 



Ranking ol 
Descendinc 

S.No. w 

1 

2 

3 

4 

5 

1 

14 .05998 

6 

7 

10 

2 

1 

2 

12.82911 

6 

7 

1 

2 

5 

3 

11.83387 

2 

1 

5 

3 

4 

4 

10.80104 

5 

1 

3 

6 

4 

5 

10.57432 

13 

10 

11 

3 

1 

6 

10.46123 

3 

1 

4 

10 

13 

7 

10-34287 

12 

11 

10 

3 

13 

8 

10.17931 

10 

1 

3 

12 

13 

9 

9 .40440 

8 

1 

9 

4 

10 

10 

3.83647 

4 

1 

10 

6 

7 

11 

7.87334 

9 

11 

12 

13 

10 


12 
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Ferent Locations of ?3S 


F Different Location of ?33 in 
3 Order of Effectiveness 


6 7 8 9 10 11 12 13 


5 9 3 4 13 12 11 8 

3 4 12 13 9 11 8 10 


7 10 6 11 13 


9 12 8 


2 11 12 13 

7 5 12 6 

11 2 5 6 

7 4 2 6 

7 11 4 2 

10 2 5 12 

5 11 8 3 

8 4 6 7 


9 7 8 10 

2 3. 4 9 

9 8 7 12 

15 8 9 

6 5 8 9 

6 7 3 13 

12 2 9 13 

3 2 5 1 


1 10 13 3 


5 8 


6 7 
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locations, give almost similar results. Ho*vever, the eigar- 
val'ja sensitivity method gives acditional inf ormaticn 
regarding the sign of damping orovided to any mode. 

Table 5.4 sho'.vs the ranking, in the descending ordar of 
ef Tc-ctiveness , of the PSS location for da.aping of each mode 
of rotor oscillation. Table 5.4 is obtained from the eigen- 
value sensitivities given in Table 5.3. 

5.4.3 Discussion of Results 

The two proposed methods for the selection of effective 
PSS locations give closely agreeing results as shown in 
Tables 5.2 and 5.3. Fro© the results shown in Table 5.3 it 
is found that for some of the PSS locations eigenvalue sen- 
sitivities are positive. This shows that providing PSS on 
such machines will result in negative damping of certain 
modes, while improving damping of other modes. This is 
expected in large power systeins [17 ]. However, in the 
exaaiple considered, their numerical values are negligibly 
small in comparison with the sensitivities of other inodes 
for the same PSS location. 

An inspection of the elements of the first column of 
Table 5.4, that is, PSS locations with ranking one, shows 
that no FS3 location occurs twice except location 6. Loca- 
tions 1 and 7 are not present in this list. This shows that 
each node of oscillation can be controlled by equipping PSS 
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or- 3 liffsient machine. However, the list of r3S locations 
•vith ranking t”'o, that is, the second column of Tatle 5.4, 
shows that location 1 occurs six times, locations 7 and 11 
occur twice. This suggests that equipping machine 1 with ?G3 
>/ili control six modes. In other v/ords, the most effoctivo 

location of PS3 is machine 1 . Next in the order of effective- 
ness are machines 6,7 and 11 . Looking into the nacerical 
values of eigenvalue sensitivities for P33 locations 6 and 7, 
it is found that both locations provide almost the same damp- 
ing to first and second modes, suggesting that either location 
can be chosen for equipping F33. Thus, by providing PS3 on 
machines 1 , 6 (or 7) and 11 all the modes except fifth can be 
sufficiently damped. The fifth mode can be damped fay providing 
PSS on machine 13 . These results are based on the eigenvalue 
sensitivities obtained at zero feedback gain of controller as 
Vj and rfh in eqn. (5.24) are eigenvectors of the open-loop 
system matrix. The sensitivities change with feedback gain. 
Consequently, the open-loop sensitivity data can be used for a 
preliminary choice of PSS locations and analysis of the system 
after incorporation of PSS at these locations will indicate 
whether any modification in the choice is to be done. 

To verify the results of Table 5.4, closed-loo? eigenvalues 
of the system, described by eqn. (5.20), are obtained for varia- 
tion in the feedback gain using the control law as 
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figures 5.1 to 5.4 show the loci of cicsec-lcop sicenvaiuos 
for variation in feedback gain considering one nachine on 
control at a tirce. 

Inspection of Figs. 5,1 to 5.4 confir/as the results 
obtained f ron Table 5.4, Loci of eigenvalues show tht^t each 
PS3 location provides large damping to only one of the nodes, 
in accordance with the list of P3S locations with ranking one 
given in the first colunm of Table 5.4. however, locations 
1,2 ,4, 6, 7, 9 and 11 provide damping to more than one mode. 
Locations 6 and 7, both give similar results for first and 
second modes, but location 7 provides more damping to first 
mode in comparison with location 6 (see Fig. 5.2). An 
inspection of Figs. 5.1 to 5.4 shows that by providing PSS 
on machines 1,7 and 11 all the modes, except fifth am’ ninth, 
can be sufficiently damped. Although, location 1 has rank 
two for the ninth mode in Table 5.4, the mode is unaffected 
by providing control on machine 1 as seen from Fig. 5.1, 

This can be explained from the relative values of eigenvalue 
sensitivities of this mode for locations 1 and 8. 

From the foregoing discussions and on the basis of 
results of Tabic 5.4 and Figs. 5,1 to 5.4, follovdng recommen- 
dations are made, as shown in Table 5,5, for providing PSS in 
the example considered. 
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Machin® S on control 


Modtin® 7 on control 


LOCI OF EI6ENVAUUES FOR VARIATION IN CONTROLLER 
GAIN (MACHINES ON CONTROL 4,5,6 AND 7} 
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FIG 5.3 LCX:i Of" EIGENVAIUES ¥Cm VARIATION IN CONTROLLER 
GAIN (MACHINES ON CONTROL 8,9AN010). 
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T 2 ''le 5.5 effective Locations of PS3 in the 13-..achine 
Systen 


I’o, of Machines to be 
Zcuipoed 7dth PoS 


?SS Locations 


3 

4 

5 


1,7 and 11 
1,7,S and 11 
1,7,8,11 and 13 


5.5 CONCLUSIONS 

A state space inodel has been develooed in this chapter for 
large niultimachine power systems. The model is based on the 
linearized machine and network equations around an operating 
point. Freedom is available to acconmodate any degree of 
detail for machine representation. The model can be used for 
eigenvalue analysis and design of PSS to improve dynamic 
stability of systems. 

It is found that proper location of PSS is itiiportant for 
the effective damping of rotor oscillations in multiraachine 
systems. Two methods are suggested for this purpose. The 
first method is independent of the feedback control structure. 
The second methexJ is based on eigenvalue sensitivity and takes 
the feedback control structure into account. It is shovr* that 
both methods give closely agreeing results. 
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A case study of a 13-machine system is presents; for 
aigenvelue analysis and for the selection of effective P33 
locations, effective ?S3 locations are obtained from the 
sucoested methods and validation of results is obtained oy the 
study of loci of eigenvalues. It is found that it is sufficient 
to ecuip relatively few machines with P35 for acequate improve- 
ment in dynamic stability. At the most five locations are 
sufficient to danp all the twelve modes of rotor oscillations 
in the example considered. 



CH.i^TER 6 


.^PLICATION CF DECBITRALIZHD OOU’TKlL ¥CA T/.‘ jcSlZl Or 
?33 Vii f.iULTL'.i;.CHIA3 ?C^ZR S'/SIBIS 

6.1 ii:tp.uix'ction 

Coordinated application and design of power syster. 
stabilizers in multimachine po«#er systems have dra’ivn much 
attention recently [17,20]. In multimachine systems the 
problem of dynamic stability can ’oe divided into two parts : 
the selection of generating units to be equipped with PSS 
and the design/tuning of stabilizers for the selected units 
such that the entire system is dynamically stable. The 
selection of PSS locations in a taultioachine system has been 
discussed in Chapter 5. For the selectee machines PSS can be 
designed by any one of the methods discussed in Chapter 2 , 
Direct application of these techniques is not possible for 
liiultimachine system, because of the practical problems of 
infoiroation transfer between various machines of the system. 
These techniques are to be modified to consider the decouoled 
nature of control for various machines. This leads to the use 
of decentralized control theory. 

In decentralized control, the entire system is divided 
into various subsystems having local control stations. At 
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each ccrtrol station, the controller utilir.-js only local 
systsn outouts end controls only local inputs. Ho "ever, 
all tha controllers are involved in controlling the sane 
large syste.r. Lefebvre et al. [77] have suggested a technique 
for the design of such controllers for eiganvaluo cssignr.er.t. 
In the first stage of this technique tha controllers are de-- 
signed for stabilizing the individual subsys teas neglac ting the 
coupling between theza* In the second stage, the interconne- 
ctions are introduced in saall incressents and the controller 
parameters are updated to retain the eigenvalues at their 
preassigned locations. Lefebvre [43 J has recently applied 
this technique for the P33 design in a three machine system 
with P3S designed on all the machines. 

Sequential algorithms have also been suggested for the 
design of PSS [20,47J using pole-assignment techniques* Un- 
fortunately the seauential -.addition of stabilizers disturbs 
previously assigned eigenvalues. 

* 

In this chapter, an investigation of PSS design in a 
large multimachine system is carried out. As described in 
Chapter 2, PSS design based on (i) eigenvalue assignment and 
{ii) the minimization of a performance index are considered. 

In the first case, coordinated design of PSS is carried out 
using the technique given in [77 j. In the second case, the 
method given in Chapter 2 is modified to consider the 
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decir.traiizec structure of ?S3. 3oth th^se rcsthots are cc::xc:.red 
using tuc- exaraoles of a 3~machine and a IS-nachina syste:*. rcr 
the 13’'-nachine systsit, the ?33 is designee for a fev; sslectac 
riachincs deterruined from the anjlycis prosontad in Chaotsr 5. 

6,2 PS3 DHSIGi B^SED ON cIGE'-r/.^a-UE ^3IGR,*£NT 

The objective of eigenvalue assignment by output f^ecloack 
is used in this section for ?S3 design. The design problen can 
be stated as : Design stabilizers for the selected niachines in 
a nultima chine system with dacentralized output feedback to 
assign closed-loop eigenvalues to preassigned locations. 

6.2.1 Formulation of Control Problem 

The state-space model of a multimachine system with n 
generating units, developed, in Sec. 5.2, can be represented as 


p 4i “ % 


vjhere 


r t t 

••• 

and 

Hi 

"" Ha2 

• * • 3 

m*' 


( 6 . 1 ) 


''ai 


represents the state variables of ith generating unit and 


its excitation system, .^s Is a full matrix and is a block 
diagonal matrix as explained in Chapter S. The dynamical 
equation of the ith generating unit has the form 
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^ -T.ii 


i « j X j 4“ 


(e.xi 


v;hore and are subraatrices of corrcspondin ; to ir.o 

ith unit and is the ith diagonal block of 3.^ . ^inilarly 
the measurable output (viz., rotor speed, electrical power or 
bus frsquerxy) of the ith unit can be expressed in terms of the 
state variables of different units as 


"^mi -mi 

A 


(6.3) 


The summation terms in eqns . (6.2) and (6.3) represent the 
interconnection of ith generating unit with other units of the 
system. The interconnections are parameterized, with 
0 < r < 1 as follows : 


P ~ ^ ^ ^i ^i C6.4) 

*J 

ri 

and 

ymi ^ -^-mi ^ 

^i 

With the introduction of parasi^ter r, the degree of interconne- 
ction in eqns. (6.4) and (6.5) can be varied. Clearly at r - 0 
the system is completely decoupled and at r = 1 it is fully 
connected. 



I?-' 


Irit^aily stabilizers for sel:;Ctcd nachi-.oc are 
no^lectir.g tne intorconnection , that is, v/ith r = C, /.-ith out- 
put foacback, using local output variables to aseior. tr.e 
closGC-loop eigenvalues using auy one of the techniques 
presented in Chaotcr 2. The interconnection is introduced in 
steps to uDdate the parameters of PSS in such a way that the 
closed-loop eigenvalues remain unchanged. Let the c’ynardcal 
equation and control law for ?33 of the ith generating unit 
be 


and 


P .?i = [Sil + Si Y„.i 
'^nu = t-iJ H * '’.rd 


( 6 • 6 ) 


(6.7) 


where represents the state variables of PSS. Ccwjibining 
cqns. (6.4) to (6.7), the dynamical equation and control law 
for the augronted system (including PSS dynamics) of the ith 
unit can be expressed as 

P »i = [ iiiJ *i + r + [3^] (6.8) 




J=X 

M 


and 

Vi = 

n 

_x. 4- r E [Cjdj ^4 

* 1 *1 )tJ 

(S.9) 


le 

11 

li 

(6.10) 
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where 


X. = [x“. 

-i ffii 1-* 


"ii 


A . , 0 

iTiii 

[0 0 




B, 


0 • * 
11 


'c„ . . 

in li 


I 

0' 


Lo 


tit 


- Ki 


A 


ij 


A O 

niij 


K. 


di 




Lft. 


'®ij 


Lo 




If ith machine has no stabilizer, then = ^mi* “ ^ai 
matrices A^^^, ^ij» ®i* ^ii will be equal to matrices, 

^m.i.i» \if ^mi» ^mij respectively. will be a 

scalar with zero value* 

For the composite pov/er system, eqns. (6.8) to (6.10) will 
take the form 


p X 

= [Aq + 

r A J X 4- [Bp] j 

y 


(6. 

11) 

I 

= [ia + 

r C„3 X 

0-* — . 




(6. 

12) 

and 








U 

= [%(r)] Y 




(6. 

13) 

where 

= Block 

Diagonal 

[ *11. 

'^22» 

**•» ^nn^ 



®D 

ss Block 

Diagonal 

C Bi , 

®2 * 

3 




Block 

Diagonal 

f Si- 

^22» 

•••» Cj^nl 




= Block 

Diagonal 

t K^l- 

*^d2 

-- J<dni 
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and the structure cf matrices and C axe as fciicA's t 

Q O 


0 

''12 

^13 * • • 

^In 


c 

'-12 

"'13 

-1 
•• “lr.j 

^21 

C 

^23 ••• 

Ao., 


^21 

c 

^23 

# 

•* -2n 

^31 

^32 

0 • « • 

^3n 

s [f-o] = 

^31 

^-32 

c 

*• "an 

• 

• 

0 

"^02 

^n3 

* i 

0 i 

• • 

0 

» 

♦ 

L-nl 

^n2 

^n3 

* 

• 

» ♦ G 


The objective is to find a decentralized output feedback 
control law, as given in eqn. (6.13), such that the composite 
power system depicted by eqn. (6.11), has specified eigenvalues. 
Using the parameterization introduced in eqns, (6.4) and (6.5) 
and defining Kp(r) to be the decentralized feedback matrix, the 
eigenvalues of the composite power system are the eigenvalues 
of the closed-loop matrix 

F{r) = Aq + rA^ + K^Cr) {C^ + r C^) 

At r = O the matrix Kjj(o), which assigns the specified eigen- 
values, is easily computed because the system is decoupled. 

The idea is to slowly introduce the connections using discrete 
increments Ar and to modify K^Cr) in such a way that the eigen- 
values remain unchanged in presence of the interconnections. 

The advantage of this approach is that updating K^Cr) to 
obtain iC|j(r +Ar) is much easier than trying to directly 
compute Kpd), the final answer. K^Cr + Ar) is obtained frc» 
K|^(r) as follows.; 



The ith eigenvalue equation of the cicsed-loop systex 


for any value of r may be written as 


[ - F(r)3 ViCr) = C 




where ’-/^(r) is the right eigenvector corresponding to the 
eigenvalue X^(r). The system also has a left eigenvector 
with the properties 


^j(r) [ - F(r)] = C 


( 5 . 15 ) 


and ^ fl for j = i 

V.(r) = d 

^ ■; I U f or J F X 

* 

where W^Cr) is the conjugate transpose of the left eigenvector 

Wj(r). Differentiating eqn. (6.14) with respect to r and pre- 
multiplying by W^(r), we get 

H*(r) ^ [\{r)I - F(r)j V.(r} -f wj^(r) [X^Cr)! - 


F{r)3 Iy ° 


(6.16) 


The second tenn in eqn. (5.16) v;ill vanish as evident from 
eqn. (6.15). Rearranging eqn. (6.16) and suppressing indica- 
tion of r-dependence of X., V., and Kri, we get 

# * 

dF ^i Hr-"- - ^ Woi I'i 

(6.17) 




In order to make the eigenvalues insensitive to r, in 

eqn. (6.17) can be set equal to zero. Then eqn. (6,17) can be 
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dK^ 

solved for which is, in turn, utilized for updating 

as follows 


dK-»C r) 

KqCx -r tx) = KjjCr) ^ — -| Ar 


( 6 . 18 ) 


Solution of eqn. (6.17) is obtained by using the tensor 
product operation. Suppose M ~ ® natrix and Lit 

rxs, then the tensor product^ of by L (<.* L) ,is defined 
to be pr X qs matrix, isi iSL ~ Cm. .LI . The *Vec’ operation on 

X ^ 

the matrix li = [Mj^, U 2 * •••» » where is a column of M 

and defined as 

M] 

Vec(M) M. 


M. 


It can be shown [77] that the set of equations (6.17) for 
i = 1 to M, where N is the total number of states in eqn. 
(6.11), is equivalent to t 


[d(r)] Vec 


dKpCr) 

B'r 


P(r) 


where the matrix ^(r) 


$ (r) 


[(Cq + xC^) VJ * 8 2 I 3d 

[ (Cq -f rC^) ¥23*“ ® I2 Bp 


[(=0 + xc^) a Bo 


(6.19) 
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and fre vector fir) is 



1 

df * 

# 

-^1 




pit)- 

i'l? 

Wi 

# 

--'?2 

[Aq 


i'2 


• 

• 

d^v 

dr" 

«• 

-* W., 

t^O 


• 


A 

As Kjj is a block-diagonal raatr-.x, coluians of ^(r) correspond- 

.. to .e.o e„t.eo o. Vec C^, ate a.ete. ao. .0. 

can be written as 


[^(r)] 


dK^, 

Vec(g^) 

Vec(^) 


Vec(^~) 


= Hr) 


( 6 . 20 ) 


dK -i 

Solution for Vec(-g 7 -=)» i = 1 to n can be obtained by computing 
the pseudo-inverse of the non-square matrix ^(r) in eqn.(6.20). 

Solution of eqn.' (6.20) during each iteration provides an 

dICQ 

update of which is utilized in eqn. (6.18) to update the 
gain matrix K^. The conditions that guarantee the existence of 
a real solution of eqn. (6.20) are [ 77 ] i 
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(a) Each subsystem (machine and PSS) is controllable and 
cbservafale, 

(b) T!;e eigenvalues of each subsystem can be assigned by 
local dynamic output feedback. 

(c) At most one subsyster; has an augmented state space wnich 
is of odd order (machine and PSS), This condition 
guarantees that the feedback gains, that is, matrix 

Kjj. are real. 

(d) -tiatrices 3-j and (Cp + rC^) have full coliaan rank and 
full rov/ rank respectively. 

6.2.2 Design Algorithm 

The following algorithn can be used for PSS design in 

multimachine syste.ms using the procedure under consideration : 

Step 1 ; Select the machines to be equipped with PSS. 

Step 2 : Obtain state-space representation of each machine as 
given in eqns. (6.4) and (6.5). 

Step 3 : Augment state equations of the selected machines to 
include stabilizer dynamics as given in eqns. 

(6,8) and (6.9). 

Step 4 ; Design feedback matrix K^j_(0) for each of the sele- 
cted machines using local output feedback, to assign 
closed-loop eigenvalues at desired locations, by 
setting the parameter r » 0 in eqns, (6.8) and (6.9). 
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Step 5 ; 


Step 6 : 

itep 7 : 

Step 8 : 

Step 9 ! 
Step 10: 


Construct sqns. {6.11), {6.12) and (6.13) for the 
composite power system from eqns. (6.6), (6.9) and 
(6.10) of individual machines. 

Calculate right and left eigenvectors cf closed-loop 
matrix F(r). 

Calculate elements of 9(r) and p(r) and solve ec,n. 

dK,. 

( 6 • 20 ) for • 

Update KqCt) from eqn. (6.13) by selecting a suffi- 
ciently small value of £t. 

Increase the value of r to r+ Ar. 

Go to step 6 if r is less than or equal to unity. 


6.3 PSS DESIGN BASED CM .’ll. 4 IZATICi^ CF PERFOR! AIXE IMDEX 

The objective of 133 design here is to minimize a perfor- 
mance index of the form : 

^ n r* 

J = / ( S a. Apf) dt (6.21) 

o i=l ^ ^ 

where Ap^ represents the fluctuation in power output of the 
ith generator and is the v;eightage to be attached to^P^^. 

The technique for PSS design based on the above objective 
is given in Chapter 2 for a single machine system. This is 
extended here for a multimachine system. As stated in 
Chapter 2, this objective is equivalent to minimizing the 
fluctuations in the power outputs of various machines following 
a disturbance. 



6.3.1 "emulation of the Control Probler. 

Eqns. (C.ll) and (6,12) for 'Ce cenposite cuv.or syster. 
(including stabilizer dynamics), with intercor.n-"»c iicnt, can 
be expressed as 


p X = [A] X -»■ [Bf 


( 6 . 22 ) 


X = [C] X 


(6.23) 


where 


A = Af» + A 
D 0 


C = Cjj + 


A and C are full matrices and 3 q is a block-diagonal matriac. 
The decentralized output feedback control law can be defined 


•Ck3 y 


( 6 . 24 ) 


#here K is a block-diagonal feedback matrix and can be expre- 

■■ . . - - stabilizers. The design 

,ne the optimum value of K 
eqn. (6,21) is met. 


.'ssed as 


(6.25) 
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Incremental change in complex power output of Generators, 
in a muitisachine system, can be expressed in terms of the 
system state variables as given in eqn. {5.16) of Chapter 5. 


v/here 


and 


63., 

= 4 



( 6 . 26 ) 


= - ''gg 

Dp)-^ J, 

0 


'Sg 

= 

-.t 

^g2 ••• 

gn-* 


“^gi 

= [ ^P, 

1 





and Dg are 

defined 

in Chapter 3. 

6p can be 


obtained from the vector of changes in complex power as 

ip = [Tj] 4, (6.27) 

whore matrix T2 is obtained from the matrix T^. System state 

variables can be expressed in terms of the state variables 
of the composite power system as 

4 = [43 i ( 6 . 28 ) 

From eqns. (6.27) and (6.28), vve get 

= [P] ^ (6.29) 

where 

Substituting the value of AP from eqn. (6.29) into eqn. 
(6.25), the cost function, J can be expressed as 



13 - 


where 


J 

A 


/ 3 X) dt 

0 "" 


Q = [I J [aj [T] 


16. 30 


The optiffium value of K that nininizes the perfcrr.ar.ee 
index is obtained using a slightly modified form of the itera- 
tive procedure described in Appendix C. The modification is 
incorporated, to account for the decentralized control, as 
indicated below : 

In step 5 where K is updated according to 


j^{i+l) 


= K 


(i) 

~ ^ dK 



will be a full matrix even though is a block dia- 

gonal matrix. The block diagonal structure is imposed on 
K' by neglecting the elements that do not constitute dia- 
gonal blocks. 


When the solution for the optimum K is obtained in the 
final iteration, in the block diagonal form, the parameters of 
the PSS for the ith machine are obtained from 




where % i = the ith diagonal block. The elements of 
define the structure of PoS as given in eqns. (6.6) and (6,7). 
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6.4 N'UliERICAL HXA.V.?LE3 

rt a-rr.achine and a ISHSiachine power system exarples are 
considered tc illustrate the techniques, presented in 3ecs, 
6.2 and 6.3, for ?SS design. 


6.4,1 A Three-Machine Power Systen 

A sample multimachine povier system given in Ref, [83], 
consisting of nine buses and three generators, is considered. 
The single-line diagxan of the system alongwith system data 
and operating conditions are given in Appendix J. Each machine 
with its single time constant excitation system is represented 
by a fourth order model, (see Appendix D) . Machine 1 is con- 
sidered as reference for the measurement of rotor angle. Thus, 
an eleventh order model, describing the open-loop system, is 
obtained. Rotor speed is used as feedback signal to design a 


second order P3S with transfer function of the form 


PCs) 


®o® ®1® ®2 


(6.31) 


for each machine of the system. Open-loop system matrices are 

as follows : 

-0,13138 0 


m 11 


A„ 

m 


12 


1.83564 
-1.20741 0 0 

134.05135 0 -50. 

’ 0.02351 0 0,00<X55 

0.90183 0 11.81610 

-21.25201 0 20.98774 


0 


O 



lie 


C.C3342 
A^^13=' C. 35459 
L-3C. 96391 


0 0.05CS8 C"| 

i 

0 7.06217 0 1 

I 

0 -36.89218 oi 

.1 




^22 ~ 


^in21 ~ 


^1523 


-0.52193 0 

-10.95584 0 

0 1 
-155.30564 0 

0,00700 0 

—1 • 93066 0 

0 -1 
j-69. 22754 0 

““ 0.19452 C 

2.27457 0 

O 0 

— -79.79709 0 


-2.64893 

-65,97776 

0 

-67.39815 

0 : 

0 
0 ' 
o_ 

0.95153 

20.67036 

0 

-75,68095 


1.39121 

0 

0 

-50 


O' 

C 

0 

0- 


b. 


m2 

'022'* 


[O 

[0 


A 

m33 


\31 


— -0.62866 
-26.83562 
0 

^14.43926 
— 0,02403 
-1.49427 
O 

_r74. 62402 


0 -2.18978 

0 -119.75087 

1 0 

0 -321.13379 
0 0*' 

0 0 
■1 0 

O 0.-. .■ 



C 25CO] ^ 

1 c 3 


0 O 25CX)]^ 

10 0] 
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‘ C.I414G 

0 

0.97C74 

01 

-r.3 

C C 2500] “ 

. i 5.03265 

^32= 1 

0 

0 

52,64938 

0 

1 

C J 

1 

0 : 


ICC] ; 

'-54.39683 

0 

130,20166 



C„ . = 0 for 

Hi j “ 

i 

j- 





Table 6.1 shows the open-loop eigenvalues of the systen. 


Table 6.1 Open-Loop Eigenvalues of the 3-Machir.e 
System 

-0.34411 + j 8.C7210 
-0.81744 + jll. 51613 
-2.01443 
-3.78028 
-8.36729 
-41.40810 
-46.18024 
-47.20853 

0.0 

6.4. 1.1 Eigenvalue Assignment 

In the first stage of design, pole assignment technique 
based on modal control theory with partial state feedback, as 
given in 3ec. 2.4.1, is used to design stabilizers for each 



rrjachir.r, r.cgleccin^ interconnections. The tecnnique sceci- 
fies tr.e zeros of Pao. For the system considexeo, they are 
specified as -2.C and -2,5 for each PSS, Three eigenvalues 
for each machine can be assigned with this technique. 

Table 6.2 shows, at r * C, the open-loop eigenvalue^ idestred 
locations of three closed-loop eigenvalues and the remaining 
closed-loop eigenvalues obtained from the design. The desired 
locations of the closed-loop eigenvalues are obtained by 
defining the sector fea] given in Fig. 6.1. This ensures a 
damping ratio greater than a rainimura value, say I , where 
5 = cosQ , The values of a and 9 chosen are respectively 
-1 and 75 degrees. The real parts of complex eigenvalues are 
selected such that they lie on the boundary of the sector and 
the locations of real eigenvalues are specified inside the 
sector. 

Table 6.2 Eigenvalues of the 3-Machine System at r * O 



Machine 1 

Machine 2 

Machine 3 

Open-loop 

eigenvalues 

0.0 

-0.32150 -r 
j 7. 96260 "■ 

-0.93839 + 
jl0.8063C ” 

at r=0 

-5.68401 

-4.70536 

-3.13048 


-44.44736 

-45.17356 

-45.62081 

Desired location 
of 3-closed-loop 
eigenvalues 

-2,0 

-3.0 ,-4.0 

-2.14 + i&,QO 
-3.0 

-2.94 ± Jll.C 
-3.C 

Remaining 

closed-loop 

eigenvalues 

-1.16588 

-44.59163 

-6.87699 

-13.37328 

-48.13837 

-3.83532 

-33.76433 

-65.58955 










Table 6.3 shows the parameters of the three stabilizers 
obtained xn the first stage of design. 


Table 6.3 PS3 Parameters for the 3->=.achine dysten 
at r =0 


Machine 

9 

0 

h 

®2 


Y2 

1 

0.04504 

0.20267 

0.22518 

4.62613 

5.25227 

2 

0.09199 

0.41397 

0.45996 

25.14671 

67.47599 

3 

0.15194 

0.6S372 

0 .75969 

61 .44055 

196.59640 


The feedback gain matrix for ith machine can oe 
obtained in terms of ?S3 parameters, defined in eqn. (6.31) as 
( Jippendix B) : 

IQ. 


where 

and 


A 


•di 


io 

0 


9 * 

'^li 


0 1 
“”^21 "^lll 


© 


ii * ” ®oi 


©A- 

2r 


8 «, 9' Y 

2i oi ^2i 


In the second stage of design, for updating the decentra- 
lized feedback gain matrix Kj^, a step of &r = C.C3 is chosen 
to introduce interconnections. Zeros of the stabilizers are 
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kept fixe;, at the specified values in updating i.. TnU'S, 
the r.ur.hir of unknowns in eqn. (6.20) are nine. The or-oer cf 
ciosec-ioop systen> is seventeen and, therefore, eqn. (6.2C) 
represents a set of seventeen sinultaneous linear algebraic 
equations in nine unknowns. Pseudo-inverse of a 

dK- 

17 X 9 .latrix, is determined for obtaining a solution of >--■ 
in eqn. (6.20). Table 6.4 shows the final values of PSS 
parai.ieters with full interconnection, that is, with r = 1. 





Table 6.5 Closed-Loop Eigenvalues of the 3-'.ac^ir,-e 
Sys tern 


.%ssigned at 
r=G 

;;ith PS3 of 

Table 6.3 

Obtained with ?S3 
of Table 6.3 at 
r » 1 

Obtainsc! '.vith ?33 
of Table 6.4 at 
r * 1 

-2.14 + j 3.00 

-1.16640 ± j 3.03543 

-2.11541 + j 7.S64C7 

—2*94 -r jll#(X) 

-2.78779 + j 12 .414^ 

-2.92745 + j 12 ,26552 

- 1.16588 

- 0.09647 

- 0.78921 

- 2.0CXXX3 

- 2.00000 

- 2.00839 

- 3.00000 

- 2.M994 

- 2.12231 

- 3.00000 

- 2.69879 

- 2.59251 

- 3.00000 

- 3.05555 

- 3.23480 

- 3.83530 

- 3.45737 

- 3.70848 

- 4.00000 

- 3.70396 

- 5.46943 

- 6.87702 
-13.37271 

-12.62314 + j 3.41233 

-13.95100 + j 4.98277 

-33.76433 

-33.26061 

-32.95496 

-44.59163 

^44.49973 

-44.62495 

-48.13835 

-48.21465 

-48 .70835 

-63.58957 

-65.70617 

-65.79412 
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A c orTi"} arisen of Tables 6*3 and 6*4 shows that the r-.odifi- 
catior.s in ?SS parameters with interconnections are not signi- 
ficant except for the poles of P3S of machines 2 and 3. The 
poles of ?SS for the second machine have changed significantly 
by considering the interconnections # 

Aesults shows in Table 6.5 emphasize the need for coordi- 
nated design of PS3 in multiiaachine system. The second column 
of Table 6.3 gives the eigenvalues of the system with ?3S 
designed by neglecting the interconnection. It is seen that 
one of the real eigenvalue is close to the imaginary axis and 
the dandling of the rotor oscillation mooes are reduced. 

It is to be noted that eigenvalues given in coluMn 3 are 
not identical to those given in coiWiH 1, although the design 
algorithm is aimed at maintaining the system eigenvalues at 
the locations given in colUiOn 1 in the presence of interconne- 
ctions. However, this is not surprizing as only nine para- 
meters (of the PS3) are used to control seventeen eigenvclucs. 
.Notwithstanding the above, the results are quite encouraging 
in that the eigenvalues of column 3 are close enough to those 
of column 1. 

6. 4. 1.2 ration of Performance Index 

The objective is to rainiaiize the performance index given 
in eqn . (6.21) with 

“1 = 1 


for i .* 1,2 and 3 



by designing the PSS for each machine. For ccnv2riienc3, tr.e 
initial choice of the PSS parameters is made fro.n ihe res -Its 
obtained in the previous section, using the eigenvai-e assign- 
ment technique. The optimal values of P3.> parameters are 
given in Table 6.6. Table 6.7 shows the closed-icop eigen-- 
values obtained with stabilizers of Table 6.6. 

Table 6.6 PS3 Parameters fox the 3-r.iachine System 
Obtained from Parameter Optimization 
Technique 


iiochine 

0.. 

0 

S- 

1 

®2 


^2 

1 

O.C5973 

0.25637 

0.29749 

4.58555 

5.19168 

2 

0.16013 

2.33035 

9.09180 

30.56508 

148.62720 

3 

0.16731 

1.33959 

2.77806 

61.20167 

188.05472 


Table 6.7 Eigenvalues of the 3-Machine System 
with PSS of Table 6.6 


-2,41665 -*■ j 9.30942 
-2.45847 + jl3. 68310 
-1,CX)378 ** 

-1.48954 
-1.80292 
-2.12409 
-2.96189 
-3.91336 
—6 .45344 
-11.18142 
-15.24140 
-33.53356 
-44.39199 
-49.45159 
-65.85462 
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The performance of the stabilizers, obtained from the 
two design techniques presented, is compared by taking the 
tisi'e-- response of fluctuations in power output and relative 
rotor angle of the machines following a disturbance. 

Figures 6.2 to 6.7 show the variations in power outputs of 
the machines and the relative rotor angles for the three 
cases, namely, for a unit step change in the reference 
voltage of machines 1,2 and 3. 

Response curves shown in Figs. 6.2 to 6.4 indicate the 
significantly small power fluctuation obtained with P33 of 
Table 6.6, that is, with stabilizers designed using the 
parameter optimization technique. The overshoot and 
settling time of these fluctuations are considerably 
smaller than those obtained with the stabilizers designed 
on eigenvalue-assignment basis. This is as expected since 
the objective of the stabilizer designed in the former case 
is the mininization of fluctuation in power outputs. It is 
interesting to observe that the relative rotor angles are 
well damped (with no oscillations) for the case of stabi- 
lizers designed on the basis of minimization of this per- 
formance index. 

A comparison of Tables 6.4 and 6.6 shows the parameters 
of the stabilizers designed using the two approaches differ 
meinly for machines 2 and 3. Sven here, the change is mainly 
in the gain and zeros of the P3S. it vtos mentioned earlier 
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tl'iat parar.eters of PSS fox the machine 1 are net such 
affected by the interconnection. Thus it '.‘/ouid appear test 
satisfactory design of PSS is required in this systeri for 
riachines 2 and 3, taking into account the interconnections 
between machines. 

Reference [43] also considers the design of PSS for a 
three machine system (different than what is given here) 
based on the technique given in [77 J. However, unlike the 
approach given here (of trying to control all the eigenvalues) 
the author of [48] has considered the assignment of only six 
eigenvalues using nine parameters of P3S. The work reported 
in this chapter was carried out independently and published 
[ 79 ] before the author came to know of the work reported in 

[ 483 * 

6.4.2 A Thirteen Machine Power System 

The case study of a 13~raachine power system described in 
Chapter 5 is taken up to illustrate the application of the 
design techniques, presented in Secs. 6.2 and 6.3, for large 
scale systems. The analysis presented, in Chapter 5, for the 
choice of P>S locations suggests that by equipping on 
three machines, namely, machines 1,7 and 11, nost of the reter 
oscillation modes could be adequately damped, stabilizers 
are, therefore, designed for machines 1,7 and 11. These 
machines ere represented by fourth order medeis. The remain- 
ing iaachines are represented by the classical models- The ■ 
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crsSiiX zf the cpen-locp and the closed-icop systt.r.s are, 
rc-spectiveiy, 31 and 37, Eigenvalues cf the oter.-lccp syste.-* 
are- given in the third colisan of Table 5.1 in -Uhaptex 5, The 
systera is dynamically unstable without stabilizers. 

6. 4. 2,1 Eigenvalue Assiganient 

As in case of the 3~»iiachine system, stabilizers arc 
first designed for the selected machines, using the algorithm 
given in 3ec. 2,4.1, by neglecting the interconnections. Zeros 
of each PSS are specified as -2.C and -2,5, Table 6.8 shovis, 
at r = C, the open-loop eigenvalues, desired locations of 
three closed-loop eigenvalues and remaining eigenvalues of 
the closed-loop system obtained from the design. Parameters 
of the stabilizers obtained in the first stage of design are 
given in Table 6.9. 


Table 6.8 Eigenvalues of the Selected Machines 
in the 13-Machine System at r - 0 



Miachine 1 

Machine 7 

Machine 11 

Open-Loop 
eigenvalues 
at r « 0 

0.08082 + 
j 9.07955 ■" 

-5.96516 

-44,64544 

0.07421 + 
j 13. 16434 “* 

-6.65288 

0.10407 4 
j8.17£27 

-6.16746 

-44.41907 

Desired Loca- 
tions of 

3 Closed-Loop 
Eigenvalue 


-2.41 + j9.ee 

-3.00 

Remaining 

-6,77244 + 

-12.93iCX) + 

-5,33129 * 

Closed-Loop 

j 4. 88345 

Jl. 68711 

j4.9C926 

Eigenvalues 

-47,54431 

-57,30393 -46.55321 












2C4 


Ts’-ile 6.9 PiS Purassc-ters for the Seioctoy 

in the 13--i.achine Systam at r - C 

ii'achine 9 

o 

i 0.29135 1.31110 1.45677 19.00323 4£.i4C2C 

7 0.34159 1,53715 l,7C794 42,76278 119.96623 

11 0.18836 C. 84761 C.94178 15.65740 33.43923 


In the second stage of design, a step of 6r * 0.G5 is 
considered in updating the gain and the poles of the PSS, 

Table 6. 1C shows the final values of ?o3 paraueters with full 
interconnection. Table 6.11 shows the closed-loop eigenvalues 
of the systen. First coluinn of Table 6.11 shows the eigenvalues 
of the decoupled system with P33 parameters of Table 6.9, that 
is, with stabilizers designed in the first stage of design. 
Second and third coliaans of Table 6.11 shows the eigenvalues of 
the fully connected system obtained with P33 parameters of 
Tables 6.9 and 6.10 resoectively. 


Table 6.10 PSS Parameters for the Selected Machines 
in the IS-Machine 3ysten at r » 1 


Machine 

8 

5 

2 

^1 

^2 

1 

C. 23609 

1.28739 

1.43044 

24.94678 

119.56558 

7 

0.32659 

1.46966 

1.63296 

42.51388 

119,47225 

11 

0.17572 

C.7%76 

0.87S62 

17.17517 

51.16336 
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Table 6,11 Closed-Loop Eigenvalues of the 1 e 

3y£ ;era 


Assigned at r = 0 
with P3S of Table 6.10 

-3.75cx:c jK.ccxicc 
C.C .-r jl3.C141C 
C.C + jl0.501C4 
C.O + jlO. 24194 
C.O + jlO. 18118 
-2,680CX3 + jlO.OCXXX) 
C.O + j 9,86085 
0.0 + j 9,68993 
0.0 + j 9.34472 

- 2.41000 + j 9.0(XXX) 

0.0 + j 8.13542 
0.0 + j 7.31242 
0.0000 
-3,CXj(X5 
-3.0C00 
-3.0000 

- 5.83123 + j 4.90926 

- 6.77243 + j 4.88344 
-12.93102 + j 1.68753 

-46.55321 

-47.54433 

-57.30651 


Obtained at r » 1 
with PCS of 
Table 6.10 

- 3.C839C + J 14. 12752 

- 0.18909 ± j 13- 393% 

- 0.40267 + jii. 79212 

- 0.40256 + j 10.86772 

- 0.C4167 + jlC.SCX3 

- 0.03537 ± j 10.51028 

- 0.68569 t j 10.47896 

- 0.03828 ± j 10. 27158 

- 0.03207 + j 9,41924 

- 0.61724 + j 8.97292 

- 0.53791 .t j 8.00539 

- 0.65133 + j 6.43692 

- 0.620C8 

- 2,89078 

- 2.99349 

- 3.CX)116 

- 6.29856 + j 5.17551 

- 9.00956 ± j 7.47702 
-10,49361 
-13.65274 
-46.58312 
-47,54057 
-57.32446 


Cbtainei at r = 1 
v.ith ?Ji c-f 
"able 6.11 

- 2.92532 + jl4.iC79C 

- 0,13524 + jU.j7CS2 

- 0.2X79 j 11.77149 

- 0.13630 + jlC.7C338 

- C. 06316 z jlC.566S7 

- C.C68C7 + jlC, 50823 

- 0.57632 j 10. 40899 

- C. 05326 7. j 10, 26468 

- O.C2639 7 j 9.41090 

- 0.43355 + j 8.77155 

- C. 43318 7 j 7.90272 

- O. 55035 7 J 6.17615 

- 0.38637 

- 2.98356 

- 3.5C351 

- 4.53175 

- 7.22301 + j 5.86688 
-10.970<% 7 j 8,82632 
-10.19024 
—15.53490 
-46.52583 
-47.92226 
-56.93663 
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A ccRcarison of Tables 6.9 and 6. 1C that tho 

ncdificaticn in the pararaeters of P3o are net significant fer 
machine 7, The first column of the Table 5,11 sno"s tnat 
system is critically stable as the real parts of the eigen- 
values of the rotor oscillation modes are zero for machines 
represented by the classical model. If the eigenvalues are 
to be retained at these locations, as per the objective of th~ 
design algorithm, the design will not be acceptable. However, 
the second and third columns show that the system is stable 
with the PSS designed for the three selected machines, al- 
though the damping of some of the rotor oscillation modes are 
small* 

6. 4, 2. 2 Minimization of Perfomance Index 

As in case of the 3-machine system, parameter optimiza- 
tion technique, of Sec, 6*3, is used to design PSS for the 
three selected machines, that is machines 1,7 and 11, The 
objective is to minimize the pei'formance index given in eqn, 
(6.21) with the weightage parameter defined as follows i 
for i * 1,7 and 11 

«! = i 

1^0 for all other values of i 

For the sake of convenience, the initial choice of PoS para- 
meters is made from the results obtained in the previous 
section. The optimal values of PSS parameters are given in' 
Table 6.12 and the coritesfsonding closed-loop eiger.vr- lues are 
given in Table 6.13, 
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Table 5.12 FJ3 Parameters for the Jelecteb -Viachir.cs 
in the i3-?!achlne System Cbtairei frcn 
Parameter Cptimization Technique 


'/.achine 



9 

^^2 


^2 

1 

C. 37164 

2,70569 

12.05420 

25.C4285 

119.54951 

7 

0.32432 

1.19484 

1.45^1 

^2,ABBC7 

119.43035 

11 

0.42671 

4,87906 

13.33516 

17.16283 

5i.l6X'S 


Table 6.13 Eigenvalues of the i3-4lachine ayster. with 
P33 of Table 6.12 


-2,89734 + 
-1.61382 + 
-0.1906C ± 
-0.37997 t 
—0.35992 + 
-0.34208 ± 
-0.12740 + 
-0.53875 i 
-0.22389 i 
-0.35956 + 
-0.28609 + 
-0.27365 ± 
-10.76684 X 
-1.81087 ± 
-12.81823 + 
-2.85494 
-3.11521 
-5.86962 
-7.10297 
-48.46174 -I- 

m-m 

-56. 90^^39 


j 13. 80945 
jl3. 61650 
j 13. 39934 
j 11. 79321 
j 10. 86744 
j 10. 59898 
jlC. 51131 
j 10. 26825 
j 9.44422 
j 9.39171 
j 8.65416 
j 6.90138 
j 9.05954 
j 2.84426 
j 2.07492 


J'. 0.071491 


*-ht two 


The effectiveness of stabilizers, designed jitn 
techniques presented, is corr.pared by taking tice respcnsr of 
fluctuations in power output and the relative xctor ar.rle of 
the selected generators following a di sturbarxe. Three caoes 
are considered : 

(1) a unit step change in reference voltage of machine 1. 

(2) a unit step change in reference voltage of raach.».ne 7. 

(3) a unit step change in reference voltage of machine 11, 

The variation in the power output and the relative rotor 
angles for machines 1,7 and 11, for all the three cases are 
given in Figs. 6.8 to 6,13. 

From these figures it is clear that the PS3 design based 
on minimization of performance index results in superior pex>- 
formance particularly in terms of variation in the generator 
outputs. Hot only the magnitudes pf oscillations are negli- 
gible, but the settling time is also reduced, 

A comparison of Tables 6.10 and 6.12 shews that the 
poles of the stabilizers, designed using the two techniques, 
remain almost unchanged. The change is mainly in gain and 
zeros of stabilizers of machine 1 and 11. 

6.4,3 Oiscussions 

A comparison of the two approaches for PiS design can be 
made based cn the results of the two examples presented. It is 
seen that t.;e method based on the ainiuiization of performance 
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index (the fluctuations in the power out'^ut ; cf -vi*; ‘»el'.-ctc3 
generators) gives excellent results. The a.'prcacn Pasrd on 
eigenvalue assignment (utilizing the algorithm giver in [77 j| 
is not so satisfactory in cooparisoAt particularly for the 13- 
machine system. 

An n machine power system, without an infinite bus, hns 
(n-l) modes of rotor oscillations. For large systems many 0; 
these modes are quite close. This can be seen from the 
results of Table 5,1. The objectives of the stabilizers in a 
system are primarily to damp certain modes of rotor oscilia~ 
tions which can be troublesome. Analytically this can be 
found from the inspection of the open-loop eigenvalues, 
utilizing suitable system model. Hcmever, most often the pro- 
blems are encountered during the system operation under various 
conditions where oscillations of power in transmission lines 
are experienced. It can be quite difficult to identify the 
particular modes that give rise to these oscillations in a 
large system. 

Even if the problerj of identification of unstable mode# 
can be solved analytically, it is difficult to associate a 
particular mode vdth a particular machine in all the cases. 

This is more so in a large system as deroonstxeted for the 
13-machine system in Chapter 5, Because of this problem, 
it is difficult to specify the eigenvalues when the inter- 
connections are neglected (as required in the first stage 









cf tha algorithm of [77]). This problem is oggr.^vAted if ro.5 
IS to be provided on all the machines as ono ( ret'erer.ce j ma- 
chine will be left out which has no rotor mods associated 
■..'ith it in the absence of interconnections (see th® 3-r,achine 
example given above). Ref. [48] considers a 3-r:achlne oxa-r,!®} 
with an infinite bus which has three modes} the incluox(,r, of 
infinite bus, however, is not satisfactory In general. 

The above discussion shews that there is a problem asso- 
ciated with the design approach based on the algorithm given in 
[77]. Moreover, additional problems cropup in the second 
stage of the algorithm where the controller parameters are 
modified. These are (1) inadequacy of limited controller 
parameters in influencing all the dominant eigenvalues and 
(2) numerical problems caused by the drift of eigenvalues v^hich 
make it impossible to distinguish between different eigenvalues 
ivhen they are quite close. 

As the power system stabilizers are designed to overcoE'.e 
the problem of low frequency poi>?er oscillations, it is con- 
venient to choose the objective of minimizing the power flu- 
ctuations. Although, the niaaber of examples considered axe not 
extensive, the results based on this approach are excellent. 
However, the parameter optimization technique requires large- 
computational effort in comparison with eigenvalue as5inr.r.:ant 
method. The minimization of performance index is an iterative 
procedure tidiich requires the solution of two Liapunov equations 
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in each iteration. The tolutlon of Liapunov equation itself 
is obtained through an iterative technique [50]. This 
increases the coniputational effort of this .Tiethod in conpari- 
son with the method based on eigenvalue assignaent. The 
computational effort of parameter optimization technique is 
about six to ten times, in terms of computer time, than 
required in the eigenvalue assignment technique* Sine© the 
design is done off-line, the coiiputationai burden involved 
does not pose a serious limitation - 

6.5 OONCLIBIOIIS 

The detailed investigation into the design of PSS in a > 
large multimachine system using docentralized control is i 

presented in this chapter. A co^arative study of the two j 
approaches, one based on eigenvalue assignment using the i 

algorithm given in [77] , the second based on ainiosization of | 

performance index is carried out using two numerical examples. I 

( 

The case study of a IS-machine system presented here is typi- I 
cal of a largo power system writ^re PSS are designod on selected | 
machines. i 

The results of the case study indicate the superiority i 
of the design technique based on paraiaotcr optimization with | 

■ f 

the objective of minimizing tlie fluctuations in the power I 

output of generators. Although the cofaputationai burden is | 

I 

greater in this isethod, the results justify its adoption | 

for largo scale systems . The inherent ireaknossos of tho | 

algorithm given in [77] for the design of PSS make it I 

unattractive for its wHeetion in large systoos. I 





CHAPTER 7 


OC»ICLUSiaB 


7.1 GENERAL 

This thesis has dealt with several aspects in design 
of power system stabilizers, naiaely, (i) choice of design 
objectives and algoritlms for PSS design as dynaoic coc^enia- 
tors, (ii) stiKiy of adequacy of system wjdels and effective- 
ness of control signals and (iii) selection of the location 
and design of PSS in large multiraachine power systesas usii^ 
decentralized control tiwoiry. Specific contributions of the 
thesis are reviewed below. 

7.2 OBJECTIVES iWD TECW^IC^S FOR PSS DESI»I 

The design of PSS in power iixtustry has been based on 
the objective of la^roving the d^Bpin^ torque. The design 
techniques are derived from classical control theory and 
utilize froqiasncy domain techniques. The drawbacks in 
this approach are (i) design freedom is not fully utilized 
and (ii) it is difficult to exteml it to multinachlne system 
model. With this in view* desi^p techniques have boon dovo- 
loped in Chapter 2 based c«i the objectives of polo assigrawint 
and minimization of a perfoiaaiK:© Index. 
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Pole assignment is a natural objective in power system 
stabilization. It has been shown [44] that the improveaent 
of camping torque is related to the shifting of eigenvalues 
corresponding to the rotor oscillations to the left in the 



design freedom cm be utilized. The oxamplos given for the 
illustration show that, the transfer function of PSS in aiany 



mine the PSS parai^tors* 







220 


Assignment of closedoioop poles anywhere in a sector in 
the left half of the complex plane guarantees stabilization 
of the system with a prescribed degree of damping. However t 
the design algorithms require exact specif ication of the 
closed-loop poles. This, in general, is difficult and 
arbitrary location of poles may not result always in a 



earlier, to relate the performance index to an objective 
related to the dating of oscillations used by power engine rs. 
The work reported in Chapter 2 overcomes both these limitations 
and is novel in terms of PSS design. Although, conputationai 
roquireaient is more in this case, the system response folioi*- 
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7.3 CHOICE OF MACHINE MODEL AND EFFECTI'/BiESS OF CWITROL 
SIGTiALS 


The siniplified model of synchronous nachins 1 neglect- 
ing the effect of aaortisseur windings) has been extensively 
used for the design of PSS, howevexi it* adequacy has not 
been analysed. The adequacy of the simolified machine 



is adequate for PSS design* 
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Chapter 4 has presented a coaparative study cf the 
control signals for the design of PSS with the objective 
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(corresponclifig to th« ntxlt of otciilations bot.'^etn tht tw3 
nachines) is negatively dansped as the gain Is increased. It 
was experir^ntally observed that the use of average speed 
did not have such limitations » however, no explanation it 
given in {25] to explain this pbrnm^m* It is shorn in 
/gspeiidix F that for identical machines and with identical 
loading, the local mode is not controllable when average 
speed is used. Thus, it is not surprising that increase 
in PSS, gain has no effect on local mode and hence the un- 
damping of local mode is not a factor in limiting th© gain, 

7 .4 PSS DESIGN IN MH-TIMACmNE 

Coordinated design of PSS in multioachine system has 
gained much attention recently. Design of PSS in multlmachin® 
systems is not straightforward i^causo of the co«piexities 
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ninimization of the perforraam:e Index. Decentralized con- 
trol theory has been used for the design of stabilizers. 

Two examples » a 3 machine and a 13 machine systems have 
been considered. In case stabilizers are designed for 
limited number of machines, it is found that the pole assign- 
ment technique given in [77] is not effective. The design 
method based on the minimization of the performance index 
(related to the minimization of power fluctuations) is found 
to superior as it gives »«11 danped responses with little 
overshoot. 

7.5 SUaZSIIOMS FOR FURTHER 

Methods presented in this thesis for the design of ?SS 


s 


- 1 
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APPS?I01X A 


CXi'.iFai3AT0H FwliCTIC-^i 

A.l CCM?2:;3AT0a OF SHCTiaN 2.3.1 

Consider the single inout linear cvnaff’.icdl syste'-'i 

P X “ [a] X b u (A.i) 

where x c 'pP and u eftl. Let x' be the vector of accessible 
states. Then 

x/= [c] X Ca.2) 

where x* c r < n and C is an r x n aatrix yliich can be 
expressed as 

c =• [0)lj,] 

being an r x r Identity matrix and 0 an r x (n-r) null 

matrix. 

Consider a second order dynamic corapensator 

p _z ss [DJ ,2 + e u CA.3) 

v/here 0 is a diagonal matrix and e Is a vector v/ith each element 
of unity. Combining eqns. (A.l) and (A. 3), 40 gc-t 






Eqn. (a. 4) represents the augronted syster,. Fcr the control 
law 

1 {k,5) 

the closed-loop system is 

p X » (S + I X (A.6) 

where 

=* [lSiC|: 1^] 

The cOTpensator transfer fynction can be obtained as 
follows. From eqns. (A. 3) and (A. 5) the state variables of 
comoensator can be obtained as 

2 * Csl - ® #cj x»(s) CA.7) 

where 

[pj =|p 4- e ^3 

Substituting the value of 2 from eqn. (A.7) into (A.o] , we 
get 





Uls) = [1 k|(sI ~ ej k5 xHs; 

For a second-order coLpcnsatcr 
“d 0 

i'*! aaw# ' 

ijf ■*** 

-0 d, 

i 

Substituting the values of e and in eqn. (A. 3), got 

F(s) (s-dj) (s-dj) kJ 


and 


.,t 






where 

U(s) = F(s) x’(s) 

and 

, 2 

^gVs} — s — ^ 1'^^^22 ^ ^ ^^1^2 

^1^^S22'*'^2*''21^ Ca«9) 

Thus, zeros of F(s), the coa^jensator transfer function, are 
eigenvalues of D. Siciiarly, for higher order dynamic cora- 
pensators also, it can be shown that zeros of the compensator 
are the eigenvalues of D. 

A *2 ca.iPE:-JSATCa O? SECIIOfl 2.3.2 

The transfer function of a single input ^th order dynamic 
compensator as given in Sec. 2.3.2 is 




where y and u are input and output of tm compensator and 

... 0 

... 0 

• 

♦ 

* 

0 1 

1 

Q = [9^ ®k-I •** ^ ® 

and 

* S|^ — Bq 

consider the single input sii^lc- output linear systen as 

p X = [ Aj ^ + b u ( A * 13 ) 




y =* [c] X 


(A. 14) 
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Co^bininc eqr.s . (A. 11} to (A.l'), th-e 

and control law of the aug.nented systec {systen anr cc.'csnsa- 
tor) can be given as 


P 1 = i.nj X 4* [jj u 


A f A 

1 = X 


and 


where 


A A 

u = lAj y 




^ -..It: 

( -. 16 } 

i.‘..i7) 



A 

K is the gain output feedback :tatrix of the augraented systeiai 
given by eqns. (A. 15) to (A.17) . For a second orcer con^ensa- 

A 

tor, matrix K will take the foxi.'. 



1 



?CL£ .;33ia-;..£NT ^ITK OUTPUT F£iD-^U:< fef. 


The t3chnic;ue given here is based on Vt'fy »iorX cf ,'.>-rro :,n , 
Hirbod fc6 ] . The n^tnod presenta.* in [66 ^ orovidss a 
for the design of full-tank co.nocnsators for .t.u 1 tivarif^- : > 
linear systems. Here the technique is restricted to single 
input linear systems 


and 


3 X = [a] X r b y 

X = [C] X 


C.l) 


iB.2) 


where, x e R^, u £ R ^ and ye?/ are respectively vectors of 
state, input and output variables. The r x 1 open-loop trans- 
fer function matrix G(s) is 

[Nj^CsT] 


G(s) = C(sl - a)"^ b 


i s) 


NJs) 


(b.3l 


where 

Aq(s) =» ^ ... * \ 

(3.4) 

and 

N^is) S' * ^12®” ^ -t ••• 4- 

(L.3) 

^(s) 

is the characteristic poiyncKaial of G(s). 

Censioer t*"'. 

output feedback law 




uCs) = u^(s) - F{s) . y(s) 

where c represents the reference input. 

The problem can be define:’, as the csteminaticn cf 1 4 r 
dynamic feedback conoensator, F^s), such that the clcsed-lccn 
system 

H(s) [I t G(s) F(s)]“^ GCs) Cfa-") 

has a desired set of eigenvalues# The k.th order compensator 
F(s) has the form 

F(s) ” MyCs)! 19 #8) 

0 

where 

A^(s) = ^2*^"*^ ♦ ,.# HI- Yjt CB#9) 

and u y % 

* ®io®* + ... + ®|ic (3.10) 

A (s) is the characteristic polynomial of coK^ensitor transfer 
function F(s). 

The resulting closed-loop system characteristic polynomial 
A^(s), defined as : 

m f*'**'^ d, (o.llj 

d ' * 

c an be v/ritten as [66 ] 



(5.12) 





3.1 Cu.i?L£T3 POL! .ASSISI £;;! 

Tha prcblea of pole assign:.ier.t can oe csfina'- 
^3(5)1 and a desired set of clust'l'i” tp r.Di':s, 

find A^(s), M^{s), of lowest cegree t'^Lzr f 

eqa. (3.12). equating coefficients of like poviars of s if. 
eqn. (B.12), we get 



k colurms Ie-*'! colwffis ^+•1 columns 



14 ; 


P-ic r2 ’•* '^ki Ko 

and 

= CCdj^-a^) 1 (d^-^)*... 



;\c *■*• 







jiquation {B.13) is a set of (n*»-k) equations in :.^r- 

known parameters of FCs). The difference vector 6. cont-iirs 
the coefficients of the polynomial 


A^(s) = ^(s) - ^g{s). 

A necessary and sufficient condition for the existence of solu- 
tion of eqn, (B.13) is 

Bank I'A^ ] * Rank 6j^3 (3.15) 

and a solution of eqn. (3.13) is 

JEk ** t^k^3 ■ (B.i6) 

where is the generalized inverse of such that 

[Xjj] l\} * i\} 

B,2 PARTI*^ POLE ASSI®£.iS?iT 

Consider the case urtien q poles are to be assigned, 
t = (n -r k - q) poles art ailonaid to assur*o arbitrary loca- 
tions, For this case the clCNied-loop system characteristic 
polynomial i^Cs) has tlwi fora 



Bk - I2 


• • • 
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and 




”k ~ { {d2--G:2 ^ J • • • 


IS ^ I ,, _■, 

^ikJ :-rc ••• \.<* 


^quation (B.13) is a set of (n-^k) equations in Likt-l;r»':j gn- 
known parameters of F(s)* The difference vector 6,^ contains 
the coefficients of the polynomial 


AJjCs) = iy(s) - ^qCs). (B.14) 

A necessary and sufficient condition for the existence of solu- 
tion of eqn* (B«13) is 

Rank [■*^ 1 ^ J * Rank Cji*15} 

and a solution of eqn. (3*13) is 

£k = (a. 16 ) 

where xf Is the generalited inverse of X,, such that 

[X^] [xg^] [X^] =■ [Xj^] 

Be 2 PrlRTI.U. POLE .ISSIC2C.13NT 

Consider the case when q poles are to be assigned, 
t = (n + k - q) poles art allmutd to assume arbitrar/ loca- 
tions, For this cast th« ciostd-loop system characteristic 
polynomial l^Cs) has ttwi for® 



A^^(s) = ( + ••• + cf^) ( - ... ^ e/, 

ih.l'^j 

where d 2 , ...» d^ are specified and ^ 2 . ...» t‘-^ 

be determined. Here the difference vector can '-'e coiaAnoG 

by equating the coefficients of like powers of s in iC.17) and 
(3.12) as 

6 . = ^ Oi ^ •*• ,»• + e*. (B.18) 

where the vectors 6^ and (i » 1# ...» t) contain res- 

pectively the coofficlents of the oolynomials 

A^(s) =* A^(s)*s'^ — A^(s).s^ (»3.19) 

and 

Al' (s ) ** A { s)«s'^^» (i*lf,,»»t) (B .20} 

* 

Hqn. (B,13) in this case will take the form 

[X|,] Ek = * [0] s (B.21) 


v/here 



.22: 


Eqn. (B.21) can bo rearranged as 

£k = 5k 

where 

Cxy* [x^, -0] 

and (i.23) 

• - t-it 

= LP}j> e i 

Hqn. (B.22) is a set of (n+k) siaultaneous algebraic equations 
in Cr(k+1) + k + t] unknowns. Vector contains parameters of 
F(s) and the coefficients of polynomial of unassigned closed- 
loop poles. The necessary and sufficient condition for the 
existence of solution of eqn. (B.22) is given in eqn. (D.IS) 
with Xjj> and replaced by X|j,, and respectively. 
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Consider the linear, tisse Invariant, conipietciy contrc- 
llable and observable system 

P X =s [a] X + [a] u Cc.i) 

y = [cj X ic.2) 

where x cr'^, u e R® and ycR^ are respectively the state, 
control and output variables. The problem Is to design a 
controller using output feedback to miniralze the quadratic 
cost function 

J ^ I H*" ^ CC.3) 

where Q is a constant n x n positive semldef inite matrix 
and R is a constant rax® positive definite matrix. The 
control law is 

us- [kJ J (C.4) 

The closed-loop systeo with Idle control law of eqrs. car. 

be given as 

P ^ I ■ 

where 

* A •• BISC 



From eqns, (C.l), (C.2), (C.3) and (C.^), get 


J = / (Q + cVrkc) X dt ic.e; 

o 

As the solution of oqn. (C.5) depends upon ini t ini st.ste 
x(o), the cost J given by eqn. (C.6) will be a function 01 
x(o). For the completely controllable and observable system 
of eqns. (C.l) and (C.2) , the cost J is given from |C.5) and 
(C.6) by [53] 


J = x^(o) S x(o) (C.7) 

where the cost matrix S is given as 

A^ A, A 

S = / e ^ (Q + C^K%<C) e ^ dt (C.8) 

o 

S can be obtained as the unique positive definite solution of 
the Liapunov equation fes ] 

aJs + S\ + (Q + C^K^aKC) » 0 (C.9) 

If the initial state xCo) is assumed to be a random 
variable unifomly distributed over an n-dimertsional unit 
sphere, the average cost cm be sem to be fc3 ] 


J » trace CsJ 


(cac) 



^ forT.ing the Kanil Ionian corresponding to sms. e-: 

(C.IO), we get 

H = trace[>] + trace[P^(A^3 ^ SA^ ♦ 0 T<ifO j (C.ll; 

where P is an n x n matrix of Langrange multipliers. The 
necessary conditions for iiinimum J [54] ar® 

•g^ = 0 ; g-g = 0 and ® (C.12) 

The first condition in (C.12) gives eqn. CC.9). The second 
condition in {C.12) gives 

PaJ + A^P + I = 0 {C.13) 

where I is an n x n identity aatrix. The third condition in 
eqn. (C.12) gives 

^ = 2(HKCPC^ - B^SPC^) = 0 CC.14) 

The gradient obtained from eqn. (C.14) Is utilized in 
applying steepest descent riwthod for i^dating the feedback 
matrix K. Following are the steps in getting the optiaun 
value of K. 

Step 1 : Start with initial guess of ^ That 

A - SK^o^C is stable. 

Step 2 : Compute S by solving Liapunov eqn. (C.9) 



step 3 
Step 4 
Step 5 


Step 6 


: Comoute ? by solving Liapunov eqn. C*::.!:*/ 
: Compute gradient from ecn. tC.14) 

: Update K according to 






where 0 < a < i 

: Go to step 2 until iy|^i is close to zero or 
difference in the values of J obtained in two 
successive iteration is less thari a prespecified 


value . 
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SYNCHRONOUS MCHINE E'^'O^^TICS^S AT^D EXCIT’illOt: 

SYSTEM 


D .l SYNCHRONOUS MCHIN3 
D.i.l Hi#ier Order Mwiel 

Machine armature is represanted by a singlt -phast 
equivalent circuit [78] show in Fig* 0.1. Tha equivalent 
circuit consists of a dependent current source in con junction 
with a constant impedance. The magnitude and phase angle of 
the dependent current source are determined from the knoviledge 
of the rotor flux linkages and tr» rotor angle. The model 
can take into account any number of circuits in the d and q 
axes of the rotor depemiing upon the details required and the 
availability of the machine data. Hov»ver, a synchronous 
machine ’.vith two rotor circuits in each axis is considered as 
shown in Fig. U.2. Tiie subscripts f and h refer to the field 
circuit and the damper circuit in the d-axis respectively. The 
subscripts g and k refer to the daaper circuits in the q-axis. 
Ml variables and parameters are expressed in per unit of 
the machine rating. Unlike in the Ref. the generator 

convention is used. Also, the convention of q-axis leadin- 

the d-axis is folloi«ed here. 




FIG. 0.2 



R otor Circuit equations 


Park's equations for rotor circuit are glvw cv 

Vr * t'^rJ ir ® 

o 

where 

y-r - 0 0 Oj^ 

4 — r^tjtTt 

~r “ ‘■-rd -rq^ 

iLl^ 

— r - rd ~rq 

= Block diag [ I I » f 

~rd ~ * --rq * 

»r^ = [»f ! Sr, - t^lc *9^* 

and [R^3 = Diag [Rf, 5 l\^i * Dlag [R,^. RgJ 

The flux linkages in the machine are defined tjy ; 


mmtm mmm. 

Is 


[Lgl U-rf 

m 

-Jr- 



Is 

- [*d 

I 

q** 

“S 


« I’fe 
“•i«J 


w Diaf [L^i* ^^3 


ia. 2 ) 


where 





[k] = Llock diag. 

[L^] = Block diag, CCLj.^1, t I 

= [M^ 

== [M, 

[‘■rd^ = 


Mj 

-*<1 




and 


[L. ] = 

rq 



‘‘‘dh-* 

^qk 



hh 

^fh 



^kg 

hg 

L 

9 — 


The rotor flux linkages are more important than the rotor 
currents as explained in CT3 ] • Therefore, eliminating and 
i from eqns. {D.l) and (0.2), we get 


P ird “ ^ -rd 

P Irq “ ^rq * ^ 


Stator Zquati ons. 

Park’s equations f«r the stator are given by s 


(0,3) 

(0.4) 
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The armature transients are, ignortd by neglecting ; '4, 

terras in Park's eqn. (u.5) . This is norrallv jv-itif 
since uhese terms are small cockered to the speec- terns am: 
also their exclusion laakes the ariaature eguaticns ilge ;rair 
which is consistent with the steady-state reox ose. nation cf 
the network to vjhich the machine ariteture circuits are conn^^ 
cted. This results in the reduction of the orcer of the 
machine by two. 3y setting pf^ and pf^ egual to rero and 
w = in eqn. (D.5), the stator equations are reduced to 
a pair of algebraic equations given by : 


^d 



0 


f 

•H 

lJ 

■ 4 . 

0 

-I 


'“-.r 

V 

L c[ J 


0 

tiirital 




1 

Q 


f 
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Substituting the values of and from eqn. {D.2) Into 
eqn. (D.6) , we get : 


^d 

. 

-r L •* 

a q 




0 

-i*’ 

q 


MR. wiMiiM 

Id 

(4^.7) 

V 

_q J 






■-d' 

mm 

U 


-'d 



where 

4' “ '■d ■ ^ ^ i4 

Id = 

Iq = IT> Jiq !jcq 



S ubtrans lent 

Although and L” are not equal. In general, 
of duraiiiY coils can !nake tlm new values equal, relic 4nc i;*" 
general procedure outlined in lief. the effect cf tn:* 

subtransient saliency is taken into account and #*? get 


where 

... Jl 0 

and 


^r, + t 
q q 


I- = -nr 

rq** "Tq 


■q Lf 


A L" 

5 ^ 1 -rl 


{J.X 


The second tens in eqn. (D.C) represents subtransient 
saliency. L” and L" are the conventional subtransient 

Ci €| 

inductances. 


.acuations of Mechanical Motion 

The rotor mechanical motion 1$ represented by 


and 

where 

-■ fc 


m 




p6 =® w - Wq 





(m 


m I 


p -is the mechanical power input 



CD. 9) 
CD. 10) 


and 


- is the electrical powr output. 



s tate E quation s 

Unlike Hef. [7S]» the variables and I are dit'jctlv 

■%* ^ 

utilized as state variables. ThuSf replacing » -inc r** 
eqns. (D.3) and (D.4) by and respectively, tne dyr.nnic u 
equations of the generator awr given as : 


pid = Ci^d C3 * f Ciiij t c,Cfj (D.in 

pl^ = C2lq *■ ^^4 *k * ''12^q (0'12) 

p *f= ^ =7 * f * ^13^ * =ioyd 

P * k” Vq ■'■ ^8 * k * ''Ik^q (O.M) 

where 


M «»1 *•! 

=s ** ^ 

III _ ^ 

^2 kg 4 4 g k 4 ^ 

C3 = ^ [dgEhCLfi^ + 

i | | % : - H ■-■''*’^1 ^■: ^ 

C 4 = 1 § C^4:‘*g<‘-kg ‘-kC^3> " ‘*3’‘k‘‘-g * “3^ 

-1 



Cil 

'=8 = - Skil-g H. L^jd-^dj) 

iif 

C„ = u)^d, tt— 

9 0 1 .,.^j 

Rf 

'=10 = “o U-~ 

til 

^11 “ a7 ■" ‘‘^ ^ 

X 

■M 

^12 ^ “ ij " ^k/‘qg^ ^4%^“^kg"‘qk 

m 

^13 “ “ a^ %^^h^df "* ^fh^dh^ 

«i) 

=14 = - i; «lctV‘qk - 


^1 


1 1 X 

^ 4h 


^2 

“ ^k^g 

- L ^ 
^kg 





" ^dh^fh^ 

*^2 

■ »i4' 
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Linearized state-space equation of the generator car tc- 
obtained by linearizing eqns, {D.9) to CD. 14) as : 


P = [A' j -r b iiE. . t h C 6F 
^ rn-’ -m -®e fd '^%ig ' m 

where 

2^ = [Aw AI^ AI^ A¥^ 


ijA5) 




-Kp/M 

c 

0 

0 

0 

O 

1 



1 

0 

0 

0 

0 

0 



0 

0 

^1 

c 

c.. 

3 

c 


= i 

^ 0 

0 

0 

^2 

0 

-4 



0 

0 

s 

0 

Of 

0 



0 

0 

Q 

% 

0 

r 

^9 

■mmm 

b 

-me 

tl 

o 

o 

9 

0 


^10 

OJ* 


b 

~mg 

s [l/M 0 
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0 

0 

of 



and 

o o 

o o 

L_J 

II 

<=11 

c 

0 

^12 

^13 

0 

0 “]^ 
'"wJ 



Steady Stat e Re.prMeniiM51S 

In the steady state repiresentitlOR 

~ cosCw^t "tA I "" sifiCiSgt “* 6 I 





where is the current or voltage variable and S'ufcscrint 
O' ci 0 no’fc 0 s th© x^fGiTGiic© firsuB©* w C®® 4^ J I » 

speed and 6 is the relative angle. Park’s conponents f,. 

can be used to denote the magnitude of the ehasorc repre- 
senting f^. The phasor diagram is shown in Fi.j, i),3« The 
usual convention is to denote 6 as the angle between the 
quadrature axis and the reference frame. Thus, the chaser 
f_ can be written as 

fa = (fq-jfd^ ^ 

Hence, in steady-state the single phase equivalent circuit of 
the machine can be represented as shown in Fig. D.4, when all 
the phasors are referred to the reference frame. 

Expression for 


As the network is represented by its Jacobian matrix frc« 
which AS is directly available, Ai can be expressed In 
terms of AS^ and as follows. Generator power output is 
given by 



Linearizing Eqn. 


*4 - * I’W Art 


(o.ia) 




where 




r lao + 


Ido-2 W 


•(2^ - I ) 

^ do ^ dc' 

••{ 2 f i I 
' qo qo^ 


[W, 2 ] = ^dt“al 3 


-i 

qo do 
^do ^qo 


e C I C C C 


0 C C i 


= -(Iqo+ f Iqo> tlqo*2{ f-l)lqj-<ldo-21do>fldo- * 


Expression for Gener a tor Terainal V oltage 

The magnitude and phase angle of generator tenrinal vol- 
tage can be obtained from Fig. i>.5 as ; 


Vt = (vi . v2,I/2 

9 =6- tan“^ (^) 


(0.19) 


CD.20) 


Linearizing eqns. (D.19) and CD.2K))i we get 


^Vt = v^^] 




Ae = 


CD.22) 


v/here 


= [ iVj AY 1* - frf,3] Ai, ♦ tw,43 4. 


r c Ci-i n 

[' w = ^iT! ' „ 1 f-.il 


1 0 




FI6.D.5 PHASOR DIA«5RAM OF SVNCHRONOUS 
MACHINE . 



D.1.2 Lower Order Modei \ 

In lower order representation ©f a synchronous .i/c only ■ 
the effect of field flux decay is considered and the effect of 
ai^ortisseur windings is nffiectad* This xesolts in s three 



winding, ipodel of the synchronoys ©achine - d and q 
in the stator and the field winding f in the d~axis tf 
rotor. Thus, in Fig, D,1 will be replaced by ar.j 
will be no h, g and k windings in Fig. D.2. e.',f.s. {J.l) 
(D.2) will be reduced to 


V 


fd 


= Vf 



and 


Is ^ 



M 




11 




•S' 


CD. 26 ) 


^ 0# 2T J 


where 


M = O]^ 


Eliminating i^ from {D.26) and (0,27), we get 


P f , 


b'fd - ”o n *f -“o L7 “df ‘d 






CD.28) 


The stator equations for lower order machine model will take 


the form : 



(I). 29) 


where 


A 


i- 





T ^ 1 

d ~ L '’ r ”” 


and 


I„ = Si 
d q 

i = I 9 

L*. 

d 




Eliminating from eqns. {J,2B) and (D.X) we get 




YT” C“*Irf + It '^ffS 
^do ^ 


'd 


-fT" 




where 




R. 


Ci3.31| 


and 


A ^^df 


'fd 


*fd 


The dynamical equation of the machine will be 
given in eqn. (0.15) with state variables as 
Various matrices in eqn. (0.15) will take the fora : 

"-Kp/M 0 C 

1 0 0 

. 0 0 -1/TJ^ 


as 
and M 


d* 


[a;] 


-me 


[0 0 1/CLJ T^^)3 


mg 




[1/M C 0] 

0 0 {- fr* C 

d© 


1#^ ** ^ 

a_j)) 


0 



Expression for ; Matrices and in ec/., CiJ.lT) 
vail take the form : 


2{ §-l) i 


[«al 5 = — 


{|0 




[W^j] = [W^j] 


0 0 

0 c 


v/here 




t i 




m 




=4 = -28 (8 - 1 ) 


Expression for^v^ 

«*||| 


[Waj] = xA 


= *d 

Expression for 


in eqn. 

(D.23) 

(1-1)1 


0 _j 


(1-1 ) *1 


! 

c J 

1 


■c 0 O' 
.0 0 1 . 


Matrices [D j and [0^] in eqn. (0.24) will take the forr.i 

P 



D.1.3 Calculation of Operating Condition 

fhe values of Park's cwaponents of voltage and c-jzrz-rz 
at the given operating point of the syste© are calciilate’^ frw. 
the results of a load-flow analysis. Let P„ and w pe res- 
pectively the active and reactive power supplied Py the genera- 
tor at the tenainal bus and ^ 1.0 end 6 ^ represent 
the magnitude and phase angle of the terminal bus voitacje. The 
power factor angle and the stator current are, ressectivtiy. 

»o = 

^mo = Vf''to »o> 

From the phasor diagram shown in Fig, 0,5» the angle between 
the q-axis and the axis of bus voltage is given by 



Angle of lead of machine reference, q-axit, from the sy»t« 
reference is given as 


^do 

= ''to 

sin 

- 

V 

= v^. 

cps • 


qo 

to 

*^0 



"" ^mo 

sinip^ * 

f ) 

iqo 

a i 

mo 

COfC^0 

t ) 
o' 



and 



D .2 EXCITATION SYSTEM [9 ] 

D.2.1 Continuously Acting Regulator and Exciter (iSHE Tyce 1 

IEEE type 1 excitation systea is shown in Fig. D.6. 7 r 
linearized equations for the excitation tysttr; can iie written 

P -?e = * ^eu 

'^^fd = ^e ^<3.321 

where 



- (Kg + i 


1 



f 

^E 



0 

[A^] = 

0 


- VTa 



-K3(K£+ 

L 


% 

*E *S 

- VTj 


^ev = 


|> as [0 

-eu \ 

D,2.2 Single Tiae Constant Excitation Systeo {IcEE Type iSj 

A special cast of typ« 1 is » iysttro tiRplcying in 
excitation source froB temlntl veltsft with controlled 
rectifiers only. In general, 'th« constents for this type 
of system are such that Kg * 1» Tg » 0 end ^ I 9 ) * The 




> 



linearized equation for this excitation syster. tf 
as given in eqn. (D.32), where 




e - -1/T^ } » 


•Kg/j 


ty 


■ V'R 


Combined dynamical equation of machine and excita tl Ofl 
system can be obtained from eqns. (D.25I and (0.32) as s 


P = [A] X + [bJ ^ + b iu * (i3.33) 

p g m rmT 

where 

X 


t t-t 


m O'* 


[A] = 


A b C 

m ~me e 


b W , A 
-ev a6 e 


CBp] = 


B 


ap 


b W c 
ev a5 


^ = C^9 

b = [0 


W p. and W . are obtained froa matrices D„ and respectively 

0,0 ■ w 

of oc|ri« 3S ^ liii# in. 



APPE|.?OIX E 


PSS DESIGN BASED (M (XMtSICAL CC^IROL 


E.l DESIGN PR(X:EDURE 


The design of PSS using classical contrci theory given 
here is based on the work reported in [26], As iiscvised in 
Chapter 2, the design of PSS using the guide lines given In 
[26] requires the determination of and of the PSS 

transfer function 


PSS(s) 


KgCl + sT^)^ 

(1 + sT^/nJ^ 


The method suggested in [26] utilizes the following terof. 

Plant Transfer Function : GEP (s) 

This is defined as the transfer function from the stabi- 
lizer output to the component of electrical torque which can be 
modulated through excitation control. The stabilizer operates 
through GEP and therefore PSS should compensate for the gain 
and phase characteristic of GEP which is strongly influenced 
by the voltage regulator gain, generator power level and ac 
system strength. 

Compensated Phase lag • 

■ 

This represents the i^ase angle ©f the conoirte stasiilser 
path, i.e., phase angle ftm the Inptit signal [usuaily speed; 



to the electrical torque. This can be obtained as 


= yw) 

The procedure for the design of PSS consists Oi two ftasic 
steps [26]. 

(1) Adjustment of parameters n^. and Tj^ such that 

(a) Maximum frequency at vrtiich the coapensated phase i«g 
passes through should b® in the range of 3 to 

3.5 Hz. 

(b) Compensated phase leg at local node frequency should be 
between 0° and 45®, preferably near 20®. 

(2) Adjustment of gain Kg based upon the deteroination of gain 
which causes instability (say K^) in the stablliiex loop. 
Gain should be adjusted at 

V #*1 If* 

Kg « 3 K® 

E.2 EXPRESSIO!'’ FOR GEP{s) 

From eqn, (2.4), electrical torque in per unit is 

Tg * ^d V S 

Linearizing eqn. {E.2) and eubstitutifif th® values cf and 

Ai^ from eqn. (2.10), m get 

■ ..€1 



and 


= =‘dtKi{WSldoWl-Si<3>i,c3 

■% = *d [M^do-'W - ''4 5V 


Taking the Laplace transfora of the state equations given *.f; 
eqn. (2.11), we get : 


SA {0 s ^ 

+ a^2*« + a^gA 

\ ft .4) 

sA6 =Aw 


(£.5) 

s ^1^ ~ 832^®^ 

*33 *34 

(£.8) 


*43 ’*’ *44 '^^fd ^4 “ 

CE.7) 


and 


From eqns. (E,3) to (E.6), can be expressed as J 

a^p(s-a^^) + a34®42 ®34^4 . t 

= C 7^^— TTTirTT^ ^Ti''*3T>TJ%Fi^34S43^ 

(8-333) (8-8441-834843) 33 44 J 4 4 J 


(£. 8 ) 


Substituting the value of frcxs (1.8) into @qn» CE.3)t 
we get 

At^ = Gj(s)A 6 + GgCs) « 

33^ 


; =^' \3% 

e 1 
where 


r ro= K ^ 

'i 6 C 8-833) C *-*44) " *34*43 


62(5) ^ 'GEp(s) * *34*ir 


(2.10) 



Block diagram of the system giving relationship of clc-ctr-cil 
torque with speed and angular dispiacrr.ent is giver, in rig. E 
E.3 DESIGN OF PSS 

Speed signal is utilized for the ?Sa design hy cia-isic 3 l 
control. As recommended in [26], P5S utilizing speed signal 
should be designed for strong ac system and at full load 
conditions. Thus, the value of considered is equal to 
0.2 p.u. for the system given in Chapter 2, Fig. £.2 shows 
the phase versus frequency plot of phase angle of GEP and 
compensated phase given by eqn. (E.l), for different viiues 
of the parameters n^ and T^. It is found that curve 'a* of 
Fig, E.2 with n^ = 3 and = 0.5 sec. fulfills the design 
criteria given in Sec, E.l. 

Fig. E.3 shows the loci of closed-loop eigenvalues for 
variation in the PSS gain K,. it is found that the value of 
K_ which causes instability of the systea is equal to O.l&S. 

O' 

Thus, as per the criteria given in # gain of C*C5 

c=-4 K* ) is selected for the PSS. 

o S 










^P^DIX F 


AN.^YSIS OF RESULT OF CJIAPTSR 4 


The state-space equation (4,26) for the t\io ni.'^chin^ 
system can he expressed as 


hi hi 


hi hi 



"" 





Defining new state variables 



^1 

A 

= Ll 

- 





A 

= 

+ *2 



we get the 

transformaticwi as 



r - ,-Ti 


■ ii- ,n 


r 


^1 


1 1 


% 



1 

= 1 

j-I £ 


h 

lUw w 4mm 


0 Isg 


u 


I 


(F.i) 


(F.2) 


where I is an identity laatrix. aith the defined trinsforrsi- 
tion, eqn. (F.l) will take mn form 



For identical machines and with ec|uai loading i v® hs70 


A22 = All 

‘^21 "" ^12 


and thus eqn- CF.3) will be reduced to a pair of dscoipled 
state equations as 









"b 

0"“ 


4 

p 


= 1 




4- 






_A_ 




2/«i 


[? 

b 

^ 1 

;1 


(F.4) 


where u£ = U2^"'U2 and u^ — For plant signal®!' that is» 
average speed, average power and bus frequency and with identi- 
cal PSS on each machine U2 will be equal to u^ ind thui, 


u^ = 0 * 


(F*5} 


and eqn. {F.4) can be expressed as 


^1 

Z2 




0 iAjj+. 



*^1 

32. 


4 


(F.&) 


which proves that the systTO i> not eawletaly cootrolUoU. 



APPBiDIX G 


ME WORK REPRESeHMICH 

The network can be represented by its Jacobian matrix in 
the polar form. For an *N* bus power systeia network, the 
Jacobian matrix is of (2H x 2H) diaensicm. By representing 
the network by its Jacobian the identity of all the buses 
is retained. The elements of the Jacobian matrix of the 
network can be obtained by one of the following ways : 

The Jacobian of the network, excluding the slack bus, it 
available from the load-flow analysis using Newton's saethod 
in polar form. The elements of Jacobian corresponding to the 
slack bus are obtained using the loirf-flow results as desert* 
bed in the alternative TOtbod given here* Thus, the Jacobiah 
available from the load-flow analysis can be conveniently 
utilized in obtaining the full Jacobian of the network* 

Alternatively, if the Jacobian is not available frees the 
load-flow analysis, the Jacobian for the cowplete network can be 
constructed using load-flow results as follows s 

Eqn. (5.1) can also be expressed as 



where 


= [ APj AQj]t 
AZ. = [ AS. iVj]* 

J 


11 


^±3 


IS a 2x2 matrix representing the diagonal block 
of [j] corresponding to ith bus in eqn. CS.i) ' 

is a 2x2 matrix representing the off-diagonal biccks 
of [j] corresponding to ith bus in ©qn. (5.1). 


[Jij] = 


[j. .] = 


where 



Vt ) 
*io' 


— ^^io “ %i 

V? ) 
xo' 

(Qio - 

^io ^jo %j 


%o 

__“Ho ^jo ^ij 


^io %J 


Hi 


Ci- = G. . cos(e^^ - 0.^) + 3^^ sinca^^ . ip 
°ij ~ %j ■* ®jo^ “ ®ij ” ®jo^ 


Y. . = G. - + IB. . -ijth element of bus adnsittance natrix. 
ij i3 ^3 

P- » Q- » and 9 .^ are respectively iictlve power, reactive 

XO ^ , ^X0 .X'C#' xP 

power, voltage magnitixie voltage angle at the ith du'‘/ 

obtained from load-flcw results* 



APPENDIX H 


JUSTIFICATION OF CMTERIA IBBD FOR CHOICE OF ?S5 

LCCATiai 


The linear mathensatical Hiodel of a ayltiRiachini power 
system used for the selection of PSS location can fct expressed 
as 

P .^ = [Aj ^ + [b] y iMA) 


Y = [Cj X 


CH.2) 


and the PSS structure can be expressed as 


% = ^i ^i 


(K.3) 


Derivation of these equations are given in Chapter 5. ‘♦hen 
the ith machine is chosen for PSS, eqns* CH.i) and CH»2) 
will take this form 


p X = [a] 21 


and 


y. s Ci X 


CH.4) 

|H.5) 


where is a row 


vector, satrix C in «cn» Ct.2l) 


Jus tiJ[icj_tion j>f_FirstJ^^ 

Consider the transforation 



23 : 


which reduces the system of equation (H.4) to 
p Z = [r^ A t3 Z + 

= [-A.] z + b|_ Uj (H.T} 

where -A. is a diagonal matrix (assiaiing distinct eigenvalues/ 
whose elements are the eigenvalues of the A matrix. The ro-'*- 
in eqn. (H.7) corresponding to the critical eigsnvaiu-j Ky 
can be expressed as 


n T. = A. i. + b • • u» 
^ 3 3 3 31 “i 


CH.8) 


where bt. is the jth element of b! . Assuming that the control 

law for u. is going to be selected in order to shift the eigen- 
^ • 

value A. to the left in con^lex plane, it can be observed that 
the control effort would be small if the raagnitudt of element 
b’.^ is large. It is to be noted that b^| is the scalar product 
of the eigenvector M- of the matrix A^ corresponding to 
and the vector b|, or 

=«j-V 

■TiiRt.i f1 nat ion of Second }Aeth^ .,lTL.§^.Sj^»Sj3 

****^'***^“^ ■ " . 

From eons. {H.3), (H.4) .n« <K.5). tl» olo«<i-lcop syst.-; 


can be represented as 

p X = (A + % % 





The value of the feedback gain, is so chossR -as tc shi''"*- 
the critical eigenvalue to the left in co.colox clans. Tni- 
is to be achieved with smallest possible value of giis. Tne 
closed-loop eigenvalue A. is given by 

. c 

~ *^i Higher order terms (H.IO) 


If higher order terms are neglected, the -^est control is deter- 
mined by maximizing the absolute value of 

* i 

eigenvalue sensitivity is given by 


s.. ^ 


31 


<iy 


(K-ll) 


where V. and VI . are respectively the eigenvectors of A and A 
3 3 

corresponding to A^. 



APPENDIX I 


13-4-IACHINE 71-BlB POWER SYSTEM 

The single line diagram of a 13-machiiii, 71-by$ ps.ver 
system is shown in Fig. i.i. Table I.l shows the Dus data 
and load-flow results. Generator, lint and shunt c;»! 3 acitor 
data are given respectively in Tables 1.2, 1.3 and 1.4 in per 
unit on 200 MVA base. 

Table I.l Bus Data and Load Flow Result 


Bus 

Generation 

Load 

Po’mr 

Bus 

Voltage 

NO m 

(MW) 

(MVAR) 





Apes. P (Otg » .1 

1 ' 

2 

3 

4 


6 


1 

630.667 

169.722 

0.0 

0.0 

1.03D 

0.0 

2 

0.0 

0.0 

0.0 

0.0 

1.057 

-5.931 

3 

506.0 

149.492 

0.0 ' 

0.0 

1.025 

-i .406 

4 

0.0 

0.0 

0.0 

0.0 

1.054 

— S .061 

5 

0.0 

0.0 

0.0 

0.0 

imi 

-5.064 

6 

98.0 

31.986 

0.0 

0.0 

1.025; 

0.397 

7 

0.0 

0.0 

12.8 

8.3 

1.045 

-5.723 

8 

297.0 

124 .227 

0.0 

; : 0*0 

1.025 

' -C.2w4 

9 

0.0 

0.0 

185.0 

130.0 

1.0433 

• -4.976 

10 

0.0 

0.0 

804) 

30.0 

1.026 

-7.413 



1 

2 

3 

11 

0.0 

0.0 

12 

0.0 

0.0 

13 

0.0 

0.0 

14 

0.0 

0.0 

15 

184.08 

114.0 

16 

0.0 

0.0 

17 

0.0 

O 

• 

o 

18 

0.0 

0.0 

19 

0.0 

0.0 

20 

0.0 

0.0 

21 

0.0 

0.0 

22 

0.0 

0.0 

23 

0.0 

o 

# 

O 

24 

0.0 

0.0 

25 

0.0 

o 

6 

26 

0.0 

O 

• 

o 

27 

304.0 

76.287 

28 

0.0 

0.0 

29 

261 .0 

70.506 

30 

0.0 

0.0 

31 

0.0 

0.0 

32 

0.0 

0.0 

33 

0 

0.0 


4 

5 

6 

155.0 

96.0 

1 .023 

0.0 

0.0 

1 .& L '3 

ICO.O 

62.0 

1 . C 07 

0.0 

0.0 

l , w 04 

0.0 

0.0 

i . C 05 

73.0 

45.5 

1.024 

36.0 

22.4 

I . CX )5 

16.0 

9.9 

1. C 09 

32.0 

19.8 

0.901 

27.0 

16.8 

1. CXJ 3 

32.0 

,19.8 

1 .016 

0.0 

0.0 

1. C 08 

75.0 

46.6 

1.036 

0.0 

C.O 

C .985 

133.0 

82.5 

1.039 

0.0 

C .0 

1.016 

0.0 

0.0 

1.025 

X .0 

2 C .0 

1.055 

0.0 

0.0 

1 .025 

1^.0, 

74,5 

1. C 57 

160.0 

99.4 

1.014 

0.0 

C .0 . 

1.025 

0.0 

0.0 . 

■ 0.993 




7 

- a .- t ; 

-13. S3-) 
-11.934 
-9.413 
-16.682 
-14.460 

-n .881 

-19.786 
-21.838 
-21.376 
-17.625 
-19.315 
—16.193 
-17.441 
-14.382 
-5.453 
—12 • 342 
■ 41.519 
-19. 9-’'' 
- 22 . 31 '' 
“ 1 "; .'26 
-25.515 



1 

2 

3 

34 

0.0 

0.0 

35 

0.0 

0.0 

36 

0.0 

0.0 

37 

0.0 

0.0 

38 

0.0 

0.0 

39 

25.0 

30.383 

40 

0.0 

0.0 

41 

0.0 

0.0 

42 

0.0 

0.0 

43 

0.0 

0.0 

44 

180.0 

55.037 

45 

0.0 

0.0 

46 

0.0 

0.0 

47 

0.0 

0.0 

48 

341 .0 

256.0 

49 

0.0 

0.0 

50 

0.0 

0.0 

51 

0.0 

0.0 

52 

0.0 

0.0 

53 

0.0 

0.0 

54 

0.0 

0.0 




4 

5 

c 

112.0 

69.5 

1 . 3 C 6 

0.0 

0.0 

0.966 

50.0 

32.0 

C .994 

147.0 

92.0 

0,967 

93.5 

88.0 

C .972 

0.0 

0.0 

I . C ‘25 

0.0 

0.0 

1.018 

255*0 

. 123.0 

1.007 

0.0 

0,0 

C .979 

O 

• 

o 

0.0 

C .946 

0.0 

0.0 

1 .025 

0.0 

0,0 

1 .043 

78.0 

■ 33.6 

1.020 

234,0 

145.0 

0.989 

0 .0 

0 , 0 . 

1.005 

295,0 

183.0 

0.997 

40,0 

24.6 

0.977 

227.0 

142,0 

1 , 0)5 

0.0 

0.0 

C .957 

0.0 

0,0 

0.972 

1 C ©.0 

■ ■ 68,0 



r- , 


- 2 - 9 . 6 :-' 
- 21 . 3^1 
*V ,CK 
*- 37 . 9 a 2 
•■37.63£ 
- 3^.166 
-36.913 
-3^ .773 
-31 .484 
-27.335 
-74 .448 
-30.514 
-32.678 
-29 .053 
-25.731 
-30. MC 
-32 .056 
-35. X9 
-26.562 
- 32.345 
—35 • 31 .5 





1 

2 

3 

4 

55 

0.0 

0.0 

25.5 

56 

0.0 

0.0 

0,0 

57 

0.0 

0.0 

55.6 

58 

0.0 

0.0 

42.0 

59 

0.0 

0.0 

57.0 

60 

0.0 

0 .0 

0.0 

61 

. 0.0 

0.0 

0.0 

62 

0.0 

0.0 

40.0 

63 

0.0 

0.0 

33.2 

64 

300.0 

72.882 

0.0 

65 

0.0 

0.0 

0.0 

66 

96.0 

25.636 

0.0 

67 

0.0 

0.0 

14.0 

68 

89.0 

26-689 

0.0 

69 

0.0 

0.0 

0.0 

70 

0,0 

0.0 

11.4 

71 

0.0 

0.0 

0.0 


H tf *S( »i*l' (. 


5 

6 


48.0 

0.985 


0.0 

1.013 

- 29.253 

35,6 

1.016 

- 29.957 

27.0 

1 .013 


27.4 

1.013 

-30.715 

0.0 

1.008 

-29.974 

0.0 

1.019 

-27 . C 74 

27.0 

1 .044 

-28.940 

20.6 

1 .042 

-27.695 

0.0 

1.025 

-19.401 

0.0 

1.057 

-24.662 

0.0 

1. C 25 

-20.610 

6.5 

1.055 

-25.641 

0.0 

1.025 

-19.029 

0.0 

1.050 

-25.363 

7.0 

0.998 

-34.79 C 

0.0 

0.927 

-35.329 



Table 1,2 Generator Oata 


Gen. 

No. 

Bus 

No. 

Rated 

MVA 

Inertia 
Constant 
H secs. 

Open Circuit 
Field Time 
Constant 

'*do 

Direct jir«ct 

Axis Axis 

Reactance Trar.sie:.*. 
X. c.u. Reattar 

1 

1 

1175.0 

12.9CX3 

7.C 

C.36 

c.c48e 

2 

3 

665.0 

9. OX) 

5,5 

0.535 

O.C69 

3 

6 

110.0 

1.923 

6 itO 

1.638 

C.5278 

4 

8 

335.0 

6.648 

6.C 

0.567 

C.2C1 

5 

15 

275.0 

2.590 

6.0 

0,6 

0.244 

6 

27 

400.0 

2,550 

6.0 

0.9 

0.1455 

7 

29 

300.0 

2,700 

6.0 

1,578 

0.1355 

8 

39 

51.75 

1.035 

6.0 

2.705 

l.C«2 

9 

44 

220.0 

5.430 

6.0 

0.9118 

0.227 

10 

48 

500.0 

7.800 

6 *0 

0.74 

0.«»95 

11 

64 

376.0 

7.140 

6,0 

C.532 

0.13 

12 

66 

127.8 

2.176 

6.0 

s, ** 

0.486 

13 

68 

103.5 

1.473 

6.0 

2.28 

C.599 



Table 1.3 Line and Transfonner Data 


Line 

No, 

From 

Bus 

To 

Bus 

Series Impedance 

R X 

1 ¥, 

^ charge 


2 

3 

4 

5 ■ 



1 

9 

8 

O.C 

0.0570 


2 

9 

7 

0.0320 

0.0780 

0.0090 

3 

9 

5 

0.0660 

0.1600 

C.C047 

4 

9 

10 

0.0520 

0.1270 

0.0140 

5 

10 

11 

0.0660 

0.1610 

0.0180 

6 

7 

10 

0.0270 

0.07CX) 

C.007C 

7 

12 

11 

0.0 

0.0530 


8 

11 

13 

0.0600 

0.1480 

0,0360 

9 

14 

13 

0.0 

0.03C0 


10 

13 

16 

0.0970 

0.2380 

0.27C 

11 

17 

15 

0,0 

0,0920 


12 

7 

6 

0.0 

0.2220 


13 

6 

4 

0.0 

0,0800 


14 

4 

3 

0.0 

0,0330 


15 

4 

5 

0,0 

C.16<X) ' 


16 

4 

12 

0.0160 

0.0790 

0.071C 

17 

12 

14 

0.0160 

0.0790 

0.0710 

18 

17 

16 

0.0 

O.OSCX) 


19 

2 

4 

0.0 

o.c^ao 


20 

4 

26 

0.0190 

0.0950 

0.1930 

21 

2 

1 

0.0 

0,0340 


22 

31 

26 

0.0340 

0.1670 

0.15CX3 


Turns 
Hat 1 C 

i,C5 


0,95 

i.oo 

1X5 

1X5 

l,OC 

1,05 

l.CX) 


0. 95 

1. CC- 

1X5 





2 .!“ 


1 

2 

3 

4 

5 

6 


23 

26 

25 

0.0 

0.0800 


C.95 

24 

25 

23 

0.2040 

0.5:^ 

0.013C 


25 

22 

23 

0.0 

0.0800 


C.95 

26 

24 

22 

0.0 

0.0840 


C.95 

27 

22 

17 

0.0480 

0,2500 

0.C505 


28 

2 

24 

0.0100 

0.020C 

0.3353 


29 

23 

21 

C.0366 

C.1412 

0.0140 


30 

21 

20 

0.0720 

0.1860 

O.CXJSO 


31 

20 

19 

0. 1460 

0.3740 

C.OlOO 


32 

19 

18 

0.0590 

0.1500 

C.004C= 


33 

18 

16 

0.0300 

0.0755 

0.G08C 


34 

•28 

27 

0.0 

0.0810 


1,05 

35 

30 

29 

0.0 

0.0610 


1,05 

36 

32 

31 

0.0 

0,09X 


0.95 

37 

31 

30 

0.0 

0.0800 


0.95 

38 

28 

32 

0.005 

0,0510 

0.6706 


39 

31 

33 

0.0130 

0.0640 

C.C58C 


40 

31 

47 

0.0100 

0.07» 

0.1770 


41 

2 

32 

0.0158 

0,1570 

0.5100 


42 

33 

34 

0.0 

0.08C3C 


C.95 

43 

35 

33 

0.0 

G.m4C 


C.95 

44 

35 

24 

0.0062 

0.0612 

0.2012 


45 

34 

36 

0.0790 

Q.WIC 

O.C22C 


46 

36 

37 

0.1690 

0.4310 

■ O.CllO 






1 

2 

3 

4 

5 

6 

7 

47 

37 

38 

0.0840 

c.ieeo 

G.C2iC 


48 

40 

39 

0.0 

C,38CC 


I.C5 

49 

40 

38 

0.0890 

C.217C 

C.C25C 


50 

38 

41 

G.109C 

0.1 96C 

C.C22: 


51 

41 

51 

0.235C 

0.6CCC 

o.ciec 


52 

42 

41 

0.0 

0.C53C 


C.96 

53 

43 

42 

0.0 

0.0S4C 


^ 0.95 

54 

47 

49 

0.0210 

C.103C 

C .09X 


55 

49 

48 

0.0 

0.0460 


l.C» 

56 

49 

50 

0.0170 

0.0840 

C.C760 


57 

49 

42 

0.0370 

C.1950 

C.039C , 


58 

50 

51 

0.0 . . 

C.Q53Ci, 


C.^ 

59 

52 

50 

0.0 

C.C840 : 


■ C.95 

60 

50 

55 

O.C290 

0.1520 

0.0X0 


61 

50 

53 

O.ClOO 

0*0520 

0.039C 


62 

53 

54 

0.0 

0.C3CO 


0.95 

63 

51 

54 

0.022C 

C.C54C 

0.CC6C 


64 

55 

56 

C.016C 

C-0S50, 

c.onc 


65 

56 

57 

0.0 

O.CSCC 


l.CC 

66 

57 

59 

0.C28C 

0.072C 



6.7 

59 

58 

C.0430 

G.1240 

v^' * V 1. 


68 

60 

59 

'' 0*C : , 

O.OSCO 


A * '*» V 

69 

53 

60 ' " ^ 

, '' 0.0360 

. 0.1840 

m *5 ^ 

1 ** ^ 


70 

45 

44 



0.12CC 


1*05 




1 

2 


71 

45 

46 

72 

46 

41 

73 

46 

59 

74 

60 

61 

75 

61 

62 

76 

58 

62 

77 

62 

63 

78 

69 

68 

79 

69 

61 

80 

67 

66 

81 

65 

64 

82 

65 

56 

83 

65 

61 

84 

65 

67 

85 

67 

63 

86 

61 

42 

87 

57 

67 

88 

45 

70 

89 

70 

38 

90 

33 

71 

91 

71 

37 

92 

45 

41 

93 

35 

43 

94 

32 

52 


4 5 


0.0370 

0.09CO 

0,0830 

0.1540 

0.1070 

0.1970 

0.0160 

0.0830 

0.0 

0.0800 

0.0420 

C.1C8C 

0,0350 

0.0890 

O 

• 

o 

C.222D 

C.C230 

C.116C 

0.0 

0.1880 

0.0 

0.C^30 

0.0280 

0.1440 

0.0230 

0.1140 

0.0240 

0.0603 

0.0390 

c.ow 

o,02:» 

0.2293 

0.0550 

0.2910 

0.1840 

0,4680 

0.1650 

0.4220 

0.0570 

0.2960 

0.0 

o.woo 

o.i5:m) 

C.^0 

0.0131 

ca^ 

0-0164 

0,1632 




5 

O.OiCC 

C.C17C 

C.C2iO 

C.016C 

0. 95 

C.002C 

CXQ9C 

C.i04C 

1. C» 
1.05 

0.0290 

0.024C 

0.0950 

0.0100 

O.C3695 

0,0070 

0.0120 

0.0110 

0.05» 

0.95 


G.OICC 

0.4293 

C.5350 




Table 1.4 Shunt Capacitor Data 


S.No. 

Bus No. 

Shunt Load 

■ Afcittance 

1 

2 

0.0 

#4275 

2 

13 

0.0 

0,1600 

3 

20 

0.0 

C.CBCX5 

4 

24 

0,0 

-0.2700 

5 

28 

0.0 

-0.3375 

6 

31 

0.0 

o.«» 

7 

32 

0.0 

-0.8700 

8 

34 

0.0 

0.22SO 

9 

35 

0.0 

-0.3220 

10 

36 

0.0 

0.1000 

11 

37 

0.0 

C.3500 

12 

38 

0.0 

0'.2000 

13 

41 

0,0 

0,2000 

14 

43 

0,0 

-0.2110' 

15 

, 46 

0.0 

0.1000 

16 

47 

0.0 

0.^»D0 

17 

50 

0,0 

O.ICXXJ' 

18 

51 

0*0 

0.0 

C.175C 

-C.27Ce: 

19 

52 


20 

54 

0.0 

o.isa’ 

21 

57 

, 0.0 

0 * 

22 

,59 . ' 

, 0*0 , 

V « V ^ ’rfV' 

23 


; 0.0''''; 
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APPB\JDIX J 


3-MCHINE 9*-BUS Pa'iER SYSTEM 

The single-line diagram of a three machins Hint hy.f ptwtr 
system is shown in Fig, J,l. Line data are given in TaWt J,l. 
Generator and Voltage Regulator data are given in Table J.2. 
Load-flow results with line flows are marked on the diagrar. 
given in Fig, J.2, Data given in Table J,1 and Table J*2 are 
in per unit on 100 MVA base. 


Table J.l Line and Transformer Data for the S-Machin# 
System 


Line 

From 

To 

Series Impedance 

Half Line 

lUJWIt 

No. 

Bus 

Bus 

R 

j 

Charging 

Adteittanct 

Rati® 

1 

1 

4 

0.0 

0.0576 


1.0 

2 

2 

7 

0.0 

0.0625 

■m 

l.C 

3 

O 

9 

0.0 

C.0586 

- 

1.0' 

4 

4 

5 

0.0100 

0.CB50 

C.CiSO 

« 

5 

5 

7 

0.0320 

C.1610 

’ 0.150)' 

- 

6 

6 

9 

0.0390 

0,170) 

0.1790 

- 

7 

7 

8 

0.0085 

0.0720 

0.0745 


8 

8 

9 

0.0119 

O.IOCB 

'0.1045 

im 

9 

4 

6 

0.0170 

0.09^^ 

c.oim 



2 :- 


Table J.2 Generator Data for the 3-Machine oysterfi 


w m m * t'v 


Generator 

1 

2 

3 

Rated MVA 

245.5 

192.0 

1* • -a* V 

120.0 

KV 

16.5 

18.0 

13.8 

Power Factor 

1.0 

0.85 

• D , 

Type 

Hydro 

Steam 

Steas 

Speed (rpm) 

180 

36C0 

3600 

(p.uO 

0.1460 

0.8958 

1.3125 

(p.u.) 

0.0608 

0.1198 

0.1813 

Xq (p.u.) 

0.0969 

0.8645 

1.2578 

X’ (p.u.) 

0.0969 

0.1969 

. 0.»30 

4 

Xj^ leakage (p.u.) 

0.0336 

0.C521 

0.0742 

^do 

8.96 

6.0C 

5.89 

T’ sec 

0 

0.535 

0-600 

qo 

H (MWS/lOO MVA) sec 

23.64 

6.40 

3.01 

VR Gain 

50.0 

50.0 

50.0 

Trr, sec 

0.020 

0.020 

C.020 
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